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ABSTRACT 



Motivated by a lack of explosive test data on nu cl ear- 
submarines , the Navy has sought other means to quality in 
stalled equipment in submarine shock environments. The 
currently used method for non -shock testable items is the 
Dynamic Design Analysis Method (DDAM) developed by the 
Naval Research Laboratory in the early 1969's. With the 
adv en t of 1 arge-sc ale c omp u ting p owe r , n ewe r nume r i c a 1 
methods have become available to predict equipment re- 
sp on se s . Th i s i n vest i gat i on i s a c omp ar a t i v e s t u dy of 
DDAM and ELSHOK ; a new generation numerical shock res p on s- 
code. The limitations and strong points of both methods 



are examined using illustrative 



ex amp 1 es . 
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I . INTRODUCTION 



A. BACKGROUND 

Modern combat vessels must be designed with the 
capability to survive moderately severe shock loadings 
i nduced by the u nderwa ter e x p 1 os i on of c on v e n t i on a 1 or 
nuclear weapons. In the scope of this investigation, 
survival includes the mission survival of the platform. It 
must maintain its ability to utilize all machinery and 
we ap on s sys t ems to c ar r y ou t its p r i mar y m i ss i on f o 1 1 ow i n g 
an underwater explosive attack. Although the vessel may 
withstand considerable hull damage and maintain its 
structural integrity, failure of internal equipment can 
render it useless as a weapon and thereby el iminate its 
value in time of crisis. This investigation examines the 
shock loading of internal equipment in submarines and some 
of the methods available to analyse shock response. 

In view of the fact that a submarine can not be 
expected to survive a direct hit by any modern, 
h i gh- i n tens i ty weapon, many attempts have been made to 
specify a level of shock loading for the purpose of design 
e v a 1 uat i on . The current sp e c i f i c a t i on s f or bu i 1 d i n g 
submarines contain the shock requirements to be met by the 



builders end vendors of installed equipment. Generally, 
all equipment is to p ass a series of sh oc k tests, i n t h e 
installed configuration as outlined by MIL-S-901C [Ref.l], 
where shock testing is practical. This document specifies 
standard explosive and mechanical test fixtures and 
procedures for conducting the tests. The suitability of an 
equipment design or installation is evaluated according to 
its ability to function as intended during and after each 
shock impulse. Equipment testing is very expensive and 
requires a great deal of preparation. In the case of 
equipment design, it is often not practical to evaluate an 
installation prior to the final design. In some instances, 
the direct testing of installed components may not be 
possible from the aspect of prohibitive size or weight. In 
these cases, a design analysis is specified. 

Th e de s i gn an a 1 y s i s f or sh oc k r e sp on se was f i r s t 
proposed in the United States during World War II. Hull 
damage reports indicated the degree of sucep t i b i 1 i ty of 
installed equipment to shock loadings from bombs and depth 
charges. Based on limited testing and observations, a 
shock design factor was specified. A series of curves were 
used to depict the variation of shock design factor with 
equipment weight to yield an equivalent static design 
force. This force was used to specify mounting hardware 
and main structural members for the equipment installation 



and varied with motion input direction. The method proved 
to be simple in appl i cat ion but reflected a lack of real ism 
with experimental tests. 

In the late 1950's, nuclear power plants were 
incorporated into submarine designs to afford therri greater- 
independence from the frequent surfacing requ i remen ts 
inherent in conventional submarine operations. Nuclear 
p owe r plant tech n o 1 o gy was new an d t h u s h ad n o p r e v i ou s 
history of shock resistance. The acquisition cost of 
n ucl e ar c omp on ents is hi gh an d the n ew dan ger s assoc i ate d 
with their damage from underwater explosion sparked a 
renewed interest in underwater explosion testing and design 
analysis. From the late 1 950 ' s through the middle 1960's, 
mos t of the currentl y u se d s t an dar ds f or su bmar i re- 
installed equipment shock design were adopted. 

The Dyn ami i c De s i gn An a 1 y s i s Me t h od < DDAM ) was p r op ose d 
by the Naval Research Laboratory and accepted by the Navy 
as a design evaluation requirement in 1963 [Ref 23. The 
details of this specification will be discussed more fully 
later; h owe v e r , the basic principle is stated here. DDAM i s 
a simplified miodal analysis method which utilizes shock 
i nputs wh i c h we r e emp i r i c a 1 1 y de r i v e d f r cm u n de r wa ter 
explosion tests of realistic ship and submarine instal- 
lations. It is assumed that the equipment and its foun 
dation together make up a system which responds as a linear- 
elastic structure to the input which is described by a 
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des i qn shock spectrum. The successful util ization of D DAM 
to evaluate equipment response to a design shock input is 
dependent on the ah i 1 i ty of the design shock spectrum to 
reflect accurately the structural environment existing in 
given installation. This will be influenced by the size 
and weight of the submarine test subject, the structural 
ays t enri emp 1 oy e d in the c on struct i on of the su bmar ine, an d 
the interaction between the equipment, foundation and hull 
structure used to derive the design shock spectrum. 

B. THE STATE OF CURRENT SHOCK DESIGN ANALYSIS 

Since its ad op t i on , DDAM has bee n u se d w i t h ou t 
modification to qualify ship and submarine equipment 
installation designs. The original design shock spectrum 
[Ref. 23 remains intact without revision. Shortly after 
t. h e adop t i on of DDAM as a de s i gn r e q u i r erne n t , u n de rwa ter 
explosive testing of nuclear submarines was banned due to 
the inherent r i sk s i n v o 1 v e d . This ban was r e 1 ax e d in 1 983 
for the low level explosive test of a 3SN 693 class 
su bmar ine. Me c h an i c a 1 and e x p 1 os i v e tests of equ i pme n t 
installations on shock simulation platforms are the only 
regular check of designs nrieeting the standards imposed by 
DDAM. A danger exists with this kind of verification in 
that the tests miay be designed to accurately reproduce the 
input motions provided for in the design shock spectrum 



utilized with DDAM . 



In recent years the increasing availability of large 



computers and the rapid development of numerical methods 
have provided the engineer with new tools to use in the 
analysis of submarine shock response problems. 

F i n i t e -E 1 erne n t/F inite -Difference methods all ow f or the 
analysis of structure and fluid responses and the work of 
Geers [Ref. 33 provides a means to incorporate f 1 u id- 
structure interaction effects. Several computer codes have 
been developed to analyze installed equipment response to 
shock waves based on these formulations. One of them is the 
El as t i c Sho ck (ELSHOK) code developed by LJeidlinger 
Associates under Defense Nuclear Agency and Office of Naval 
Research funding [Ref. 43. The motivation for this work has 
been to investigate modern submarine underwater explosive 
shock response in conjunction with a testing program using 
small to large scale models and shaped, explosive, tapered 
charges. Throughout its development, ELSHOK has been 
validated using these highly controlled tests. In the low 
level explosive test of the SSN 688 class submarine, ELSHOK 
was used to predict the level of equipment responses prior 
to the tests and to increase the level of confidence 
associated with the tests. This code has now been made 
available to the Naval Postgraduate School by the Defense 
Nuclear Agency along with support from UJeidlinger Associates. 
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c. 



PURPOSE FOR THIS: INVESTIGATION 



DDAM i s a de s i gn tool • Its tunc t i on is to prov i de a. 
numerical method by which the engineer can check his design 
as to its adequacy for installation in a submarine. It is 
intended to be conservative and s t r a i gh t f or war d t o ap p 1 y . 
When DDAM was first proposed [Ref. 4], its authors warned 
that the design shock factors specified [Ref. 2] were not 
absolutely determined and should be reviewed as explosive 
testing progressed and data was accumulated. Although 
testing procedures for installed equipment components have 
become very sophisticated, these simulators may not accur- 
ately reflect the response of installations in modern day 
submar ines. W i t h ou t the ab i 1 i t y t o c ar r y ou t h i gh - 
intensity shock tests on real submarines, verification of 
the design shock spectrum used in DDAM can not be 
accomplished as intended. The evolution of numerical codes 
capable of simulating the response of installed equipment 
to a specified shock loading makes available a means to 
check the applicability of the originally proposed design 
shock factors provided in DDAM to modern submarines. There 
is no flexibility built into DDAM to all ow de s i gn f or o t h e r 
than the originally prescribed shock input magnitude. DDAM 
makes no separate account for equipment/structure inter- 



action effects. 



i n v est i gat i on is t o e x am i ne the 



Th e p u r p ose of this 
response of several equipment models installed in a present 
day scale, 6?d 0 long t on ( LT > su bmar i n e u s i n g DDAM an d 
ELSHOK . The e q u i pme n t mode 1 s r an ge i n we i g h t f r orri 
1,000-20,000 pounds. Simple models are used to make 
c or r e 1 a t i on s p oss i b 1 e be twee n the two me t h ods and wor s t 
case results are obtained. Additionally, equ i pmen t/hu 1 1 
interactions are investigated tor their possible influence 
on overall equipment response. 
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DYNAMIC DESIGN ANALYSIS METHOD 



I I . 



A. DESCRIPTION OF DDAM 

As. men t i on e d in the i n t r odu c t i on , DDAM i s a s i mp 1 i f i e d 
moda 1 an a 1 ys i s me t h od wh i c h utilizes sh oc k i n p u t s wh i c h 
i/i e r e emp i r i c a 1 1 y de r i v e d f r om u n derwa ter e x p 1 os i on test s of 
realistic ship and submarine installations. Although the 
analysis, is identical -for both types. of vessels, all refer- 
ences in this document will be to submarines. In order to 
utilize DDAM to evaluate a design, several assumptions must 
be sa t i sf i ed . 

1. ) The equipment and foundation make up a linearly 

e 1 ast i c system . 

2. ) The structure can be represented reasonably 

well using a lumped parameter model . 

3. ) The structure is hard mounted to the rest of 

the submarine. No bottoming, or damped mounts 

are all owed . 

In Navy shock requ i r emen ts , equipment shock classi- 
fications are based on a graded system. These grades run 
from A for mission critical itemis, to C for items that 
could indirectly affect the ability of the vessel to 
function as intended. For instance, a locker which becomes 
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adrift and could injure personnel or surrounding equipment 
could be grade C. The grade of the i t em is specified by 
contract and determines the degree of shock qualification 
it must undergo to be considered safe for installation. In 
this study, all equipment models were considered grade A. 

The ideal way to compute the elastic shock response of 
an installed equipment is to consider the entire structure 
•'submarine and equipment) as an elastic system. The normal 
modes and natural frequencies of this system can then be 
determined and the response of each normal mode to the 
water applied pressure loadings computed. The modal 
responses can then be superimposed to get the resultant 
equipment response. In an everyday design method such as 
DDAM, this approach is impractical . To simpl i f y the above 
procedure, DDAM utilizes the concept of a shock spectrum. 
Any point on the complete structure can be designated as a 
reference point which can then be considered to be a fixed 
base for the equipment on one side of it. The dynamic 
response of this portion can then be computed by 
superposition of normal mode responses to the motion of the 
fixed base. In DDAM, this motion is specified as a step 
input to each normal mode and is calculated by empirical 
equations based on underwater explosion tests of models and 
submarines [Ref. 53. A detailed treatment of the 
mathematical basis for DDAM is given by Butt [Ref. 63. 
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In order to use DDAM, a lumped parameter model of the 
equipment must be formulated to obtain the matrix 
express i on : 

CM] 66 + CK3 OO = 03 > <1) 

which has the associated eigenvalue problem 

CM] CX ia3 [ CO 2 3 = C K 3 CXj*3 < 2 ) 

where C X j a ] is the matrix of eigenvectors and CtO £ 3 is the 
diagonal matrix of natural frequencies. The advantage of 
this procedure is that it allows uncoupling of the 
equations of motion and subsequent solution of the dynamic 
problem in terms of its component modes. Careful 
formulation of the equations of motion to capture all the 
dominant characteristics of the equipment leads to good 
results with only a small number of modes. When DDAM was 
developed, it was intended to be used. pr i mar i 1 y in hand 
computations. The major difficulty encountered in the use 
of DDAM was the so 1 u t i on of the e i gen v a 1 u e p r ob 1 erri wh i c h 
can readily be handled today with small computers. The 
normal mode theory upon which DDAM is based is well 
established in shock and v i br a t i on practice. Th e area o f 
concern is the adequacy of the adopted design inputs to 
reflect present day submarine building practices. 

B. DESIGN SHOCK INPUTS 

The specification of design shock inputs for use in 
DDAM calculations determines the usefulness of this method 



17 



in the e v a. 1 u a t i on of equ i pme n t i n s t a 1 1 a t i on d e ■=• i qn s . DDAM 
is a standard. The design input values must provide 
consistent fixed-based equipment excitations for the broad 
base of submarine platforms in use today to simulate a 
"standard" underwater explosion. To simulate this require- 
ment here, charge weights and standoffs were selected so the. 
each case corresponded to a constant value of energy flux 
density; i . e . , 

W/ R 2 = C o n s t a n t < 3 > 

In this equation, W is the equivalent weight of the charge 
in pounds of TNT and R is the standoff distance from the 
hull in feet. The submarine platforms in use today range in 
size from about 4 , 60 0 -13,700 LT . Each submarine hull is 
characterized by its own modal properties (mass distri- 
bution, natural frequencies, mode s hapes) . The s am e piece 
of equipment as is installed in a small submarine will 
demonstrate a different response when installed in the same 
configuration, in a large submarine. The effect of hull/ 
equipment interaction becomes an important consideration in 
large pieces of equipment tuned to one or more natural fre- 
quencies of the hull structure. The design shock input used 
with DDAM must incorporate these factors so that when it is 
said that a piece of equipment meets the specification for 
qualification, the qualification level is the same for all 



su bmar i ne cl asse s . 



In the p u b 1 i c a t i on s r e gar ding DDAM , little is said 
abou t t h e or i g i n s of the p ar t i c u 1 ar sh oc k i n p u t sp ectr um 
utilized to evaluate equipment response. It has been 
empirical 1 y de r i v e d t o p r ov i de v slues c on s i s t e n t w i t h the 
data upon which it is based. In the event that it is 
decided to design submarines to resist underwater explosion 
of a different shock in ten si ty than the one chosen for the 
present method, a major effort would be required to 
construct a data base upon which to base the new input 
v a 1 u e s . 

C. USING DDAM 

S i nee the dev e 1 o p me n t of DDAM , n ew t oo 1 s h av e replaced 
the slide rules of engineers. Among them are the readily 
available de sk t op c omp uters. In the c ou r se of t h i s wo r k , i 
was necessary to use DDAM to analyse simple structures of 
fewer than 50 degrees of freedom. To this end, a computer 
program has been written in the BASIC computer language to 
carry out the required computations. A listing and short 
users' manual can be found in Appendix A. The program was 
verified by c omp ar i n g results w i t h h an d c a 1 c u 1 a t i on s an d 
published sample problems. 
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III. E.-PLUSH-'E SHOCK RESPONSE ANALYSIS USING ELSHOK 



A . GENERAL PR I NC I PLES OF OPERAT I ON 

The ELSHOK computer code [Ref. 4] calculates the 
transient response of a submerged, ring-stiffened shell of 
revolution, with or without internal structure, to an 
u nde rwa ter sh oc k wav e eman a t i n g -f r om an e x p 1 os i v e sou rce 
placed at an arbitrary location away from the shell. The 
structure is considered to be 1 i nearly elastic, and the 
surrounding fluid is treated as an infinite acoustic 
medium. Modal structural analysis is used in all phases of* 
the calculations. Internal equipment response is treated 
by coupling the free-free modes of the empty ring-stiffened 
shell and the fixed-base modes of the internal equipment 
through use of dynamic boundary' conditions [Refs. 7,83. 

The structure -fluid interaction is approximated by means of 
the Dou b 1 y Asymp tot i c Ap p r ox i ma t i on < DAA ) due to Geer s 
[Ref. 33. The form of the DAA used is that obtained when 
the normal fluid displacement of the structure -flu id 
interface is expanded using surface expansion functions 
which are orthogonal over the wet surface of the submarine. 
At frequencies of zero and infinity, the pressure -ve 1 oc i ty 
relations are exact so that in transient analysis, the DAA 
yields exact results at early and late times, and by the 
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nature of its f ormu 1 a t i on p r oy i de s a srrioo t h t r an s i t i on 
be twe e n the se two limits. In effect, the DAA a 1 1 ows f or 
uncoupling of the fluid field problem from the structural 
field probl em . 

The structural problem solved by ELSHOK is separated 
into two parts. A modal subs t rue tur i ng procedure is used to 
solve the dynamic response problem for internal equipment. 
The advantage of this is to el iminate the need to handle the 
modes and natural frequencies of the combined structural 
problem, as well as the requirement for a combined system 
stiffness matrix. Interaction forces and moments, and 
compatibility of deformation at the she 1 1 -subs true tur e 
attachment points are used to solve for the dynamic response 
of the component parts. 

Referring to the ideal case of DDAM formulated without 
the use of a shock input value; ELSHOK is a numerical means 
to arrive at the input to the substructure without depend- 
ence on an explosive testing database and with the added 
advantage that interaction effects between the hull and 
substructure are taken into account. ELSHOK performs a 
transient analysis whereas DDAM utilizes a simplified 
"front end" to arrive at the maximum forces and deflections 
in a given response problem for a single magnitude of 
1 oad i ng . 
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B. 



ORGANIZATION AND IMPLEMENTATION OF ELSHOK 



The ELSHOK code is implemented as a series of seven 
programs. The major components are: 



1. B0S0R4 - structural analyzer for shell [Ref . ? ] 

2. ACESNID - virtual mass processor 



3 . 


P I FLASH 


- sh e 1 1 - +' 1 u i d pr ocf ssor 


4. 


SAP IV 


- structural analyzer for substructure 
CRe-f. 16] 


5 . 


PICRUST 


- substructure processor 


6 . 


USLOB 


- time i n t e gr a t i on p r oc e ssor 


7 . 


PUSLOB 


- plotting processor 



Figure 1 depicts the general relationship am on g the 
code modules. The numbers above correspond to the 
chronological execution order to be followed in a given 
analysis. In order to carry out an analysis using ELSHOK, 
a lumped parameter model of the submarine being analyzed 
must be formulated or available as well as the equipment 
model. More will be said about models later. The compu- 
tations are carried out in four phases. 

B . 1 . Phase I - Shell and Fluid Analysis 

B0SQR4 is a finite difference code used to 
capture the effects of the submarine hull structure on the 
overall response. Two models are actually used: a full 
model containing information about the major configuration 



at . 

2 uj t 

nr 
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Figure 1 - Organization of the ELSHOK Computer Code [Ref. 4] 



of the submarine and its general dynamic properties, and a 
compartment model representing the local i zed area of 
interest on the hull. The -full model is used to capture 
the gross effects of the response such as rigid body 
translation and whipping ( N = 1 modes), and rolling and 
twisting <N = S torsional modes) of the entire structure. 
The compartment model is used to capture local details of 
the response in the area of interest <N = 0 breathing and 
N l 2 modes). N is the integral number of c i r cumf eren t i al 
waves or harmonics in the circumferential distributions of 
the surface expansion functions. For each N, a separate 
B0S0R4 calculation is required. The B0S0R4 code provides 
the in -vacuo free -free modes and natural frequencies of the 
submar i ne hull. 

The second part of phase I is the use of ACESNID 
to compute the virtual mass array, which is used in the 
late time contribution of the DAA . The virtual mass array 
is determined from the solution, based on simple sources, 
of a steady-state low frequency problem in which normal 
displacements corresponding to the surface expansion 
functions are applied to the surface of a cavity of 



revol u t i on i n an 


i n f i n i t e 


+ 1 


u i d 


h av i n g 


the 


same size 


an d 


shape as the wet 


surface 


o+ 


the 


submar 


i ne . 


Only one 


run 



this module is required for all values of N considered. 



The -final step in phase I is the execution o-f 
PI FLASH to compile the data from all B0SGR4 runs and 
ACESNID into a she 1 1 -and-f 1 u i d tile tor later use. The 
analysis done tor the work in this thesis considered only 
the hull re sp on se o-f on e su bmar i n e so on 1 y on e 
she 1 1 -and-f 1 u i d tile was needed tor all subsequent 
calculations. 

B . 2 Ph ase I I - Subs tructure An a 1 ys i s- 

SAP I { J is a commonly available -finite element code 
used in the ELSHOK code to perform the modal analysis of 
each substructure (equipment model). The equipment 
i n s t a 1 1 a t i on is -f i r s t cast i n t o an ap p r op r i a t e -finite 
element representation and then SAP I L 1 is used to solve tor 
the desired mode shapes and -frequencies. This information 
is placed in a substructure mode tile along with the 
geometry, stiffness and mass information. The mass, mode, 
and frequency information were also used in the DDAM 
an a 1 ys i s f or c omp arat i v e p u r p ose s . 

The final step in phase II is the execution of the 
PI CRUST code. This code prepares the input for phase III 
calculations by placing the abov e moda 1 inf orma. t i on on a 
substructure file. Additionally, it calculates influence 
coefficients cor r espond i ng to the base or support motions- of 
the substructure and accounts for the user-specified 
c on n e c t i v i t y . 
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Phase III - Time Integration 



Following the completion o+ Phases I and II, the 
she 1 1 -and-t 1 u i d +ile is merged with the substructure tile 
to form the input tile tor USLOB which is ELSHQK'* s time 
integration processor. USLOB allows spec i t i cat i on ot the 
charge weight and geometry. For the problems ot concern 
here, an exponentially decaying, empirical model was used 
to represent the explosive shock source. It i s ot the 
t orm : 



P 1 < R ) = Ki(W 1/3 /R>^ 2 ex p <-t/$o> <4) 



where ; 



Ki 

Ka 

t 

UJ 






K 3 



K 4 



the incident pressure at a radius R trom 
the explosion (on the hull) 

a multiplicative c on s t an t t or i n c i de n t 
pressur e 

a spatial decay constant tor incident 
pr essur e 

t ime at ter arrival ot shock wave at point 
ot interest 

weight ot spherical charge 
Ks id 1/8 < UJ 1/S / R >K 4 

multiplicative constant -for time constant 
o-f exponential decay 

spatial decay constant -for time constant 
of exponential decay 



In all the cases, the geometry was such as to maintain con- 
stant energy flux according to eq.<3). The integration in 
time is done using a modified version of the Runge-Kutta 
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n umer i c a 1 i n t e gr a. t i on sc h erne . Du r i ng the i n teqr a. t i on , 

LI S LOB ou t p u t s v e 1 oc i t y h i s t or i es f or u ser-sp e c i f i e d point s 
on the hull and substructure. 

B . 4 . Phase ID 1 - Plotting 

Ph ase I c ap ah i 1 i ties we r e not utilized in this in 
tigation. However, ELSHOK does have the capability to gen- 
erate plot files containing v e 1 oc i t y h i s t or y i n f or ma t i on f or 
p r oc e ss i n g on a graphics term i n a 1 . 

B . 5 . De f 1 e c t i on Ca 1 c u 1 a t i on s 

As mentioned pre v i ou s 1 y , ELSHOK p r odu c e s v e 1 oc i t y 
histories at the points of interest on the substructure. 

DDAM produces maximum deflections or forces at the same 
points for similar models. In order to make comparisons 
between the two methods, it was decided to convert the 
v e 1 oc i t y h i s t or i es f r om ELSHOK to de f 1 e c t i on h i s t or i e s u s i n g 
an . i n t egrator based on S i mpson x s 1/3 Rule. For each 
e q u i pme n t case, a separate p r ogr am was wr i 1 1 e n in BAS I C t o 
accomodate the variations in configuration between each 
model and the differences in output from the USLOB code due 
to point of interest specification. These programs are 
listed in Ap pend i x B . 



MODELS USED IN THE ANALYSIS 



IU. 



A . FINITE ELEMENT/FINITE DIFFERENCE MODELS 

The i dee. of re presen t i ng a 9 i yen dome. I n as. a. col 1 ec t i on 
of discrete elements or values was not -first appl i ed in 
■finite e 1 emen t/f ini t e difference methods. Early mathemati- 
cians estimated the value of Tt to nearly 48 decimal places 
by representing the circle with a polygon of a finitely 
large number of sides. The polygon used by the mathemati- 
cians is a model of the circle they wished to represent. 

In modern numerical methods, it is realized that very few 
problems, in the real world can be analyzed without the use 
of simplifying mode 1 s . B0S0R4 , the finite difference code 
in ELSHOK used to describe the dynamic behavior of the sub- 
marine hull, requires a greatly simplified model of its 
real wor Id counter p ar t . Similar! y , bo t h DDAM an d ELSHOK 
use a substructure model to capture the fixed-base 
properties of installed equipment. 

B. SUBMARINE HULL MODELING USING B0S0R4 

Ac t u a 1 su bmar i n e hulls ar e v ery c omp lex structures. 
However, they do possess certain characteristics that can 
be used to simplify their modeling. From the aspect of 
submarine hull shock response, the hull is a free -free ring 
stiffened cylinder with suitable end closures. 



The 



internal weights of machinery and other outfitting must be 
accounted for but these items can be smeared into their 
adjacent hull structure. The model thus becomes a cylinder 
whose cross sections reflect the density of the submarine 
at a given frame. Since a cylinder is a surface of 
revolution, the properties of each cross section need only 
be specified at one point. The reader will notice that any 
hard structure serving to stiffen the hull in a non- 
ax i symmetric sense, such as non-circular frames or massive 
machinery foundations, can not be represented in this 
fashion. If these structures are in the area of interest 
in the hull, they should be modeled as internal structures 
in phase II of EL3H0K. 

The B0SGR4 code allows the submarine hull model to be 
constructed by specifying the properties of the submarine 
hull along a line of revolution. A numerical complication 
enters the problem when dealing with models of submarines 
because the structure often repeats itself regularly along 
its axis of revolution. With such a configuration, there 
are many modes in which the motion of the stiffeners is of 
smal 1 ampl i tude compared to that of the shel 1 . The bays 
between the rings vibrate at frequencies which may 
approximate those corresponding to a s i mp 1 y-suppor ted 
cylinder of the same geometry as the bay. Multiple Or- 
el ose 1 y-sp ace d eigenvalues correspond to modes in which 
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one or more of the bays is vibrating while others are 
unaffected. To avoid eigenvalue difficulties, the pru de n t 
measure is to analyse as small a segment of shell as 
p oss i b 1 e in or de r to avoid numerical d i f f i c u 1 tie s w i t h 
multiple, or closely-spaced eigenvalues as indicated by 
Bushnell CRe-f. 9 3 . In submarine shock analysis, the 
interest is in the gross response contributions of the 
entire hull to a given section of interest. The authors of 
ELSHOK have allowed for this through the use of two models. 

A full model is used to obtain the gross contributions 
of the entire hull structure in the N = 0< torsional) and 
N = 1 < tr ansi at i onal and whipping) modes. This will include 

such things as bulkheads, heavy ring stiffeners, and end 
closures. A compartment model is then used to obtain more 
detailed information about the particular section of 
interest. In this way, responses away from the section of 
interest do not have to be carried through the remaining 
calculations since they do not affect the response of this 
area anyway. 

Due to security restrictions placed in information re- 
garding the construction techniques and arrangements of U.S. 
nuclear submarines, it was not possible to model an actual 
submarine tor use in the analysis work contained in this 
thesis. With the aid of Weidlinger Associates of New York, 
however, a model resembl ing a general submarine hull struc- 
ture of 6908 long tons displacement was obtained. This 
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model is available through the U.S. Naval Postgraduate 
School but is not included here so that restrictions do not 
have to be placed on the circulation of this work. 

Figure 2 is a schematic depiction of the full subma- 
rine. The submarine being modeled has been cut into 14? 
slices along the longitudinal axis. In the model, in for- 
mat ion regarding Young's modulus, Poisson's ratio, mass 
density, thermal expansion coefficient, plating half- 
thickness and whether or not stiffeners are smeared is pro- 
vided for each one of these points. Twenty-six discrete 
rings are represented along the axis of the model , twenty 
of which are in the local area of interest. These rings 
represent f ramies in the real submarine where T-st i ffeners 
give increased stiffness to the hull. In actuality, these 
stiffeners are spaced fairly evenly along the entire model, 
but smearing them outside the area of interest will not 
appreciably affect the vibrational characteristics in the 
area of interest. Using the B0S0R4 code, modal analysis 
using finite difference techniques is carried out. Table 
1 is a listing of the natural frequencies and wave numbers 
for each mode retained in the analysis. Seven N = 0 
torsional modes, and thirty-eight N = 1 rigid body 
d i sp 1 acemen t and whipping modes have been retained for the 
full mode 1 . 
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Figure 2 Full Model Used to Capture Gross Response of Submari 



MODE 


N 


FREQ 

(Hz) 


1 


1 


3.597583E-02 


3 


1 


2 . S772W5E+ 0 0 


5 


1 


1 . 29831 8E+ 01 


z -8 


1 


1 .301 155E+81 


Q 


i 

i 


2 . 305725E+0 1 


1 1 


1 


2 . 992289E+ 0 1 


1 3 


1 


3.41 1510E+01 


15 


1 


3 . 583254E+0 1 


1 7 


1 


3.71 80 95E+ 0 1 


1 9 


1 


3 . 847244E+0 1 


21 


1 


4 . 0 2 7 5 8 3 E + 0 1 


23 


1 


4 . 295375E + 0 1 


25 


1 


4 . 492591 E+0 1 


27 


1 


5.281 1 50 E+0 1 


zy 


1 


5 . 544798E+0 1 


31 


1 


5 . 994389E + 0 1 


33 


1 


8. 31 2 191 E+0 1 


35 


1 


8 . 5741 09E+0 1 


37 


1 


8 . 781 0 98E + 0 1 


3? 


0 


4.01 3475E-0 5 


41 


0 


1 . 485899E + 0 1 


43 


0 


3.021 883E+01 


45 


0 


4 . 57980 4E+0 1 


TABLE I 


- Wa 


ve Numbers ar 



iODE N FREQ 

< Hz ) 



•y 


i 


4 . 71 8345E-0 2 


•4 


1 


7.8901 12E+00 


& 


1 


1 . 793831 E+0 1 


o 


1 


1 . 897728E + 8 1 


1 0 


1 


2 . 72091 SE+0 1 


1 2 


1 


3 . 228725E + 0 1 


1 4 


1 


3 . 498258E+ 0 1 


1 8 


1 


3 . 590094E+0 1 


18 


1 


3 . 30 8358E + 0 1 


20 


1 


3 . V82V85E + 0 1 


22 


1 


4 . 0591 42E+0 1 


24 


1 


4 . 390 227E+0 1 


28 


1 


4.918127E+01 


23 


1 


5 . 289882E+ 0 1 


30 


1 


5. 634701 E+0 1 


32 


1 


6. 131 432E+0 1 


34 


1 


8 . 52 1 988E + 0 1 


38 


1 


8 . 8S2898E+ 0 1 


33 


1 


8 . 345927E+0 i 


40 


0 


7 . 584749E + 00 


42 


0 


2 . 342790E+0 1 


44 


0 


3 . 85595SE+ 8 1 



Frequencies of Full Model 



Figure 3 is a schematic representation of the compart- 
ment (local area of interest) model used in this analysis. 
In this model, all frames are modeled as discrete ring 
stiffeners. Since compatibility must be maintained between 
the two models, (compartment and full), there is a one to 
one correspondence in the area of interest between points 
in each model. The compar tmen t model is used to obtain the 
mode. 1 c on t r i bu t i on s c or r e sp on d i n q t o N > 2 < n orm a 1 , me r i d- 
ional, and circumferential displacements). As mentioned 
p r e v i ou s 1 y , de f orma t i on s or mo t i on s oc c u ring be yon d the 
region of interest should be unimportant for the N = 9 
breathing modes and modes corresponding to N l 2. Each 
bou n d i n g ring or bulkhead is v e r y s tiff in it s own plane 
and moves primarily in that plane for localized loading. 
Table 2 is a listing of the natural frequencies and 
c i rcumf erent i al wave numbers for each mode retained in the 
analysis. For the compar tmen t model, 25 modes corres- 
ponding to N =2, 24 modes corresponding to N = 3, 24 modes 
corresponding to N = 4 and 25 modes corresponding to N = 8 
have been retained. A total of 254 modes were found by 
BGSGR4 for the full and compartment models. Of these 
modes, 143 were retained for further calculations in phase 
II and III of ELSHOK. The 111 modes that were dropped 
affect responses primarily outside of the area of interest. 
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Figure 3 Compartment Model Depicting Area of Interest 
(Note: all major stiffeners included) 
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TABLE II - Dave Numbers and Frequencies o i Compartment Model 



Figure 4 is a schematic of the general shell structure 
model and notation conventions used in the remainder of the 
ELSHOK an a 1 y s i s . Two c oor d i n a t e svs t ems ar e u se d . X , Y an d 
2 refer to the submarine global coordinate frame and x,y 
and z refer to the substructure local coordinate frame. In 
both systems, the x-axes run from aft to forward and the 
z-axis orientation is related to the Z axis in the global 

system by the* angle c zq where 0 refers to the substructure . 
The variable s is used to locate points along the meridian 
of the full model , and u,v, and w are used for local shell 
displacements and motions, located circumferential ly by the 
angle 0. This concludes the discussion of the submarine 
model used in ELSHOK. 

One of the conveniences realized in using a D DAM 
an a 1 ys is to c ar r y ou t design checks- i s the absen c e of a 
su bmar i n e mode 1 in the c omp u t a t i on s . The de s i gn sh oc k 
values are intended to provide all inputs to the substruc- 
ture model. The penalty which must be weighed here is the 
lack of regard which this places on the peculiarities of a 
given submarine structure. On the other hand, the ELSHOK 
submarine model is difficult to construct requiring a great 
deal of i n f ormat i on an d sk ill top r op e r 1 y mode 1 a g i v e n 
hull, but the price is paid only once for each submarine 
class. All subsequent calculations in this thesis are 
carried out utilizing the same shell model calculations. 
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Figure 4 - She 1 1 -Subs true ture Con-figuration [Ref. 43 



C. SUBSTRUCTURE MODELING USING SAP IV 

In both DDAM and ELSHOK calculations, a substructure 
model is formulated to determine the modes, masses, and 
natural frequencies of the equipment system being design 
checked. ^ LJh e n DDAM was first de v eloped, a c omp u t e r p r ogr am 
was av a i Table wh i c h c ou 1 d so 1 ve the e i gen v a 1 u e p r ob 1 em f or 
up to twenty degrees of freedom. However, the scarcity of 
computer resources dictated that most users utilize hand 
computations for this purpose. Consequently, most early 
design chec k s we re limited in sc op e . DDAM i s we 11 suite d 
to tabular methods of computation but with the common 
a v a i 1 ab i 1 i t y of sma 1 1 c omp uter s today , h an d c omp u t a t i on s 
are no longer required, and the restrictions on degrees of 
freedom are largely removed. This is not to say that great 
numbers of degrees of freedom are required in every case. 

By careful examination of the equipment installation being 
checked, the major response contributions can be captured 
using a small number of modes. In the cases considered, in 
this thesis, emphasis was placed on using equivalent models 
f or DDAM an d ELSHOK r a t h e r t h an c 1 ose 1 y mode ling a real 
component in a submarine as it would be in each case in 
actual practice. The interest here 1 ies in how the results 
compare, for the same models. 

Phase II of ELSHOK utilizes the SAP IV finite element 
code to model submar i ne- i ns tal 1 ed equipment. Because 



ELSHOK i s wr i 1 1 e n utilizing p r ogr am modu 1 e s wh i c h execute 
independently and exchange data through output tiles, the 
SAP I 'v 1 code could be used to run the modal analysis of 
equ i pme n t mode 1 s f or DDAM an d ELSHOK t h u s e n su ring equ i v a 1 
ence in input to both methods. The models selected tor 
an a i y s i s we re kept sma 1 1 so that c omp ar i son s c ou 1 d be made 
without excessive compl i cat ion. 

Three equipment models have been examined tor this 
thesis. Case I is a cantilever foundation mode 1 with a 
1008 lb weight attached to its tree end. Cases II and III 
are a simple beam foundation supporting two weights at 
center span and having attachment points on two separate, 
discretering stiffeners. 

C . 1 . Substructure Case I Model 

The case I model development is depicted in 
figure 5. The model represents * I860 lb valve supported b 
a c an t i 1 e v e r f ou n da t i on . The f ou n da t i on i s bu i 1 t u p f r om 
A I SC C 18x30 steel channels [Ref. 113 and the valve is 
rigidly fixed to the free end. Since the intent was not t 
q u a 1 i f y any p ar t i c u 1 ar v a 1 v e design, the entire v a 1 v e h as 
been modeled as a particle (lumped mass). The beam mass i 
represented by six lumped masses. The foundation is 
designed to be fixed to a discrete-ring stiffener in the 
submarine in an upright position. Although 3-dimensional 
met o ion has been allowed, the maximum deflection occurs in 
the athwart ship direction in response to a side -on shock 
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1 Dad i n g . Th i s 



equ i pme n t i nstdl ! 3. t i on . 



C . 2 , bu bs tructure Case I I Mode 1 



The case II model development is depicted in 
figure 6 . The model represents a foundation constructed 
f r om A I SC W 27x 1 77 I -beams [Ref. 11] supporting two 1 0 , 0 8 8 
1 b we i q h ts . The f oundat i on struc tore i s represen ted by 1 8 
beam elements and 1? masses. It spans two discrete 
s t i f f e Pi e r - 5 t o wh i c h its ends ar e f i x e d . Th i s mode 1 i s 
r epre sen tat i v e of a large p ump or t u r bomach i ner y 
i nstal 1 at i on . 

C . 3 . Subs tructure Case III Mode 1 

The case III model is similar to the case II 
model except the foundation structure has been changed to 
A I SC UJ 27x 144 I -beams [Ref. 11] an d the su p p or ted we i qh t s 
have been reduced to 5,808 1b apiece. This model 
represents an intermediate weight equipment installation 
such as a main feed pump. 
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AISC C 18X30 Steel Foundation 
Properties 



I xx = 333.1 in. 4 

1 yy = 206.0 in . 4 

J = 540 .0 in . 4 



web thickness 
■f 1 ange th i ckness 
area o-f channel 



0 . 673 i n . 

0 . 436 i n . 

8.82 in. 2 





Section A-A 





REAL WORLD FOUNDATION 
AND MASS 



FINITE ELEMENT MODEL 
FOR DDAM & ELSHOK 



Figure 5 - Case I Equipment and Finite Element Model 
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AISC W 27X177 Steel Foundation 
Properties 

I xx = 6740 in. 4 
I yy = 556 in. 4 

J = 20.1 in. 4 

web thickness = 0.725 in. 

■flange thickness = 1.190 in. 

area o-f beam = 52.2 in. 2 



CASE II 



Ik 



H 






Section A-A 



AISC W 27X114 Steel Foundation 
Properties 

*xx — 4090 in. 4 
I yy = 159 in. 4 

J = 7.36 in. 4 

web thickness = 0.570 in. 
■flange thickness = 0.932 in. 
area o-f beam =33.6 in. 2 



CASE I I I 




Section A-A 




Figure 6 - Case II & III Finite Element Model 
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U . ANALYS I S 



The goal of the analysis procedure used in this inves- 
tigation is two-told. The primary objective is to analyze 
equivalent substructures using DDAM and ELSHOK so a com- 
parison of the maximum predicted deflections can be made. 
Secondly, any hu 1 1 /subs true ture interaction effects are to 
be noted. The analysis is complicated by several factors. 
All the calculations are based on constant energy flux which 
is dependent on charge weight and standoff. Since DDAM 
shock inputs ultimately represent the results of explosive 
shock tests, the geometry of the "analytical charge" is 
invisible to the user and assumed in the empirical shock 
spectrum. ELSHOK uses inputs of charge weight and standoff 
to calculate shock loading by eq. (4) so the effects of 
their variation will change the transient response of the 
substructure even though the shock intensity is constant. 

In this analysis, three sets of charge weight and standoff 
were used. A further compl i cat ion results due to the fact 
that DDAM only can be used to calculate maximum relative 
deflections or forces. ELSHOK calculates transient 
velocities incurred by the model which must be converted to 
maximum relative deflections. 
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CASE ANALYSIS PROCEDURE 



A . 

Eac h c ase an a 1 y s i s was be gu n by f i rst c on s t r u c t i n g the 
equipment models to be investigated using SAP I CL A dynamic 
analysis was performed to -find the natural -frequencies and 
mode shapes o-f the model . This information was saved along 
w i th the mass matr i x for subsequen t DOAM cal cu 1 at i on s . 

Using PI CRUST, the SAP I C* data was then reduced to a suit- 
able form to be merged with the hull structure/fluid data 
from Phase I calculations. Finally, all the resulting data 
were integrated using the USLOB code to compute the t r an - 
si ent v e 1 oc i t y p r of i 1 e s f or each equipment i n s t a 1 1 a t i on 
c on f i gu rat i on an d c h ar ge we i gh t/s tan dof f se t . 

The PI CRUST code was used to specify the installation 
configuration of each model. Model coordinate system orien- 
tation to the hull system is determined by the angle cc 0 j n 
fig. (4). The location of the equipment along the longi- 
tudinal axis is also specified. For example, in the case I 
model , the angle is 279 degrees and the base attachment is 
at frame ?5 in the submarine model (within the compartment 
model). For the case II and III models, two different 
orientations were investigated, one with the model mounted 
athwart ship and the other with the model mounted verti- 
cally, i.e., shock input from the side and bottom of the 
submar i ne . 
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The US LOB code allows input of charge weight and 
standoff. To investigate how these parameters affect the 
transient response of the substructure , three sets of val- 
ues were used tor these parameters based on three different 
charge weights. The parameters chosen are listed in the 
toll ow i ng table. 



TABLE III - Charge Weights and Standoffs for Analysis 

After selection of all factors affecting the geometry 
of the problem, the time step increment and integration 
limits were specified. In each case, enough time steps 
were chosen to capture the peak response amplitudes. This 
number was found by trial and error. 

After calculation of the velocity profiles for a given 
model, this information was integrated using Simpson's 1/3 
rule to obtain a deflection history and the maximum deflec- 
tion response of the model. As alluded to earlier, the 
models were constructed so that suitable differences taken 
between the velocity profiles of designated points would 
yield the relative deflection, at any instant in time, of 



Charge Weight 
(lbs TNT) 



S t an doff Distance 
< I nc he s/Fee t > 



5,800 

10,000 

15,000 



1 ,414/1 17.9 
2 , 000/1 66 . 7 
2 , 450/204 . 2 
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the point in quest i on . Th i s p r oc e du r e u n f or t u n a t e 1 y r- e - 
1 ies, to a certain extent, on symmetry in the model . For 
instance, in the case I model , the de-flection o-f the weight 
on the end of the cantilever in the a thwart ship direction 
can be obtained by integrating the difference between the 
velocity histories of the base and the tip of the beam. 
Figures <7> and <S) show schematically the differences 
taken to calculate deflections for each mode 1 configura- 
tion. 




SHOCK 

WAVE 



^rel — ^base , athwar tsh i p 



^ t i p , athwar t sh i p 



Figure 7 - Relative Velocities Used to Calculate 
Relative Deflection Between Mass and 
Base for Case I 
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Vertical Relative Velocity of Nodes 3 & 6 wrt -fixed base 
Nodes (shock wave from ATH4ARTSHI P) 




Vertical Relative Velocity of Nodes 3 & 6 wrt fixed base 
Nodes (shock wave from below keel) 



^3,rel — V3 )Ver t - <Vi }U ert + V2,vert>/2 
^6,rel = ^6,vert “ <V4 ( yert + V5,vert>/2 



Figure 8 - Case II & III Relative Velocity Calculations for 
Athwartship and Vertical Shock Waves 
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After t hi e E L S H 0 K resu 1 t -=■ w e r e o b t a i n e d for e a. c h model , 
DDAM was used to calculate the maximum deflections using 
the inf orma t i on sa v e d f r om the SAP I U runs. A file wa s c on - 
strutted, using an editing program, to input the masses, 
natural frequencies, and mode shapes for the model, to the 
DDAM p r ogr ami . Du ring e x ecut ion, the p r ogr ami utilize s 
p r omp t s t o q u e r y the u se r abou t the t y p e of i n s t a 1 1 a t i on 
and shock input direction. The maxi mum relative deflection 
is calculated by the NRL formula CRef. 123: 

J m = abs( J j ,max> + A j ) 2 -Ji ..max* <5> 



where ; 

abs< J j ?max > = the largest modal deflection of a 
point on the model 

= modal relative deflection at a point on 
the model 

4ri = max i mum relative de f 1 e c t i on of a p o i n t on 
the mode 1 

This formula is used to reflect the fact that all modes do 
n o t experience their max i mum de 1 e c t i on s s i mu 1 t an e ou s 1 y . 

Figures <?>, (10), and (11) sumimarize the model 

configurations and shock input directions for cases I, II, 



an d III. 


Tabl e 


s IV, V, .and VI 


are summaries of the 


i n ~ 


format i on 


used 


in the DDAM and 


ELSHOK cal cu 1 at i ons 


of the 


c a ses. A 


•full 


case I analysis 


is presented in Appe 


n d i x C 



e 



4 ? 



1414 " 



50130 lb TNT 



< b) 





c ) 




- Sh oc k c on f i gu rat i on s t or Case I 



“.M 



Figure ? 



( *. ) 







TNT 



5 TNT 



TNT 



b TNT 



I I 
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< b ) 




Figure 11 - Shock configurations for Case III 



TABLE I 1 -, 1 



C a S h I I n u t s t ci D DAM and EL S H G K 



Numbe r 


of 


Degrees of Freedom 


1 8 


Number 


of 


Con nec ted Degree s of Free dom 




Numbe r 


of 


Unconnected Degrees of Free dom 


1 5 


Nurriber 


of 


I n t e r f ac e Po i n t s 


1 


Number 


of 


Frequencies I n An a 1 y s i s 


1 0 



Sub 


s tructure 


Na t u r 


3. 1 F r e q u e ri c 


its (cps) 


i 


445.41 


5 


7608 . 71 


8 12169.6 




1654. 85 


6 


9558 . 0 5 


9 12234.0 




1 343 . 07 


~7 


10454.9 


1 0 1 80 14.2 


4 


4037.05 









Number of Time Steps In the Analysis 1*:0 1 

Size of Time Steps 0.5E-6 sec. 



TABLE U - Case II Inputs to ODAH and ELSHQK 



Number of Degrees of Freedom 34 

Number of Connected Degrees of Freedom 4 

Number of Unconnected Degrees of Freedom 30 

Number of Interface Points 4 

Numbe r of Frequencies In Analysis 5 

Sub struct u r e Na t u r a 1 Frequenc i e s < c p s > 

1 48.7743 4 551.143 



2 392.390 5 5S3.723 

3 392.598 

Number of Time Steps In the Analysis 1201 

Size of Time Steps 



1 . 0E-5 



TABLE 'v'l 



- Case III Inputs to DDAM and ELSHQK 



Number o-f Degrees o-f Freedom 3 

Number o-f Connected Degrees o-f Freedom 

Number o-f Unconnected Degrees o-f Freedom 38 

Number of Interface Points 4 

Number o-f Frequencies In Analysis 18 



Substructure Natural Frequencies Ceps) 



1 52.7625 

2 388.0991 

3 331.861 

4 537.983 



5 578 . 582 

6 353.971 

7 1213.45 



8 1213.51 

9 1 483 . 73 

10 1492.67 



Number o-f Time Steps In the Analysis 
Size o-f Time Steps 



1 30 1 

2.8E-5 sec. 
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RESULTS 



A 1 thou gh t h e number of cases an d t h e i r p e r mu t a t i on s 
e x p 1 o red in this in y e s t i g a t i on are n o t n ume r o u s , the y h a y e 
been c on t r i y e d to ac c e n t u a t e the s i m i 1 ar i t i e s an d differ- 
ences between the DDAM and ELSHOK me thods of analysis. As 
men t i oned prey i ousl y , ELSHOK proy i de s a great deal more 
information about the shock response of submar i ne- i nstal led 
e q u i pme n t t h an DDAM . Howe y e r , DDAM i s not intended t o be a 
theoretical tool but rather is simply a design check. For 
t h i s r e ason , wh ere ELSHOK w ill predict different max i mum 
deflections depending on c h ar ge we i gh t an d s t an dof f , DDAM 
will predict on 1 y on e y a 1 u e f or a 1 1 y ar- i a t i on s . 

A. CASE I (1,000 LB CASE) 

The predicted responses of DDAM and ELSHOK compare well 
f or t his c ase . The DDAM c a 1 c u 1 a t i on y i e 1 ds a max i mum de - 
flection of 0.0143“ for the end of the cantilever beam when 
subjected to an athwart ship shock input, whereas the ELSHOK 
cal c u 1 a t i on predicted a 0.0131" de f 1 e c t i on f or a 10, 0 Q 0 1 b 
charge and a 2,000 inch standoff as in fig (?b>. The 
15,000 lb and 5,6 00 lb charges both caused slightly lower- 
deflections. In this case, DDAM is more conservative 
(about 9)0 than the ELSHOK calculation. This illustrates 
that where the equipment weight is only about 1,000 1b, 



DO 



little interaction is to be expected between the hull 
resp on se wh i c h i s r e 1 a t i v e 1 y 1 ow frequency, an d the hi gh e r 
■frequency equipment response. The highest hull frequency 
is only about 188 Hz. Coupling effects between the 
responses of the equipment mass, and the hull mass are 
i ns i gn i f i gan t due to the large ratio of the masse s . If the 

mass of on 1 y segment f i v e in the full mode 1 i s c on si de r e d 
to be the entire hull mass, then the effective mass of the 
hull is about 4108 slugs. The mass of the model is less 
t h an 2 . 9 s 1 u gs . This y i elds a mass ratio of abou t 8.9 8 8 73 . 
DDAM is based upon the approximation that the hull mass can 
be c on s i de red infinite wh i c h implies a mass ratio of 
0.08080... p r o v i d i n g no accou n t f or coupling o t h e r t h an 
t h rou gh t h e safe t y f ac t or bu i 1 t into its de s i gn sp e c t r urn . 
The close agreement between ELSHOK and DDAM is probably due 
to the lack of response coupl ing in both methods- of 
a n a 1 y sis. 

B. CASE I I (28 ,880 LB CASE ) 

Two equipment configurations were investigated in the 
ELSHOK an a 1 y se s f or t h i s c ase . In the p r ob 1 em ge ome t r y i n 
wh i c h t h e su bmar i n e is side 1 oade d , figure < 9c ) , the 
maximum deflection is again in good agreement with that 
predicted by DDAM; 8.2433 inches in the former versus 
8.2744 in the latter. However, the input motion to the 
structure is not primarily in the direction of maximum 
deflection in this instance. I n t h e se c on d c on f i gu r a t i on 
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of t h i s c ase , c or r e sp on d i n g t o f i gu r e ( ?d ) , 3. mu c h gr eate r 
de-flection < 1 33'< greater) is predicted by ELSHOK. There 
ar e two p oss i b 1 e e x p 1 an a t i on s . Some h u 1 1 /s true t u r e 
interaction is to be expected because the -fundamental 
natural frequency of the substructure is 48.3 Hz which 
falls in the middle of the entire range of hull frequen- 
cies. If, as in the former case f a hull effect i v e m a = s o f 
4180 slugs is assumed, and the model mass is 67.4 slugs in 
this case, then the ratio of the masses is increased to 
0.8184 or about 1.6 percent. In this example, the influ- 
ence of coupling on the response can be expected to amplify 
the expected result for an infinite mass hull system 
cons i derabl y . Also a source of uncertainty is the test data 
u p on wh i c h the sh oc k in p u t sp ectr urri u se d i n DDAM i s based . 
Most explosive tests are set up with the charge geometry of 
figure (?c> to avoid problems associated with bubble migra- 
tion to the hull due to bouyancy effects. In the second 
configuration, figure <9d>, the charge was situated as if 
it had detonated beneath the keel of the submarine, a 
realistic p oss i b i 1 i t y wh e n c on s i de r i n g mode r n we ap on s 
t a c t i c s bu t per h a p s n o t simulated w ell by e x p 1 o s i v e test- 
ing. This particular c on f i gu r a t i on su b j e c t s the s t r u c t u r e 
to greater shock input motions in its more flexible direc- 
t i on p r odu c i n g the larger de f 1 e c t i on s , h owe v e r DDAM 
provides the same inputs for both problems. 
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C • CASE III (10, 060 LB CASE) 

T wo equ i pm e n t conti q u r a t i on s we re in v e s t i gated in the 
ELSHGK analysis for this case. The results are very- 
similar to the results for case II. In the configuration 
of Fig. 18(a), a maximum deflection of 0.2184 inch was 
predicted and in the c on f i gu r a t i on of Fig 1 8 < b ) , a max i mum 
deflect! on of 8 . 692 i nch was p r e d i c t e d by ELSHGK wh ere as 
DOAM predicted a max i mum of 8.2743 i n c h f or bo t h c on f i g- 
urat ions. Again, the same possible causes can be cited for- 
th© large <158/0 difference in the second instance. In 
this case the fundamental frequency of the foundation is 

52.7 Hz. The mass of the equipment mode 1 in this case is 

36 .8 si u gs > i e 1 d i n g a mass rat i o be twee n the hull and 
substructure of about 8.68877. 

A summary of the calculated results is presented in 
Tab 1 e M II. With the limited c a 1 c u 1 a t i on s a v a i 1 ab 1 e i n t h i s 
work, no attempt could be made to deduce systematically the 
exact c au se for d i sc r e p an c i e s be twe e n ELSHGK an d DDAM i n 
cases II and III. 
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DDAM 



MI I CONCLUSIONS 



Since its adop t i on by the n ay a 1 sh oc k c ommu n i t y , 
has proyided an easy-to-use and cony en i en t means to design 
check e q u i pme n t i n s t a 1 1 a t i on s p r op ose d for su bmar i n e s . It 
is in de p endent of an y we 11 defined su bmar i n e s t r u c t u r a 1 
input and can be applied early in the overall design pro- 
cess before such details are well developed. However, the 
same qual i ties that al 1 ow for the f 1 ex i b i 1 i ty i n DDAM may 
also c on t r i bu t e t o p oss i b 1 e in ac c u r ac i e s wh en t h i s me t h od 
is applied to submarines of radically different size or 
design than the ones reflected in the empirical Design 
Spec tr urn . 

The w or k in this thesis i s i n no way c on s trued t o be 
an a 1 1 -e n c omp ass i n g e v a 1 u a t i on of DDAM an d ELSHOK . Fr om 
the infinite n umber of p oss i b 1 e e q u i pme n t c on f i gu r a t i on s 
and sizes, three simple mode 1 s h a v e been selected wh i c h e x - 
amine only the hull mounted equipment problem. The intent 
of this investigation has been to examine the capability of 
the two methods to predict the shock response of these 
equipments and not to rank one against the other. ELSHOK 
or any var i at i on of this code i s much too comp 1 ex and 
c omp u t e r resource i n t e n s i v e t o u se a s a de s i gn tool an d wa s 
not intended for this purpose. 
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I n the case of equ i pmen t whose response i s re 1 at i ye 1 y 
high frequency compared to the hull response, DDAM and 
ELSHOK deflection predictions are in good agreement. The 
maximum hull response frequency is limited to about 196 Hz 
for the 6999 LT submarine used in this work. The good cor- 
relation is attributed to the lack of interaction and 
coupling effects between the hull and substructure. When 
the response frequencies of the substructure approached 
those of the hull and the mass of the substructure became a 
s i gn i f i c an t p or t i on of t h e ov e r a 1 1 sys tern, the c omp ar i son s 
between ELSHOK and DDAM deflection predictions were not 
good . 

The indication s of this work ar e that DDAM doe s n o t 
correctly reflect the h u 1 1 -substructure response coupling 
amp 1 i f i c a t i on wh i c h bee ome s ap p ar e n t wh e n the su bs t r u c t u r e 
response is tuned to the hull response. Additionally, the 
emp i r i cal database upon wh i ch DDAM i s re 1 i an t for shock 
input may fail to represent the heavy equipment problem 
we 1 1 . 

More basic research remains to be done before many of 
the questions or criticisms posed by this thesis can be 
considered conclusive. Due to the scarcity of publ i shed 
reports describing the input data used to generate the 
de s i gn sh oc k values u se d by DDAM , the ap p 1 i c a t i on of this 
method as a general design qualification should be regarded 
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with caution. T h e b^.s i s of t h e Design bpt*c trurri us&d in 
DDAM sh ou 1 d be av a i 1 ab 1 e -for an a 1 y s i s an d re v i s i on . 

Add i t i on a 1 wo r k r ema i n s to be done t o date r m i n e wh ere 
the d i sc rep an c i e s be t we e n EL3K0K an d DDAM first bee ome s i g- 
n if i cant. The cases analyzed here represent both ends of 
the substructure response spectrum but lend little infor- 
mation to what happens in the transition region between the 
1 ow- i nteract i on an d h i gh -interaction f requenc y r e sp on se 
regimes. 

DDAM should be retained as a design check method. 
However, some modifications are recommended here. The 
p r e se n t De s i gn Sh oc k i n p u t s used i n DDAM sh ou 1 d be u p da t e d 
to reflect the incresed size of modern day submarines. If 
the detail s of the original f or mu 1 a t i on we r e k n own or 
revealed, the current values could be updated by generating 
new analytical data based on computer simulations. This 
new spectrum could then be verified by explosive testing. 

It i s ap p ar e n t f r om the hull -su bs t r u c t ur e c ou pi in g phe- 
nomenon noted here that a modified shock spectrum should be 
f ormu 1 ated for each new c 1 ass of submar i ne wh i ch wou 1 d 
reflect the peculiarities of this class versus the previous 
ones . 

At the pre se n t time, wh e n n ew su bmar i n e s c os t billion s 
of dollars and are by no means numerous, every effort must 
be made to ensure their survival in an underwater explosive 
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shou Id include tak i ng a. new 1 ook 



shock enu i ronmen t . Th i ■= 
at DDAM and ensuring it upholds the standards set by- 
shock pol icy. There is a great deal of tlexibil i tv 
ent in DDAM and together with modern numerical techn 
it can be updated to reflect current technology. 



r i a u a. i 
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APPENDIX A ; DDAM USER'S MANUAL AND PROGRAM LISTING 



The DDAM Pr ogr am contained herein p r ov i des an e asy 
means to perform a smal 1 to medium-scale DDAM analysis of a 
given internal equiprrient design. The program is written in 
I BM PC BAS I C 1 an gu age and i s not restricted t o su brriar i n e 
analyses only but includes the standard DDAM case 
p oss i b i 1 i t i e s . The sol u t i on a 1 go r i t hm f o 1 1 ows t h e me t h od 
cited in Ref . 2 . 

A. USING THE DDAM PROGRAM 
A . 1 . Op t i on Select! ons 

The DDAM program is menu-driven and written to be 
se 1 f -exp 1 anatory , Prior to starting the program, it must 
b e a v a i 1 ab 1 e on a diskette a 1 on g with the BAS I CA . COM 
progr am . 

To start the program type: 

BAS I CA < ENTER > 

< F3> DDAM < ENTER > 

<F2> 

The first command line loads and executes the BASIC inter- 
preter and the second line loads the DDAM program. The 
third line s t ar t s execut i on of the DDAM Pr ogr am . 

Up on starting, the title of the p r ogr am < DDAM ) i s 
written on the screen and several seconds later a message 
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a p p e a r a- desc r i b i n g t h e p r ov i de d functions. Foil ow the 
directions. 

F ou r execut i on op t i on s w i 1 1 a p pear f or sou r ce of 
data and program exit to the operating system , They are: 



[ 1 ] Input Mode Sh ap e s an d Mass Ma t r i x f or the 
probl em -from Keyboard 

[ 2 ] Input Mode Sh ap e s an d Mass Ma t r i x f r om a 
d i sk file 

C 3 ] Input a Mass Ma t r i x an d 8 1 i f f n e ss Ma t r i x 

-for previously formulated equations of motion 

[5] Exit to DOS 



No option t 4 ] i s present ; h owe u e r , p r ou i s i on s h au e be e n 
m ade to install a user -def i ned subroutine at line 25800 in 
the program. 

Op t i on [ 1 ] is useful f or sma 11 p r ob 1 em s and i s n o t 
recommended for situations inuoluing more than four degrees 
of freedom. It allows all data input to be carried out 
using the keyboard. In larger problems, the entry of mode 
sh ap e s be c ome s t e d i ou s and e r r or p r on e . 

Option [2] is r ecomrriended for most problems. It 
allows the use of a diskfile for input. The file is 
user-specified and utilizes a free-format data structure. 
The information must be present in the file in a set order 



wh i c h i s de sc r i be d later. 



Selection of option C 3 3 causes a message to be 
d i sp 1 aye d stating that the JACOB I p r oqr am rriu s t f i r s t be 
run to solve the eigenvalue problem for mode shapes and 
natural frequencies, Th i s p r oqr am listing -foil ows the DDAM 
listing in t h i s ap p e n d i x . Inc or p or a t i n g the JACOB I p r oqr am 
i n DDAM wou 1 d h a v e reduced s i gn i f i c an 1 1 y the size of p r o- 
b 1 em wh i c h could be h an died. The JACOB I p r oqr ami is in- 
cluded as a convenience to the reader and was not integral 
to the analysis o-f the case studies done in this thesis. 

I t s ope r a t i on i s de sc r i bed later. 

A . 2 . Op t i on [ 1 3 Ex ecut i on 

Selecting option C 1 3 allows input of the masses, 
modes, and natural frequencies o-f the problem -from the 
keyboard. The number of modes- and masses- are first input. 

A mass input subroutine queries the user interactively for 
the diagonal entries of the mass- matrix. After input is 
completed, the diagonal elements are displayed enabling the 
user to make corrections as necessary. A zero correction 
terminates mass matrix input. Next, the program asks for- 
input of the natural frequencies of the system. It r e - 
q u i res on e n a t u r a 1 frequenc y f or e ac h miode shape . Th e 
input format is similar to the one used to input the mass 
elements. The natural frequencies are required in units of 
radians/second. Each mode shape is then input. A-fter 
initializing each mode , a correct i on op t i on is p r ov i de d . 
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The program next calculates the participation factors and 
moda 1 we i gh t s f or t h e p r ob 1 erri . 

Fo 1 1 ow i n g these p r e 1 i m i n ar y step s , DDAM r e q u i re s 
i n -forma t i on to -fix the installation details -for the equip- 
ment and type of analysis to be carried out. The standard 
cases of Ref. 2 are allowed. Once' the case is fixed, the 
Design Values are. calculated and displayed as DA< 1 ) . . . DA t n 
where n is the number of modes selected for- the analysis. 
The Design Values are the equivalent static accelerations 
calculated from empirical formulas for each mode. 

Several output options are provided. Modal force 
or deflections or both can be selected. After selection, 
the appropriate calculations are made and the results dis- 
p 1 aye d . Initially a p r i n t op t i on was se p ar a t e 1 y p r ov i de d . 
However, in the present version, selection of the display 
to screen option with printer on will cause output to go t 
the screen and the printing device. Disabling the printer- 
will cause output to go only to the screen. Following 
the output options, choices are provided to change the 
equipment installation type and repeat the analysis or to 
exit the program . 

A . 3 . Option C23 Execution 

Selecting option C 2 3 allows input of masses, nat- 
ural frequencies, and mode shapes from a disk file. After 
the data is read in, the program execution is similar to 
that in option Cl]. Option [2] prompts the user for the 



d i men s } on of t h e mass me. t r i x , the n umber of mode s i n t h e 
analysis, and the name o-f the input file. The input file 
is a free -format ASCII file containing the masses, natural 
frequencies, and mode shapes for the problem. This file 
may be constructed from structural analyzer (SAP I CO output 
or any o t h e r p r ob 1 em f or mu 1 a t i on (Tie t h od using a su i t ab 1 e 

edit or p r.ogr ami . < EDL IN, the edit or p r ov i de d w i t h the "I BM 

Disk Operating System" is not recommended for this purpose 
due to its line orient a tier.) The format f or the input 
file foil ows : 

mass< 1 ) .... iriass(m) (slugs) 

omega < 1 > .... omiega(mi) (rad/s) 

p h i ( 1 , 1 ) .... p h i ( 1 , n ) 

p h i < m , 1 ) .... p h i < m , n ) 

where m = number of masses and 

n = number of modes in the analysis 

In the above format, the data is read in sequentially from 

left to right. The data is echoed to the screen .and/or < 

printer in order to provide a check for proper input 

□per at i on . 
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B. PRGGRAM LISTING -DOAN 



10 CLS: KEY OFF : LOCATE 3,1,8 

28 PRINT"********************************************************************** 
**♦*’ 

21 PRINT"** 

**■ 



38 PRINT"** 


DDDD 


DDDD 




fttt 


MhW 


**" 

48 PRINT"** 


DD D 


DD D 


A 


A 


ffl 


m m 


**■ 

58 PRINT"** 


DD D 


DD D 


AAAAAA&i 


m 


m m 


**■ 

68 PRINT"** 

V vfl 


DD D 


DD D 


A 


A 


m 


m 


" A 

78 PRINT"** 


DDDD 


DDDD 


A 


A 


m 


tf. 


**' 

88 PRINT"** 
**’ 

98 PRINT"** 


DYNAMIC 


DESIGN 


ANALYSIS 


METHOD 



**■ 



188 PRINT'** 

**■ 

118 PRINT"** 

**' 

128 PRM" ********************************************************************** 
***** 

121 DIN MASS( 48) 

122 DIM PHI (40, 48) 

123 DIM U<48) 

124 DIM NX(48,48) ,HX2(48,48) ,SUM1(48) ,SUM2(48) ,P<48) 

125 DIM UA<48) 

127 DIM V8(48) ,A8(48) ,DA(48) 

128 DIM AA'140) ,VA(48) ,F<48 ,48) ,DEFL<48) ,DELTA(48 ,40) 

138 FOR I = 1 TO 5888 : NEXT I 

148 CLS: LOCATE 1,1,8 

158 PRINT" The purpose oF this program is to allow the user a convenient raeth 
od" :PRINT 

168 PRINT'D)' which to carry out a small to medium sized DDAM analysis oF a give 
n" SPRINT 

178 PRINT‘piece oF equipment. Familiarity with NOf&FAL MODE ANALYSIS and some u 
n-" : PRINT 

188 PRIfTT" derstandi ng oF the theory behind DDAM is prerequisite to conFident us 
e" SPRINT 

198 PRINT" oF this program. The program is prompt driven but the user should hav 
e "SPRINT 

288 PRINT" ava i 1 abl e to him / her either :" sPRINT sPRINT 

218 PRINT" A.) Mode Shapes and a Mass Matrix For the" 

228 PRINT" problem oF interest." sPRINT sPRINT 

238 PRINT" B.) A Mass Matrix and StiFFness Matrix From" 
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248 PRINT' previously formulated EQUATIONS of MOTION 8 : PRINT 

: PRINT: PRINT 



250 PRINT' ( Press ANf key to continue )':LGCATE 25,1,0 

260 At = INKErt: IF At = " THEM 260 
278 C-LS: 



280 LOCATE 5,1,0 

290 PRINT' In order to initialize the program input, choose one of the follow 
ing' [PRINT 

300 PRJNT’menu selections =—=== =>' :PRINT :PR!NT 



318 LOCATE 10,20,0 
320 LOCATE 11,20,0 

325 LOCATE 13,20,0 

326 LOCATE 14,20,3 
330 LOCATE 16,20,0 
340 LOCATE 17,20,0 
350 LOCATE 18,20,0 
369 LOCATE 20,20,0 



PRINT* 1 1 1 Input Mode Shapes and Hass Matrix" 

PRINT' for the problem from keyboard* 

PRINT* E 23 Input Mode Shapes and Mass Matrix* 

PRINT' from a disk file ’ 

PRINT't33 Input a Mass Matrix and Stiffness Matrix" 
PRINT' for previously formulated equations of' 
PRINT' motion* 

PRINT' t 51 Exit to DOS' 



370 A* = INKEYt : IF At >= CHRt(49) AND At <= CHRt<53) GOTO 488 
380 IF At = " GOTO 378 
390 8EEP : GOTO 370 



438 8 = VAL(At) : CLS 



410 ON B GOSUB 5888,15008,28080,25308,500 
500 SYSTEM :END 



688 '************************************************************************** 



601 ' SUBROUTINE TO INPUT MASS MATRIX 

683 ' ************************************************************************** 

605 INPUT 'Enter the dimension of the Mass Matrix => N 

610 IF N <= 28 GOTO 622 

620 PRINT :PR1NT:PR1NT'***** PROBLEM SIZE IS LIMITED TO 20 X 20 MATRICES **** 
* ' :PR1NT :PRINT: GOTO 688 

621 PRINT SPRINT: PRINT 

622 INPUT 'Enter the number of modes you wish to use => ";M 

638 CLS: PRINT'***** Enter the elements of the DIAGONAL Mass Matrix *»***' :PRI 
NT -.PRINT 



648 PRINT' If you make an error you will have an opportun-' 

650 PRINT' ity later to correct it.':PRINT:PRINT 

670 FOR 1 = 1 TO N 

688 PRINT ' MCjIj'.'jIj') = ';:INPUT' ', MASS(I) 

690 IF (1 = 5) OR (I = 10) OR (I = 15) THEN CLS 
700 IF MASSED <> 0 GOTO 730 

710 PRINT" ***** MASS MATRIX MUST BE POSITIVE DEFINITE — INPUT NON-ZERO VALU 
E ****' 

720 PRINT: GOTO 680 
730 PRINT: NEXT I 
740 CLS 

750 FOR I = 1 TO N 

760 PRINT 'M( ' ; I ; ' , ' j I j ' ) = ';MSSU), 

770 NEXT I 

780 PRINT: PRINT: PRINT 

790 PRINT'Enter diagonal element number to correct mistaken entries or " 

830 PRINT: INPUT ' Enter < ZERO > to continue =>', RESPONSE 
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819 IF RESPONSE = 8 THEN RETURN 

320 IF RESPONSE <= 8 OR RESPONSE > N THEN CLS : GOTO 758 ELSE 1 = RESPONSE 
850 PRINT ; PRINT :1NPUT ‘Enter correct element value =) M< I ) 

380 PRINT :PRINT'MC ;I , • j I ; 1 ) CHANGED TO ';MASS(!) 

878 FOR I = 1 TO 1500 : NEXT 1 : CLS : GOTO 758 
880 RETURN 

1880 '************************************************************************* 
1010 " SUBROUTINE TO INPUT MODE SHAPES 

1020 '************************************************************************* 
1050 CLS: LOCATE 5,1,0 

1060 PRINT'***** Enter OMEGA's corresponding to Mode Shapes *****' : PRINT 
1070 PRINT' < smallest to largest )*:PRINT:PRINT 

1080 FOR 1 = 1 TO M 

1090 PRINT ' omegaCjl;') = " ; :1NPUT' ',«<]) 

1100 IF < 1=5) OR (1=10) OR (1=15) THEN CLS 
1110 NEXT 1 

1120 CLS : FOR 1 = 1 TO M 

1130 PRIMT'QMEGAC ;1 ;") = *;W(I), 

1149 NEXT 1 

1150 PR1NT:PR1NT:PR1NT 

1160 PRINT 'Enter Frequency number to correct mistaken entries or ':PR1NT 
1170 INPUT ' Enter < ZERO ) to continue =>■ .RESPONSE 

1180 IF RESPONSE = 0 THEN GOTO 1230 

1190 IF RESPONSE <= 0 OR RESPONSE ) N THEN CLS : GOTO 1120 ELSE 1 = RESPONSE 
1200 PR!NT;PR1NT:1NPUT 'Enter correct Frequency value => ', U(l) 

1210 PRINT :PR1NT'0MEGA(* ;I ;*) changed to ';U(D 
1228 FOR 1 = 1 TO 2000 : NEXT 1 : GOTO 1128 
1230 CLS: LOCATE 5,1,0 

1240 PRINT'***** Enter EIGENVECTOR'S corresponding to Mode Shapes ***♦*■ :PR1NT 
1250 PRINT’ ( begin with Mode One )" :PRINT:PR!NT 

1260 FOR 1 = 1 TO M 
1270 FOR J = 1 TO N 

1280 PRINT' MODE SHAPE * ; I j 1 PhiCjJ;') = *;: INPUT’ * , PHI ( I , J) 

1290 IF <J=5) OR (J=10) OR <J=15) THEN CLS : LOCATE 5,1,0 

1300 NEXT J 

1305 CLS : PRINT '.PRINT 

1318 PRINT '*************** MODE NUM8ER *;1;’ ***************' 

1320 PRINT :PRINT : FOR J = 1 TO N 

1330 PRINT* Phi('iJj') = *;PH1(1,J), 

1340 NEXT J 

1350 PRINT: PRINT 

1360 PRINT'Enter element number in mode to change or’iPRINT 

1370 INPUT' Enter ( ZERO > to continue => ', RESPONSE 

1380 IF RESPONSE = 0 THEN GOTO 1430 

1390 IF RESPONSE <= 0 OR RESPONSE ) N THEN CLSiGOTO 1310 ELSE J = RESPONSE 

1400 PRINT: INPUT' En t er correct element value =) , ;PHI(I,J): PRINT 

1410 PRINT' Ph i (' ;J;*) CHANGED TO ' jPHKl ,J) 

1420 FOR K = 1 TO 2000 : NEXT K : CLS : GOTO 1310 

1430 CLS:NEXT 1 
1440 RETURN 

5000 •'I************************************************************************* 



73 



5805 ' SUBROUTINE TO DRIVE NODE/HASS INPUT 
5010 '************************************************************************** 
5020 LOCATE 5,1,0 

5030 PRINT' The program has been initialized to allow user input of the'iPRINT 
5040 PRINT'Mode Shapes and Mass Matrix. If you desire one o-f the other op-':PR 
ItTT 

5050 PRINT* t i ons , PRESS < escape key >; otherwise =>' : PRIt-JT : PRINT 

5055 PRINT' PRESS < space bar ) to continue )’ 

5960 A$ = INKEY* : IF A$ = CHRK27) GOTO 278 

5070 IF M O CHRK32) GOTO 5060 

5075 CIS : LOCATE 5,1,0 : GOSUB 600 

5080 CLS : LOCATE 5,1,0 : GOSUB 1000 

5090 CLS : GOSUB 6808 

5100 CLS: LOCATE 25,1,0 

5110 PRINT 

5120 GOSUB 8000 

5130 CLS: GOSUB 8200 

5140 CLS : GOSUB 10008 

6000 ' ************************************************************************* 
6010 ' SUBROUTINE TO FIND PARTICIPATION 
6828 ' FACTORS 

6030 '************************************************************************* 
6048 LOCATE 25,1,0 
6050 T$(l) = TIME! 

6060 PRINT' CALCULATING PARTICIPATION FACTORS = start : ';TIHE* 

6080 FOR I = 1 TO M 

6090 FOR J = 1 TO N 

7080 MX(I,J> = MASS<J> * PHKI.J) 

7810 MX2(I,J> = MX(1,J> * PHI < I , J) 

7020 NEXT J 

7830 NEXT I 
7040 FOR I = 1 TO M 
7053 SUMl(l) = MX< 1,1) 

7060 SUM2-I) = MX2(I,1> 

7078 FOR J = 2 TO N 

7080 SUNK I ) = SUMKI) ♦ HX( I , J) 

7090 SUM2(I) = SUM2(I> ♦ NX2<I,J) 

7108 NEXT J 

7110 NEXT I 

7120 FOR I = 1 TO M 

7130 P(I) = SUNKI) / SU12(I> 

7140 NEXT I 

7150 FOR I = 1 TO H 

7160 MBAR(l) = P(I) * SUNKI) 

7170 NEXT I 

7180 LOCATE 25,60,0 : PRINT' done: '{TIMES 
7190 DUMMY = M 

7200 IF DUMMY > 10 THEN KX = M/2 ELSE KX = M 
7210 L=1:J=0 
7215 LOCATE 5,1,0 

7228 PRINT'***** PARTICIPATION FACTORS PNO MODAL MASSES *****' ;PR1NT: PRINT 
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7221 LPRINT 1 ***** P FACTORS AND MODAL MASSES »****■ 

7222 LPRIfTT 



7230 FOR I = L TO KX 

7240 PR I ITT USING PP$ j I ,P<I) ,MBAR<I) 

7241 LPRINT USING PM j I ,P< I ) ,MBAR< I ) 

7258 J = J 4 1 

7260 NEXT 1 

7270 IF DUmf > 10 AND J < 11 THEN A* = INKEY* : GOTO 7280 ELSE GOTO 7311 

7280 PRINT :PR]NT:PRINT 

7281 LPRINT: LPRINT 

7298 PRINT 1 PRESS ANY KEY TO CONTINUE... 1 

7399 IF A* = 11 THEN GOTO 7380 

7318 KX = M : L = 18: CLS : GOTO 7215 

7311 PR]NT:PR2NT:PR1NT‘ PRESS ANY’ KEY TO CONTINUE. 

7328 A* = INKEY* : IF A* = 11 THEN GOTO 7320 
7330 RETURN 

3000 '************************************************************************* 
8010 ' SUBROUTINE TO CALCULATE MODAL 

8020 ' WEIGHTS 

8030 '************************************************************************* 
8835 CLS: LOCATE 5,1,0 

8840 PRINT 1 ***** CALCULATED MODAL WEIGHTS **»**■ 

8041 LPRINT 1 ***** CALCULATED MODAL WEIGHTS »****” 

8045 WJ* = "tttt tttttt.ttttStttttt 1 

8046 PRINT : PRINT 1 N WEIGHT IN KIPS’ 

8047 LPRINT 1 N WEIGHT IN KIPS 1 
8050 PRINT 

8851 LPRINT 
8078 FOR 1 = 1 TO M 

8088 WAN) =336 * SUMKJ) * P(I) / I860 
8898 PRINT USING UU$ ; I, WAN) 

8091 LPRINT USING WJ* j I,WA< I) 

8895 NEXT I 

8108 LPRINT :LPRINT :PRINT :PRINT :PRINT* Press any ke 

y to continue. . .' 

8118 A* = INKEY*: IF A* = 11 THEN GOTO 8110 
8128 RETURN 

8200 •’************************************************************************* 
8210 ' SU8R0UTINE TO DETERMINE GEOMETRY 
8220 ' OF ANALYSIS AND A <0) , V (8) , D (A) 

8230 '************************************************************************* 
8240 PRINT” Sene in-format ion is now needed to complete the details about th 
e 1 

8250 PRINT'type of vessel the structure I am analyzing is installed in and whet 
her 1 

8268 PRINT'or not the analysis is to consider elastic or elastic-plastic defor- 
ma- 1 

8270 PRINT 1 t ions. 1 



7225 PRINT 1 N MODAL PARTICIPATION 

7226 PP* =”«« #OttMttr AAA 

7227 LPRINT" N MODAL PARTICIPATION 



MODAL MASS 1 

««.«tt#««r AAA1 

MODAL MASS 1 
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8280 PRINT : PRINT 

8290 PRINT'PLEASE ENTER ONE OF THE FOLLOWING =>' 



8300 LOCATE 9,20,0: 
8310 LOCATE 11,28,8 
8320 LOCATE 13,20,0 
8330 LOCATE 15,20,8 
8340 LOCATE 17,20,0 
3358 LOCATE 19,20,0 



( hull mounted system )* 

< deck mounted system )’ 

( shell plating mounted system )' 

< hull mounted system )' 

( deck mounted system )' 

( shell plating mounted system )' 



PRINT 1 !!] SUBMARINE 

PRINT' [21 SUBMARINE 

PRINT " C33 SUBMARINE 

PRINT' [4] SURFACE SHIP 

PRINT"! 51 SURFACE SHIP 

PRINT' [61 SURFACE SHIP 

83i0 A* = INKEY* : IF A*>=CHR*<49) AND A* <= CHR*<54> GOTO 8390 
8370 IF A* = " THEN GOTO 83(50 
8380 8EEP: GOTO 8360 
8398 8 = VAL<A*):CLS 

8408 PRINT' ENTER ==> III for ELASTIC ANALYSIS' :PRINT:PRINT 
8410 PRINT' !2I for ELASTIC - PLASTIC ANALYSIS ' 

8420 A* = INKEY* : IF A* )= CHR*<49) AND A* <= CHRt(50) GOTO 8450 
8438 IF A* = " THEN GOTO 8420 
8448 BEEP: GOTO 8428 
8458 C = VALIA*) 

8468 IF C=2 AND <8 = 3 OR 8 = 6) THEN GOTO 8478 ELSE GOTO 8598 

8478 PRINT : PRINT : PRINT : PRINT" ***** ELASTIC - PLASTIC OPTION NOT AVAILABLE FOR ** 

***' : PRINT 

8480 PRINT'***** SHELL PLATING MOUNTED SYSTEMS ***»*' 



3431 LOCATE 18,20,8 

3482 PRINT' PRESS ANY KEY TO CONTINUE .... 
8433 A* = INKEY* : IF A* =" THEN GOTO 8483 
8490 CLS ; GOTO 8240 
3500 CLS: LOCATE 5,1,8 



8518 PRINT' ENTER - 


— > 


!11 


for motion 


in VERTICAL 


direction' :PRINT 


8520 PRINT' 




!23 


for motion 


in ATHWARTSHIP 


direction' :PRINT 


8530 PRINT' 




[31 


for motion 


in FORE AND AFT direction' 



8540 A* = INKEY*: IF A* >= CHR*\49) AND A* <= CHR*(51) GOTO 8578 
8558 IF A* = " THEN GOTO 8548 
8560 BEEP: GOTO 3548 

3570 D = VALIA*) 

3571 IF B = 1 THEN LPRINT' SUBMARINE < hull mounted system )'; 

8572 IF 8 = 2 THEN LPRINT* SUBMARINE < deck mounted system )'; 

8573 IF 8 = 3 THEN LPRINT' SUBMARINE I shell plating mounted system )'; 

8574 IF B = 4 THEN LPRINT' SURFACE SHIP ( hull mounted system )'; 

8575 IF B = 5 THEN LPRINT' SURFACE SHIP < deck mounted system )’; 

8576 IF 8 = 6 THEN LPRINT' SURFACE SHIP ( shell plating mounted system )'; 

8577 IF C = 1 THEN LPRINT' for ELASTIC ANALYSIS' 

8578 IF C = 2 THEN LPRINT' for ELASTIC - PLASTIC ANALYSIS' 

8579 IF D = 1 THEN LPRINT ‘.LPRINT ' for motion in VERTICAL direction' sLPRINT 

8580 IF D = 2 THEN LPRINT ‘.LPRINT ' for motion in ATHWART SHIP direct ion' ‘.LPRINT 

8581 IF D = 3 THEN LPRINT :LPRIMT ' for motion in FORE AND AFT direction' ‘.LPRINT 
8590 FOR I * 1 TO N 

8600 IF 8 <> 1 THEN GOTO 8680 

8610 A0< 1) = 10.4* «488*UA<I»/<2BHIA<I))) 

8628 V8< I ) = 20! * <<430+UA(]))/<10B+UA(I))) 

8638 IF <D=1 OR D=2) AND C=1 THEN MULT1= 1 :MULT2=1 : GOTO 9118 ELSE 

8648 IF <D=3 AND C=l) THEN MULT1=.4:MULT2=.4 : GOTO 9110 ELSE 
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3659 

8660 

8670 

8680 

8690 

8700 

8710 

8720 

8730 

8740 

8758 

8760 

8778 

8780 

8790 

8888 

8818 

8820 

8830 

8848 

8850 

8868 

8878 

8880 

8898 

8900 

8910 

8920 

8930 

8940 

8950 

8960 

8978 

8980 

8998 

9000 

9818 

9020 

9830 

9048 

9850 

9068 

9880 

9090 

9100 

9118 

9120 

9138 

9140 

9158 

9160 



IF (Ol OR 0=2) AND 0=2 THEN MULT 1=1 :MULT2=.5 : GOTO 9118 ELSE 
MULTI =.4 : MULT2=,2 
GOTO 9110 

IF B <) 2 THEN GOTO 8780 

ABU) = 5.2 * ( <480 ♦ UAU))/(28 ♦ UA(I))) 

V8U) = 18! * < < 480 +UA< I ) )/< 1 8 0 +WA< 1 ) ) > 

IF (Ol AND C=l> THEN MULT1=1 :MULT2=1 : GOTO 9110 ELSE 

IF <02 AND C=l) THEN MULT1=2 :MULT2=2 : GOTO 9118 ELSE 

IF (D=3 AND C=l) THEN MULT1=.8:MULT2=.8 : GOTO 9110 ELSE 
IF (0=1 AND 02) THEN MULT1=1 :MULT2=.5 : GOTO 9110 ELSE 

IF (D=2 AND 02) THEN MULT1=J2 :MULT2=1 ! : GOTO 9118 ELSE 

MULTI = .8 : MULT2 = .4 
GOTO 9118 

IF 8 <> 3 THEN GOTO 8850 

ABU) = 52 * <<488+UA(l))/(20+UA(l))) 

V0U) = 100 * U488+UAU))/U00+UAU))) 

IF 0=1 THEN MULT1=1 :MULT2=1 : GOTO 9110 ELSE 
IF 0=2 THEN MULT 1= . 2 :MULT 2= . 2 : GOTO 9110 ELSE 
MULTI = .08 : MULT2 = .08 
GOTO 9110 

IF 8 <> 4 THEN GOTO 8958 

ABU) = 20 * (<37.5 ♦ UAU ) ) *< 1 2 ♦ UAU))/(6 + UAU) A 2)) 

U8U) = 68 * <U2 + UAU))/'<6 + UAU))) 

IF (D=l AND C=l) THEN MULT1=1 :MULT2=1 : GOTO 9110 ELSE 
IF (0=2 AND Ol) THEN MULT1=.4:MULT2=.4 : GOTO 9110 ELSE 

IF (0=3 AND 0=1) THEN MULT1=.2:MULT2=.2 : GOTO 9110 ELSE 

IF (0=1 AND 02) THEN MULTOl :MULT2=.5 : GOTO 9110 ELSE 
IF (0=2 AND 0=2) THEN MULT1=.4:MULT2=,2 : GOTO 9118 ELSE 
MULTI = .2 : MULT2 = .1 
GOTO 9110 

IF 8 <> 5 THEN GOTO 9858 

ABU) = 10 * ((37.5 4 UAU))*<12 ♦ UA(I))/(6 + UA(1) A 2)) 

00(1) = 30 * ((12 4 UA(l))/<6 4 UAU))) 

IF (Ol AND 0=1) THEN MULT1=1 :MULT2=1 : GOTO 9110 ELSE 
IF (0=2 AND 0=1) THEN MULT1=.4:MULT2=,4 : GOTO 9110 ELSE 

IF (03 AND 0=1) THEN MULT1=.4:MULT2=.4 : GOTO 9110 ELSE 

IF (0=1 AND 0=2) THEN MULT1=1 :MULT2=.5 : GOTO 9118 ELSE 
IF (0=2 AND 0=2) THEN MULT1=.4:MULT2=.2 : GOTO 9110 ELSE 
MULTI = .4 : MULT2 = .2 
GOTO 9110 

ABU) = 48 * ((37.5 4 UA(I))*<12+ UAU))/<6 4 UA(1)) A 2) 

V8U) = 128 * ((12 4 UA(I))/(6 + UAU))) 

IF Ol THEN MULT 1=1 :MULT2=1 : GOTO 9118 ELSE 
IF 02 THEN MULT 1 = . 2 :MULT2= . 2 : GOTO 9110 ELSE 
MULTI = .1 : MULT2 = .1 

V0(I) = MULT2 * 018(1) : ABU) = MULTI * ABU) : NEXT 1 
FOR I = 1 TO M 

AAU) = ABU) * 386 
VAU) = V8U) * U(I) 

IF A8S(AA(D) > A8S(VAU)) THEN DAU) = WU) ELSE DAU) = AAU) 

IF ABS(DAU)) ( 2316 THEN DAU) = 2316 
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9170 NEXT 1 

917! CLS: LOCATE 5,1,0 

9172 FOR 1 = 1 TO Ms PRINT "DA("jl ;") = ' ;DAd) , 

9173 NEXT 1 

9174 PRINT sPRINTsPRINT" Press any key to continue . .." 

9175 A* = INKEYIsIF Al = "" THEN GOTO 9175 ELSE CLS : LOCATE 5,1,0 

9176 FOR I = 1 TO M : LPRINT : L PRINT ’DAC;I;‘) = ‘ ;OA<D: NEXT I : LPRINT 
9180 RETURN 

10090 '************************************************************************ 
10010 ' SUBROUTINE TO CALCULATE MODAL 
10920 7 F 0 R C E.S 0 R 0 E F L E C T I 0 N S 

10030 '************************************************************************ 
13040 LOCATE 5,1,0 

10050 PRINT" Almost done.... Please ENTER the number of the type of informat 
ion" 

10060 PRlNT:PRINT"you desire as output ==)’ 



10070 LOCATE 10,1,0 


:PR1NT" 


m 


MODAL FORCES" 


10088 LOCATE 12,1,0 


:PR3NT" 


121 


MODAL DEFLECT ISIS" 


10090 LOCATE 14,1,0 


SPRINT" 


[31 


FORCES & DEFLECT I ONS 1 


10095 LOCATE 16,1,0 


SPRINT" 


[41 


CHANGE ANALYSIS TYPE 


10096 LOCATE 18,1,0 


SPRINT" 


E5I EXIT TO DOS" 



10100 Al = 1NKEYI: IF Al >= CHRK49) AND Al <= CHRK53) THEN GOTO 10140 

10110 IF A* = ** THEN GOTO 10100 

10138 BEEP : GOTO 10100 

10140 FLAG1 = VAL(AI) : CLS 

10150 IF FLAG1 = 2 THEN GOTO 10220 

10155 IF FLAG1 = 4 THEN GOTO 16160 

10156 IF FLAG1 = 5 THEN SYSTEM 

10160 REM ****** CALCULATE MODAL FORCES ******** 

10170 FOR 1 = 1 TO M 
10180 FOR J = 1 TO N 

10198 F< 1 ,J) = P(I) * MASS(J) * PHI < I ,J) * DAd) 

10200 NEXT J 
10210 NEXT I 

10220 IF FLAG1 = 1 THEN GOTO 12080 

10230 REM ****** CALCULATE MODAL DEFLECTIONS ********* 

10240 FOR 1 = 1 TO M 
10250 FOR J = 1 TO N 

10260 DELTAd , J) = PHld.J) * Pd) * DAd) / Ud) A 2 

10290 NEXT J 

10300 NEXT 1 

10310 FOR J = 1 TO N 

10320 BIGGEST = 0 

10330 FOR I = 1 TO M 

10340 IF BIGGEST < A8S< DELTAd ,J>> THEN BIGGEST = ABStDELTAd ,J)) : K=1 

10350 NEXT 1 

10360 S = 0 

10370 FOR 1 = 1 TO M 

10380 IF I = K THEN GOTO 10400 

10390 S = S + DELTAd, J) A 2 ■ 

10480 NEXT I 
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10410 OEFL(J) = ABS(D£LTA(K,J)) ♦ SQR<S) 
13420 NEXT J 



12000 REM ********* OUTPUT RESULTS *********** 

12018 CLS : LOCATE 5,1,0 

12020 PRINT ' ENTER 113 TO ECHO RESULTS TO SCREEN" :P 

RINTs PRINT 

12030 PRINT " 123 TO ECHO RESULTS TO PRINTER"! 

PRINT:PRINT 

12848 A* = MET* 

12050 IF A* >= CHR*<49) AND A* <= CHR*(50> THEN GOTO 12878 
12068 IF A* = "" THEN GOTO 12048 ELSE BEEP: GOTO 12340 
12070 FLAG2 = UALtA*) 

12080 IF FLAG2 = 2 THEN GOTO 12300 
12090 IF FLAG1 = 2 THEN GOTO 12190 
12100 FOR I = 1 TO M 

12118 CLS : LOCATE 5,1,8 : PRINT' *** N 0 D A L FORCES FOR MODE 
< " ; I ; ' ) libs] ***' : PRINT 

12111 LPRINTsLPRINT" *** M 0 0 A L FORCES FOR MODE < ";I;" ) tlb 
s3 ***' : LPRINT 



12112 

12114 

12115 
12120 

12130 

12131 
12140 
12159 
12168 
12178 
12180 
12190 
12208 
12218 
I N T 
12211 
N T 
12215 
12220 

12230 

12231 

12232 

12233 

12235 

12236 
12248 
12250 
1225! 
12268 
12270 
12288 
12290 



LPRINT " N LPRINT 
PRINT " N " : PRINT 

ff$ = •#* iu.mittr''"'' 1 

FOR J = 1 TO N 

PRINT USING FF* ; J,F(I,J) 

LPRINT USING FF* ; J,F(I,J) 

NEXT J 

PRINT -.PRINT: PRINT :PRINT 

PRINT* Press any key to continue ...* 

A* = INKEY*: IF A* = " THEN GOTO 12170 

NEXT I 

IF FLAG1 = 1 THEN GOTO 12280 
CLS: LOCATE 5,1,0 

PRINT' *** DEFLECTIONS AT MASS ATTACHMENT P 
S ***' : PRINT :PRINT 

LPRINT : LPRINT'** DEFLECTIONS AT MASS ATTACH ME 
POINTS **' :LPR1NT 
DD* = a Hfl . ttMMMMM ' 

FOR I = 1 TO N 
PRINT'DELTA < ' ; I ; ' > = '; 

LPRINT' DELTA < ' ; I ; ' ) = '; 

LPRINT USING DD* ; DEFLU) ; 

LPRINT' in.' 

PRINT USING DD* ; DEFL(I) ; 

PRINT 'in.' 

NEXT I 

PRINT :PRINT:PR1NT 
LPRINTsLPRINT 

PRINT' Press any key to continue ...' 

A* = INKEY* : IF A* = " THEN GOTO 12270 
CLS : LOCATE 5,1,0 
GOTO 10080 



0 
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12300 PRINT ' OPTION NOT AVAILABLE YET;=> USE < Pr-tsc > KEY 1 
12310 FOR 1 = 1 TO 3000 : NEXT 1 GOTO 12000 

15000 ' ******************* ***************************************************** 
15010 •' SUBROUTINE TO INPUT MASS/MODES 
15020 ' F R 0 N D 1 S K F 1 L E 

15030 ' ************************************************************************ 
15040 CLS : LOCATE 5,1,0 

15050 PRINT 1 The program has been initialized to allow user input of the*:PRlNT 
1 5068 PRINT s Mode Shapes, Frequencies, and Masses uia an input -file." :FRINT 
15070 PRINT : PRINT" If you desire one of the other options," :PR]NT:PRINT 

15080 PRINT* PRESS fescape key); otherwise =>" -.PRINT : PRINT 

15090 PRINT* PRESS <space bar) to continue * 

15100 A* = INKEY* : IF A$ = CHRf-:27) GOTO 270 
15110 IF AJ 0 CHR*<32) GOTO 15100 
15120 CLS: LOCATE 5,1,0 

15130 INPUT ‘Enter the dimension of the Hass Matrix — ) *;N 
15140 IF N <= 50 THEN GOTO 15170 

15150 PRIIfT : PRINT : PRINT'***** PROBLEM SIZE IS LIMITED TO 58 DEG OF FREEDOM **** 

*• 

15168 FOR 1 = 1 TO 1888 : NEXT 1 : GOTO 15138 

15178 PRINT: PRINT: INPUT" Enter the number of modes you wish to use =) *;M 
15175 PRINT: PRINT: INPUT* Enter the input file name =) *,FIN$ 

15188 CLS 

15198 '************************************************************************ 
15288 ' INPUT MASS, FREQ, AND PHI 

15218 '************************************************************************ 
15228 OPEN FIN* FOR INPUT AS HI 
15240 FOR 1 = 1 TO N 
15258 INPUT ill, MASS(I) 

15268 NEXT 1 

15278 FOR 1 = 1 TO M 

15288 INPUT II, U( I) 

15298 U<I> «U(I> * 6.283185 

15388 NEXT I 

15318 FOR J = 1 TO N 

15338 FOR I = 1 TO M 

15348 'INPUT 81 , PHI (I, J) 

15350 NEXT I 
15368 NEXT J 

15378 7 ************************************************************************ 
15388 ' OUTPUT MASS, OMEGA, AND PHI 

15398 ' ************************************************************************ 
15488 LOCATE 5,1,8 : PRINT’MASS MATRIX DIAGONAL ELEMENTS * : PRINT 

15481 LPRINT* ***** MASS DIAGONAL ELEMENTS ****** 

15482 LPRINT 

15485 m = *88 88.88888 AAAA 

15418 FOR 1 = 1 TO N 
15428 PRINT USING fttt ; I,MASS(I) 

15421 LPRINT USING Hi* ; I, MASSU) 

15438 NEXT I 
15431 PRINT: PRINT* 



....PRESS Ctfl KEY TO CONTINUE 



15432 AS = INKEYS : IF AS = " THEN GOTO 15432 
15434 CLS 

15449 PRINT: PRINT: PRINT 

15458 PRINT' FREQUENCIES USED IN ANALYSIS <RPS)' -.PRINT 

15451 LPRINT 

15452 LPRINT'***** FREQUENCIES USED IN ANALYSIS (RPS) 



*****’ iLPRINT 



15468 

15478 

15471 

15488 

15481 

15482 
15484 
15498 
15491 
15588 
15581 
15518 
15528 
15538 
15548 
15558 
15551 
15568 
15579 
15571 
15588 
15598 

15591 

15592 

15593 

15594 
15688 
15618 
15628 
15621 
15638 
15648 
15641 
15658 
15668 

15661 

15662 

15663 
15678 
15688 
15698 
15691 
15788 
15718 



FOR 1 = 1 TO M 
PRINT I;' " ;U(I) 

LPRINT lj' 'iU(I) 

NEXT 1 

PRINT SPRINT’ ....PRESS ANY KEY TO CONTINUE' 

AS = INKEYS : IF AS = *' THEN GOTO 15482 
CLS 

PRINT:PR!NT:PRINT 
LPRINT iLPRINT 

PRINT'MODES USED IN ANALYSIS' :PR1NT 

LPRINT'***** MODES USED IN ANALYSIS ***»*' iLPRINT 
K = M : US = 'IIII.IHIIHHI AAAA ' 

IF K ) 5 THEN GOTO 15688 
FOR J = 1 TO N 

FOR I = 1 TO M 

PRINT USING US; PHI<1,J>; 

LPRINT USING Ul; PHKI.J); 

NEXT I 
PRINT ’ ’ 

LPRINT ' ' 

NEXT J 
PRINT: PRINT 
LPRINT : LPRINT 

PRINT: PRINT* ....PRESS ANY KEY TO CONTINUE' 

AS = INKEYS : IF AS = ” THEN GOTO 15593 
CLS: GOTO 15748 
FOR J = 1 TO N 

FOR I = 1 TO 5 

PRINT USING US; PHI(1,J); 

LPRINT USING US; PHKI.J); 

NEXT I 
PRINT ' ' 

LPRINT ' ' 

NEXT J 

PRINT :PRINT iLPRINT iLPRINT 

PRINT* . . . .PRESS ANY KEY TO CONTINUE' 

AS = INKEYS : IF AS = " THEN GOTO 15662 
CLS: PRINT: PRINT 
FOR J = 1 TO N 

FOR 1 = 6 TO N 

PRINT USING US; PHK1,J); 

LPRINT USING US; PHKI.J); 

NEXT I 
PRINT ' ' 
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15711 LPR1NT * ' 

15728 NEXT J 

15738 PRINT : PRINT :LPRIMT:LPR1NT 

15731 PRINT' . . . .PRESS ANY KET TO CONTINUE' 

15732 A* = INKEY* : IF A* = " THEN GOTO 15732 

15734 CIS: PRINT •. PRINT 

15735 CLOSE 8! 

15748 ■'************************************************************************ 
15758 ' NORMALIZE MODE SHAPES 

15768 ' ************************************************************************ 
15778 FOR I = 1 TO M 
15788 BIG = 8! 

15798 FOR J = 1 TO N 

15888 IF BIG < PHKI.J) THEN BIG = PHKI.J) 

15818 NEXT J 

15828 FOR J = 1 TO N 

15838 PHKI ,J) = PHI < I , J)/B3 G 

15848 NEXT J 

15858 NEXT I 

15868 ''************************************************************************ 
15878 •' OUTPUT NORMALIZED MODE SHAPES 

15888 '************************************************************************ 



15898 PRINT'NORMALIZED MODES USED IN ANALYSIS' .‘PRINT :LPRINT 

15891 LPRINT'***** NORMALIZED MODES FOR ANALYSIS *****' 

15892 LPRINT 
15988 K = M 

15918 IF K > 5 THEN GOTO 15988 

15928 FOR J = 1 TO N 

15938 FOR I = 1 TO M 

15948 PRINT USING W; PHI < I , J) ; 

15941 LPRINT USING U$;PHI< I, J); 

15958 NEXT 1 

15968 PRINT ' ' 

15961 LPRINT ' 1 

15978 NEXT J 

15971 PRINT: PRINT' 

15972 A* = INKEY* : IF A* 

15973 CLS 



15988 

15998 



....PRESS ANY KET TO CONTINUE’ 
" THEN GOTO 15972 

GOTO 16128 



LPRINT : LPRINT 
FOR J = 1 TO N 

FOR I = 1 TO 5 
16818 PRINT USING Uf; PHKI.J); 

16811 LPRINT USING U$; PHKI.J) J 

16828 NEXT I 

16838 PRINT ' ' 

16831 LPRINT ' ' 

16848 NEXT J 

16858 PRINT :PRINT :LPRINT :LPRINT 

16851 PRINT' . . . .PRESS ANY KEY TO CONTINUE' 

16852 A* = INKEYI : IF A* = " THEN GOTO 15852 

16853 CLS 
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14040 FOR J = 1 TO N 
14078 FOR I = 4 TO N 



14878 

14888 

148B1 

14090 

14188 

14181 



PRINT USING U»; PHKI ,J) ; 
LPRINT USING UJ; PHKI.J); 



NEXT I 
PRINT ' * 
LPRINT 1 ’ 



14110 NEX 
14128 PRI! 
14121 PRINT 



NEXT J 

PRINT :PRINT iLPRINT :LPRINT 



....PRESS ANY KEY TO CONTINUE" 



14122 A$ = INKEY$' : IF A* = " THEN GOTO 14122 
14138 CLS : GOSUB 4000 

14140 CLS : LOCATE 21,1,0 

14150 PRINT 

14155 CLS : GOSUB 8888 

14148 CLS : GOSUB 8200 

14178 CLS : GOSUB 10808 

20090 CLS 

20810 LOCATE 5,1,8 

20020 PRINT 1 *** YOU MUST FIRST USE THE PROGRAM JACOBI TO SOLVE EIGENVALUE ***" 
20030 PRINT"*** PROBLEM. IT HILL CALCULATE NATURAL FREQUENCIES AND MODE ***" 
20040 PRINT"*** SHAPES THAT CAN THEN BE USED TO RUN THE DOAN PROGRAM. ***" 

20050 PRINT: PRINT: PRINT: 

20040 PRINT" ENTER <F3> JACOB I <RETURT4> * :PRIKT 

28878 PRINT" <F2>" 

200 B0 STOP 

25000 CLSiPRINT"*** USER MAY USE OPTION 4 TO INSTALL HIS/HER SUBROUTINE ***" 
25018 PRINT" *** BEGINNING AT LINE 25000.***" 

25020 FOR I = 1 TO 5088 : NEXT I : GOTO 270 
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C. PROGRAM JACOBI - INSTRUCTIONS FOR USE 

The fol 1 ow i n 9 BAS' I C program may De ussd as- i s- or 
modified by the user to solve the Eigenvalue problem once 
the e q u a t i on s of mo t i on ar e f ormu 1 a t e d . It is self e x p 1 an - 
a t or y i n e x ecu t i on with the e x c e p t i on s that the A ma t r i x 
replaces the stiffness matrix and the B matrix replaces th 
mass matrix. The program only requires that the upper tr i 
angle elements of the stiffness matrix be suppl ied along 
with the diagonal elements of the mass matrix. 



D. PROGRAM JAC08I LISTING (Adapted Iren ret'. 13) 



4 CLS 

5 REM =CYCLIC JACOBI METHOO= 

19 DIM A(10,10), 8(10),U(10,10) 

1? REM =INPUT DATA= 

28 G0SU8 70 
28 REM=INITIALI2E= 

30 GQSUB 440 
40 GOSUB 980 

45 REM =0UTPUT RESULTS= 

50 GOSUB 750 
40 GOTO 3000 
45 REN = INPUT DATA= 

70 PRINT 

80 PRINT 1 CYCLIC JACOBI METHOD' 

148 PRINT 
180 PRINT 

190 PRINT ' ENTER THE SIZE OF PROBLEM (i.e. no. ot rows)' 
200 PRINT 
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218 INPUT 'N = ’ ; N 

228 LPRIMT "PROBLEM DIMENSION = N 

239 INPUT "NUMBER OF SIGNIFICANT FIGURES = B ; Si 

240 LPRIMT 
258 PRINT 

268 PRINT " ENTER THE UPPER TRIANGULAR ELEMENTS OF MATRIX A ( by coiurans ):" 
278 PRINT 

288 LPRIMT "UPPER ELEMENTS OF MATRIX A =" 

298 LPRIMT 

328 FOR J = 1 TO N 

338 PRINT 

348 PRINT " ENTER UPPER PART OF COLIM1 ’ j J 
358 PRINT 

368 FOR I = 1 TO J 

378 PRINT " A(";I;Y;J;’) = 

388 INPUT A< I , J) 

381 LPR1NT * A<";I;Y ;J;") = "; A<I,J) 

382 LPRIMT 
398 NEXT I 
488 PRINT 
448 PRINT 
458 NEXT J 

468 FOR I = 1 TO N 
478 II = I - 1 
488 FOR J = 1 TO II 
498 A(I,J) = A(J,D 
588 NEXT J 
518 NEXT I 
528 PRINT 

538 PRINT " ENTER ELEMENTS OF DIAGONAL MATRIX B:" 

531 LPRIMT " ELEMENTS OF DIAGONAL MATRIX B " 

532 LPRIMT 
548 PRINT 

578 FOR I = 1 TO N 
588 PRINT " B< " ; I ; " ) = *; 

598 INPUT 8(1) 

591 LPRIMT * B< " ; I ; ' ) = ";B(I> 

688 NEXT I 

610 PRINT 

611 LPRIMT 
658 RETURN 

655 REM INITIALIZE TOLERANCE AND MAX NO, OF ROTATIONS = 

668 Z = 2 * SI 
678 T1 = 1/(18 A Z) 

638 PRINT 

698 CLS: PRINT "TOLERANCE =" ;T1 

691 LPRIMT "TOLERANCE = " ;T1 

788 R = 5*N*N 

718 R1 = 8 

728 T2 = ,1 

738 N1 = N - 1 
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748 RETURN 

746 REM “OUTPUT El GENSOLUTI 0M== 

758 PRINT 
768 PRINT 

778 PRINT ■SOLUTION" 

771 LPRINT 
738 PRINT 
781 LPRINT 

732 LPRINT ***************** SOLUTION ****************** 
783 LPRINT 

888 PRINT * NUMBER OF ROTATIONS REQUIRED = *;R1 

381 LPRINT 'NUMBER OF ROTATIONS REQUIRED = *;R1 

382 LPRINT 
883 LPRINT 
818 PRINT 
328 PRINT 

838 FOR J = 1 TO N 
348 PRINT 

358 PRINT * EIGENVALUE *jJ;* IS *;B(J) 

351 LPRINT "EIGENVALUE *;Jj* IS *■ ;B<J) 

852 LPRINT 

853 LPRINT *1TS EIGENVECTOR IS * 

868 PRINT 

878 PRINT * ITS EIGENVECTOR IS * 

888 PRINT 

898 FOR 1 = 1 TO N 
988 PRINT U<I,J) 

981 LPRINT U<I,J) 

982 LPRINT 
918 NEXT 1 
928 PRINT 
921 LPRINT 

938 PRINT * PRESS SPACE 8AR TO CONTINUE * 

948 IF INKEY* = * * THEN 942 

941 GOTO 948 

942 CLS 
958 PRINT 
968 NEXT J 
978 RETURN 

975 REM “EIGENPROBLEM SOLUTJON= 

988 GOSUB 1130 

985 REM “PERFORM ONE CYCLE OF ROTATIONS 
998 GOSUB 1298 

995 REM = CHECK TOLERANCE = 

1888 IF XI < T1 THEN 1118 

1885 REM = CHECK NO. OF ROTATIONS = 

1818 IF R1 > R THEN 1848 
1828 T2 = .1 * XI 
1838 GOTO 998 
1848 PRINT 

1858 PRINT * *** ERROR **** 
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1851 LPRINT '**#***#**#♦**#**#*** ERROR ********************' 
1868 PRINT 

1878 PRINT ‘ NO CONVERGENCE ATTAINED’ 

1871 LPRINT ’ NO CONVERGENCE ATTAINED ’ 

1888 PRINT 
1881 LPRINT 

1898 PRINT ’ WITH ’;R1;’ ROTATIONS’ 

1891 LPRINT ’ WITH ‘ ;R1 ROTATIONS’ 

1188 END 
1118 G0SU8 1838 
1128 RETURN 
1138 FOR I = 1 TO N 
1148 FOR J = 1 TO N 
1158 U(I,J) = 8 
1168 NEXT J 
1178 U(I,I) = 1 
1188 NEXT I 
1198 FOR I = 1 TO N 
1288 B1 = SQR(Bd)) 

1218 8(1) = 1/81 

1228 NEXT I 

1238 FOR 1 = 1 TO N 

1248 FOR J = 1 TO N 

1258 Ad , J) = Bd) * Ad ,J) * 8(J) 

1268 NEXT J 
1278 NEXT 1 
1288 RETURN 
1298 XI = 8 
1388 FOR K = 1 TO N1 
1318 K1 = K + 1 
.1328 FOR L = K1 TO N 
1338 A1 = A(K‘K) 

1348 A2 = A(K,L) 

1358 A3 = A(L,L) 

1368 X = A2 * A2 /<A1 * A3) 

1365 REM = CHECK IF ROTATION IS NEEDED = 

1378 IF X ) XI THEN 1398 
1388 GOTO 1888 
1398 XI = X 

1488 IF X < T2 THEN 1888 

1418 R1 = R1 4 1 

1415 REM = COMPUTE ANGLE = 

1428 IF A1 = A3 THEN 1478 
1438 2 = .5 * (A1 - A3) /A2 
1448 21 = 1 4 1/ <Z * Z) 

1458 T = -Z * d + SQR(ZD) 

1468 GOTO 1488 
1478 T = 1 

1488 C = 1 / SQR< 1 4 T * T) 

1498 S = C * T 
1588 S2 = S * S 
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1518 C2 = C * C 
1528 A(K,L) = 8 

1525 RBI = TRANSFOrai DIAGONAL ELEMENTS = 

1538 A0 = 2 * A2 * C * S 

1548 A<K,K) = A1 *C2 ♦ A8 ♦ A3 * S2 

1558 A(L,L) = A1 * S2 - A8 + A3 * C2 

1555 REM = TRANSFORM OFF DIAGONAL ELEMENTS = 

1568 FOR 1 = 1 TO N 
1578 IF I < K THB1 1688 
1588 IF I ) K THB1 1648 
1598 GOTO 1748 
1688 A8 = A(I,K) 

1618 A(I,K) = C * A8 ♦ S * A<I,L) 

1628 A(I,L) = -S * A8 * C * A(I,L) 

1638 GOTO 1748 
1648 IF I < L THEN 1678 
1658 IF I > L THEN 1718 
1668 GOTO 1748 
1678 A8 = A(K,I) 

1688 ACK.l) = C * A0 + S * A(I,L) 

1698 A< I f L) = -S * A8 ♦ C * A(I,L) 

1788 GOTO 1748 
1718 A8 = A< K , I ) 

1728 A(K , I ) = C * A8 ♦ S * AIL, I) 

1738 A(L,I) = -S * A8 + C * A(L,I) 

1748 NEXT I 

1745 REM = TRANSFORM MATRIX U TO GB1ERATE E1GBWECTORS 
1758 FOR I = 1 TO N 
1768 U8 = U(I ,K) 

1778 U(I,K) = C * U8 + S * U(I,L) 

1788 U(I,L) = -S * U8 ♦ C * U<I ,L) 

1798 NEXT I 
1888 NEXT L 
1818 NEXT K 
1828 RETURN 

1825 REM = NORMALIZE EIGBWECTORS = 

1838 FOR I = 1 TO N 
1848 FOR J = 1 TO N 
1858 UU.J) = U<I,J) * 8(1) 

1868 NEXT J 
1878 NEXT 1 
1888 FOR I = 1 TO N 
1898 8(1) = A(I,I) 

1988 NEXT I 

1985 REM = ORDER EIGB1SOLUT1CN = 

1918 FOR 1 = 1 TO N1 
1928 II = I M 
1938 Z = 8(1) 

1948 M = I 

1958 FOR J = II TO N 
1968 IF Z ( B< J) THB1 1998 
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1978 2 = B(J) 

1938 hi = J 
1998 NEXT J 
2888 B<M) = B(I) 

2818 8U) = 2 
2828 FOR J = 1 TO N 
2838 Z = U(J,1) 

2848 U(J,1) = U<J,M) 

2858 U( J,M) = 2 
2869 NEXT J 
2878 NEXT I 
2888 RETURN 

3888 PRINT 1 END OF PROGRAM 1 
3818 END 
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APPENDIX S: PROGRAMS USED TO CONVERT VELOCITY HISTORIES 



TO DEFLECTION HISTORIES 



18 

20 ' PROGRAM ** S1MPS1 ** USED FOR CASE 1 

30 ' INTEGRATES VELOCITY HISTORY WRT TIME TO OBTAIN DEFLECTIONS 

40 

50 CLS : DIM A*<3) 

S3 F(l)=8 : F(2)=0 : F(3)=0 : F0RHH=" 8. 81188 ’ : F0RM2*=" 88.888888 " 



70 INPUT" ENTER THE NAME OF AN INPUT FILE ",FILE1N$ : PRINT 

80 INPUT' ENTER THE ORDINATES TO INTEGRATE OVER ", NUMBER : PRINT 

90 INPUT" ENTER THE SIZE OF THE INTEGRATION INTERVAL ",H 

100 CLS 



110 OPEN FI LEIN* FOR INPUT AS 81 : OPEN "D1SPR2.0UT" FOR OUTPUT AS 82 

128 ' 

130 ' USE SIMPSON'S 1/3 RULE TO CALCULATE DEFLECTIONS 

140 ' 

150 INPUT 81 ,M,TIME<8) ,VEL11(0) ,VEL12<8) ,VEL13(8) ,VEL61(8) ,VEL62<0) ,VEL63<0) 

160 FOR I = 1 TO 2 

170 INPUT 81 ,M ,T1ME< I ) ,VEL1 1(1) ,VEL1 2( I ) ,VEL13( I ) ,VEL61 ( I ) ,VEL62( I ) ,VEL63( I ) 
180 NEXT I 

190 VIA = VEL6K8) - VEL1K0) 

200 V1B = VEL6K1) - VELIKI) 

210 VIC = VEL6K2) - VEL1K2) 

228 91 = (VIA ♦ 4*V1B + VIC) * H / 3 
230 V2A = VEL62(8) - VEL12(8) 

240 V2B = VEL62(1) - VEL12<1) 

258 V2C = VEL62(2) - VEL12(2) 

260 92 = (V2A ♦ 4*V2B + V2C) * H / 3 
270 V3A = VEL63(0) - VEL13(0) 

2B8 V3B = VEL63(1) - VEL13(1) 

290 V3C = VEL63(2) - VEL13(2) 

300 93 = (V3A + 4*V3B + V3C) * H / 3 

318 F(l) = F(l) ♦ 91 

320 F(2) = F(2) + 92 

330 F( 3) = F(3) + 93 

348 T = TIME(l) * 1800 

350 FOR 1 = 1 TO 3 

360 IF ABS(QMAXd)) > ABS(F(D) THEN GOTO 380 
370 QHAXd) = Fd) : TmXd) = T 

380 NEXT I 

390 PRINT 82, " T = PRINT 82, USING FORMMjT; 

408 PRINT 82, "msec DX = "j ; PRINT 82, USING F0RH2t;F(l); 

410 PRINT 82, "in. DY = PRINT 82, USING F0RM2$;F(2); 

420 PRINT 82, "in. DZ = *;: PRINT 82, USING F0RM21;F(3);:PRINT 82, "in." 

438 VEL1K0) = VEL1K2) i VEL12(0) = VEL12(2) : VEL13(0) = VEL13(2) 
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440 VEL6K8) = VEL6K2) : VEL62<8) = VEL62-2) : VEL630) = VEL63<2) 

450 IF M ♦ 2 < NUMBER THEN GOTO 1 60 

4(50 At(l) = ’(X) = * : Af(2) = ’(Y) = 1 : Af(3) = ’(2) = 8 

470 FOR I = 1 TO 3 

480 PRINT 82, ’TNAX‘;A*(I); : PRINT #2, USING F0RH1$;TMAX(D ; 

485 PRINT 82,* usee DMAX" ;A*<I) ; : PRINT 42, USING F0RH2t;Qmx<I) ; 

486 PRINT 02,* in. 1 
490 NEXT I 

583 CLOSE 01 : CLOSE 02 



10 ' 

20 ' PROGRAM ** SINPS2 ** USED FOR CASES II AND III 

30 ' INTEGRATES VELOCITY HISTORY URT TIME TO OBTAIN DEFLECTIONS 

46 — 

50 CLS : DIM A*<3) 

60 F<1)=8 : F<2)=8 : FORMIC’ 00,0084 * : FORM2i=“ 08.088888 * 



70 INPUT’ ENTER THE NAME OF AN INPUT FILE ’ ,FILEIN$ : PRINT 

80 INPUT’ ENTER THE ORDINATES TO INTEGRATE OVER ’, NUMBER : PRINT 

90 INPUT’ ENTER THE SIZE OF THE INTEGRATION INTERVAL ’ ,H 

100 CLS 



110 OPEN FI LEIN* FOR INPUT AS 81 : OPEN ’DISPR2.0UT’ FOR OUTPUT AS 82 
115 PRINT 82, ’H = ’;H;’ SEC’ 

120 ' 

130 ' USE SIMPSON'S 1/3 RULE TO CALCULATE DEFLECTIONS 

140 ' - 

158 INPUT 81 ,M,TIME(0) ,VEL12(0) ,VEL92(0) ,VEL52(0) ,VEL112<8) ,VEL192<0) ,VEL152(0) 

168 FOR I = 1 TO 2 

170 INPUT 81 ,M ,T IME( I ) ,VEL 12( I ) , VEL92C I ) ,VEL52< I ) ,VEL1 12(1) ,VEL192( I ) ,VEL152< I ) 
180 NEXT I 

190 VIA = (VEL12(0) ♦ VEL92(0))/2 - VEL52<0) 

280 V1B = (VEL12U) ♦ VEL92(l)>/2 - VEL52U) 

210 VIC = (VEL12C2) ♦ VEL92<2))/2 - VEL52<2) 

220 Q1 = (VIA ♦ 4*V1B 4 VIC) * H / 3 

230 V2A = (VEL1 12(8) + VEL192(0))/2 - VEL152(0) 

240 V2B = (VEL1 1 2< 1 ) ♦ VEL192(l))/2 - VEL152(1) 

250 V2C = (VEL112(2) ♦ VEL192(2))/2 - VEL152(2) 

260 Q2 = (V2A 4 4*V2B 4 V2C) * H / 3 

318 F(l) = F(l) 4 oi 

320 F<2) = F(2) ♦ Q2 

340 T = TIME(l) * 1000 

350 FOR 1 = 1 TO 2 

360 IF ABS(QMAXd)) ) A8S(F(D) THEN GOTO 380 
370 QMAX(I) = F(I) : TMAX(I) = T 

380 NEXT I 

390 PRINT 82, ’ T = PRINT 82, USING FORMAT; 

400 PRINT 82, ’msec D(5,Y) = *; : PRINT 82, USING F0RM2*;F<1); 
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418 PRINT 42, "in. 0(15, Y) = * j: PRINT 42, USING F0fti2-i;F(2) ; 

420 PRINT 42, "in." 

430 UEL12(0) = VEL12<2) : VEL92(8> = VEL92<2> : VEL52<8> = VEL52<2) 

448 VEL112<0> = VELU2<2) : VEL192<0> = VEL192(2> : VEL152C0) = VEL152<2) 

459 IF M + 3 < NUMBER THEN GOTO 140 

440 A*<1) = " < 5 ,Y) = " : Ai(2) = "(15, Y) = " 

478 FOR 1 = 1 TO 2 

430 PRINT 42, "THAX" ;Ai(I) ; : PRINT 42, USING F0RM1$;TMAX<I); 

485 PRINT 42," msec DMAX" ;AI(I) ; : PRINT 42, USING FGRM2$;(MX(1); 

434 PRINT 42," in." 

498 NEXT I 

588 CLOSE 41 : CLOSE 42 
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APPENDIX C ~ TYPICAL CASE I ANALYSIS 

As noted in Chapter V , the first step in each analys 
was the construction of an appropriate SAP I V model of th 
equipment substructure under consideration. In this cas 
the model represents a 1,000 lb value attached to the tr- 
end of a cantilever f oundat i on . 

A. SAP IV INPUT DATA FOR CASE I 

The following data file was used to generate the mod 
information for the Case I model [Ref. 10]. 



HEADING CARD 

CASE 1 — VALUE ON CANTILEVER FOUNDATION 
MASTER CONTROL CARD 

7 . 1 . 0 . 4 . 1 . 0 . 0 . 0 . 0.0 
NODAL FT DATA 

0,1 ,1 , 0 , 0,1 , 1 , 1 , 0 . , 0 . , 0 . , 1 , 0 .@ 
0,2,1,0,0,0,1,1,0.,0.,5.,1,0.0 
0,6,1,0,0,0,1,1,0.,0.,25.,1,0.0 
0,7,1 ,1 ,1,1 ,1 ,1 ,10. ,0. ,10.,0.0 

TYPE 2 - 3D BEAM ELEMENTS 

2 , 5 , 1 , 0 , 1 

MATERIAL PROPERTIES 

1 , 30 . 0ES , 0 . 3 , 7 . 29E-4 ,0.0 
ELEMENT PROPERTY CARD 

1.17.64.0. 0.0.0. 539 . 95 . 333 . 95 . 20 6 . 0 



0 . 0 , 0 . 0 , 0 . 0 , 0 . 0 
0 . 0 ,0.0,0. 0 , 0 . 0 
0 . 0 , 0 . 0 , 0 . 0 ,0.0 
0.0, 0.0, 0.0, 0.0 

BEAM DATA CARDS 

1 .1 .2.7.1 .1 .0. 0.0.0.000000.1 

5 . 5 . 6 . 7 . 1 ,1 ,0,0,0,8,008000,1 

CONCENTRATED MASS ON FREE END 

6.0. 2.591 ,2.591 , 2 . 59 1 , 8 . 8 , 8 . 8 , 8 . 8 

8.0. 0. 0 0 0 , 0 . 0 0 0,0.000,0.0,0.0,0. 0 

DYNAM I C ANA L Y S I S C A R D < I N ACCOR DAN C E WITH R E F . 4 

0 , 0 , 0 , 0 . 0 , 0 . 0 , 8 , 0 . 0 
1 ,8 

1 . 1 . 1 . 1 . 1 . 1 . 1 

END OF DATA (DUMMY TITLE CARD) 

0 , 0 , 0 , 8 , 8 , 0 , 0 , 0 , 8 , 0 

The output -from the SAP IV code is not included here but 
all the essential in-formation is included in the USLOB code 
output . 

B. FT CRUST INPUT DATA FOR CASE I 

The PI CRUST code is used to process the SAP IV output 
and calculate influence coefficients corresponding to the 
base or support motions of the substructure. Additionally, 
c on n e c t i v i t y be twe en the hull an d f ou n da t i on is specified 
and a substructure file is produced. 

The following input file was used for CASE I: 

LIST OPTIONS 

1 , 1 , 1 , 1 , 8 , 1 , 1 , 1 , 8 , 8 
0 , 1 , 1 , 1 , 1 , 1 , 0 , 1 , 0 , 0 

M 0 DA L S P E C. I F I CAT I ON 
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5.0. 0 

CONNECTIVITY 

1,5,10 

270 

MODES FROM SAP IV INCLUDED IN ANALYSIS 

1 , 2 , 3 , 4 , 5 

END 

The output 1 i sting -from the PI CRUST code is not 
included here but all essential output -from these 
c a 1 c u 1 a t i on s is included in the USLOB ou t p u t . 

C. USLOB INPUT DATA FOR CASE I (1414” STANDOFF) 

The -following input data was suppl i ed to the USLOB 
code to specify the parameters for the time integration 
and associated output. A total of 1201 time steps are 
specified f or the i n t e gr a t i on inter v a 1 , e ac h c ov e r i n g an 
interval of 5 . 6 E - 6 seconds. The XXXX parameters denote 
parameters for eq. 5 which are omm i tted due to the 
sensitivity of this e q u a t i on . 

The following input data correspond to the USLOB output 
included in this appendix. 

TIME STEP AND CONTROL CARD FOR INTEGRATION 

1201 , 1 , 1 , 9 , 11 , 1,1 

EQUATION (5) PARAMETER SPECIFICATION 

5 . 0 E-6 , XXXX ,1606,0.0 

5000 . , XXXX ,XXXX ,XXXX ,XXXX 
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OUTPUT OPTIONS 



0,1 , 0 , 0 , 0, 0 , 0 , 0,1 ,0 
0 , 0,1 ,1 , 0 , 0 , 0 , 1 , 0,0 

SHELL VELOCITY OUTPUTS AT EACH TIME STEP 
1 , 5,10 

SUBSTRUCTURE VELOCITY OUTPUTS AT EACH TIME STEP 

6 

1 3 ^ 3 ^ 3^-3^ , S , i* , 1 , S ,<£,&, S 

CTh i s n o t a t i on specifies the x , y , z t r an s 1 a t i on s o+ t h 
base an d tip n ode s of" the su bs tructure re sp e c t i m e 1 y . ) 



END 
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D. LISTING - USIOB OUTPUT FOR CASE I 



CASE 1 -- VALVE CM FOUNDATION VERT OISPL 



CONTROL PARAMETERS FOR SUBSTRUCTURE 1 



NPOINT = 7 
NCOOF = 3 
NFREQ = 10 
NSYMTY = 9 



NDOF = 18 
NUDOF = 15 
NFACE = I 
NFORCE = 13 



CONNECT I CM DATA FOR INTERFACE POINTS 

SHL SUB X Y 2 XX 

56 I 1 1 1 0 

0DEGROT = 9.BBBB08E+81 

1 EQUATION NUMBERS FOR SUBSTRUCTURE 1 



YY 22 LS6PT LFIND DEGCON 
0 0 5010 27 2.700000E+B2 



POINT X Y 

1 1 2 

2 4 5 

3 7 8 

4 10 11 

5 13 14 

6 16 17 

7 0 0 



2 XX YY 22 

3 0 8 0 

6 0 0 0 

9 0 0 0 

12 0 0 0 

15 0 0 0 

18 0 0 0 

0 8 0 0 



FIXED-BASE NATURAL FREQUENCIES FOR SUBSTRUCTURE 1 

MODE FREQ MODE FREQ MODE FREQ MODE FREQ 

1 4.45541E+82 2 1.05485E+03 3 1.34387E+83 4 4.03785E+03 

5 7.60871E+03 6 9.55805E+03 7 1.84549E+04 3 1.21696E+84 

9 1.22840E+04 18 1 .89142E+04 



1 COMPACTED FIXED-BASE MODES FOR SUBSTRUCTURE 1 



1 2 

1 1.9139E-18 -1.2402E-91 

2 3.1568E-28 8.9711E-89 

3 -1.2482E-81 -4.8588E-29 

4 3.8186E-18 -2.4745E-81 

5 6.2987E-29 1.7908E-08 

6 -2.4745E-81 -9.6785E-20 

7 5.7852E-18 -3.6971E-01 

8 9.4185E-28 2.6743E-88 



3 4 

9.5490E-18 6.2397E-12 
1.2402E-81 1 . 1 183E-12 
3.7187E-11 -1.4719E+88 
1.9853E-89 1.8036E-11 
2.4745E-01 1 . 8086E-1 2 
5.9782E-11 -2.3685E+90 
2.8465E-09 9.9880E-12 
3.6971E-01 1.7811E-12 



5 

2.4775E-11 

1.6989E-13 

-2.3785E+08 

3.9865E-11 

2.7333E-13 

-1.4497E+00 

3.9373E-11 

2.6992E-13 
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9 -3.6971E-81 -1 . 4460E-1 ? 
18 7.5646E-18 -4.9020E-61 
11 1 . 2478E-1 9 3.5459E-08 



5.8957E-11 -2.3393E400 
3 . 7743E-89 5.9093E-12 
4.9928E-01 1.0677E-12 



1 .4839E+83 
2.3492E-1 1 
1.6103E-13 



12 -4.9828E-01 -1.9173E-19 3.5172E-11 -1.3957E+80 2.3572E+90 

13 9.3830E-13 -6.3836E-01 4.6843E-09 -3.9555E-13 -1 .5718E-12 

14 1.5485E-19 4.4096E-08 6.0836E-01 -5.9949E-14 -1 .8775E-14 

15 -6.8336E-81 -2.3794E-19 -2.3536E-12 9.3383E-02 -4.23&2E-92 



1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

13 

14 

15 

16 



1 .4719E+88 
1.3628E-13 
2.4183E-10 
2.3685E+00 
2.1925E-13 
-1.5016E-10 
2.3393E+80 
2.1655E-13 
-1 .4869E-10 
1 . 3957E+90 
1.2929E-13 
2.4276E-10 
-9.3383E-82 
-8.6446E-15 
-2.2919E-12 



7 

3.5420E-14 
2.2993E-11 
-2.3685E+80 
5.3840E-14 
3.6999E-11 
1 .4719E+08 
4 . 9358E-1 4 
3.6543E-11 
1 .4538E+80 
2.7551E-14 
2.1883E-11 
-2.3754E+00 
-1.9932E-15 
-1 . 4588E-1 2 
2.2420E-02 



8 

1.2133E-20 
-1 .4719E+80 
-3.3664E-12 
1.5886E-20 
-2.3685E+00 
5.4501E-12 
1.1252E-28 
-2 . 3393E+00 
-5.4572E-12 
4.4918E-21 
-1 .3957E+00 
3.3850E-12 
-4.7535E-22 
9.3383E-02 
-2.3083E-14 



? 

1.2095E-12 
6.6535E-13 
-1 . 4634E+ 80 
-7.6500E-13 
1.0680E-12 
2.3693E+00 
-4.6836E-12 
1.0516E-12 
-2 . 3724E+ 00 
-6.5906E-12 
6.2577E-13 
1.4715E+00 
-5.2983E-12 
-4.1996E-14 
-1.0035E-02 



18 

-2 . 3705E+88 
9.7987E-11 
-1 . 8812E-12 
-1 . 4497E+ 00 
5.7223E-11 
-3.7449E-12 
1 . 4839E+00 
-6.8350E-11 
-5.6047E-12 
2.3572E+00 
-1.0703E-10 
-7.4239E-12 
-4.2362E-02 
-7.7491E-12 
-9.2163E-12 



1 SU8STRUCTURE CONTROL PARAMETERS 

SU8 NFREQ NOOF NCDOF NPTS LOCO NEQ8 N8L0CK NROUFC 

1 10 18 3 7 518 15 1 3 

N8LFQR KOL8AR N8L8AR 8EGFRQ 8EGS8L 8EGBAR 8EGF0R 
1 254 1 1 1 1 1 

CONTROL PARAMETERS FOR SOLUTION 



NTIME = 


1201 


NEQS = 


638 


NSK1P = 


2 


NRECS = 


681 


NCHRG = 


1 


NTHETA = 


17 


NQUAD = 


9 


NDELTA = 


21 


NF1NE = 


11 


NQDOT = 


264 


KOUPLE = 


2 


NSU8S = 


1 



SOLUTION PARAMETERS 



TOTAU1 = 3.541492E+04 


AREA = 


3.639044E+06 


RHOFLU = 9.452160E-05 


CFLU = 


6.008800E+04 


DELT = 5.080000E-86 


TEND = 


6.000000E-83 


XLOAD = 8.844880E+82 


STOFF = 


1 .41 4000E+ 03 


RLOAD = 1 . 696080E+ 03 


ORNRAD = 


0 . 
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DATA FOR EMPIRICAL PRESSURE-TIME HISTORY 



WCHRG = 5.388888E+83 
PZMLT = XXXXXXJOOOO (X 

pzexp = xxxxxxxxxxxx 



SURCUT = XXXXXXXXXXXX 
TH1LT = XXXXXXXXXXXX 
THEY? = XXJCCYXXXKXXX 



1VEL0C1TY DISTRIBUTION AT SHELL STATION 56 (SEGMENT 5 POINT 18) 



'EP 


TIME 


UDOT < 0) 


UDOT( 90) 


UDQTC188) 


UDOT (278) 


1 


8. 


0. 


0. 


0. 


8. 


3 


1.8888E-85 


4.9444E-81 


7.5446E-92 


1 . 7639E-8 1 


7.5445E-82 


5 


2.8888E-85 


1 . 8855E+08 


2.8722E-01 


6.5644E-81 


2.8720E-81 


7 


3.8888E-85 


3.8391E+00 


5.3428E-81 


1 . 3085E+88 


5.8418E-01 


9 


4.8888E-95 


6 . 1 484E+88 


9.3465E-01 


2.8219E+88 


9.3436E-01 


11 


5.8888E-85 


8.6881E+88 


1.3193E+00 


2 . 7669E+88 


1 .3136E+88 


13 


6.8888E-05 


1 . 1 347E+ 8 1 


1.7214E+00 


3.4983E+00 


1 . 7202E+00 


15 


7.8888E-85 


1.4145E+81 


2.1416E+00 


4 . 2275E+ 08 


2 . 1397E+08 


17 


8.0808E-05 


1 . 7083E+8 1 


2.5792E+08 


4.9559E+88 


2.5763E+88 


19 


9.8880E-85 


2.8192E+01 


3.8265E+88 


5. 6665E+88 


3 . 0224E+08 


21 


1.8080E-04 


2.3241E+01 


3.4895E+08 


6.3713E+00 


3 . 4848E+88 


23 


1.1008E-84 


2.6547E+81 


3.9697E+08 


7.0994E+00 


3.9625E+08 


25 


1.2008E-04 


2.9933E+81 


4.4593E+88 


7.8237E+88 


4.4582E+88 


27 


1.3808E-84 


3.3481E+81 


4.9651E+88 


8.5656E+08 


4.9539E+80 


29 


1 . 4888E-04 


3.7123E+81 


5.4786E+88 


9.3141E+88 


5 . 4658E+88 


31 


1.5008E-84 


4.0922E+01 


6.8854E+00 


1 .0891E+81 


5.9891E+88 


33 


1.6000E-04 


4.4771E+81 


6.5359E+88 


1 .8856E+0 1 


6.5167E+88 


35 


1.7008E-04 


4.3745E+81 


7.0787E+00 


1 . 1 647E+0 1 


7.8482E+88 


37 


1.3800E-04 


5.2721E+81 


7.6848E+88 


1 . 241 1E+81 


7 . 5788E+88 


39 


1.9088E-84 


5. 6819E+0 1 


8 . 1 394E+98 


1.3280E+81 


8.1097E+88 


41 


2.0088E-04 


6.8891E+8 1 


8.6690E+00 


1 . 3957E+8 1 


8 . 6353E+ 00 


43 


2.1888E-84 


6.5836E+81 


9.1899E+08 


1 . 4723E+0 1 


9.1519E+88 


45 


2.2008E-84 


6.9161E+81 


9.7044E+00 


1 .5459E+0 1 


9 . 6618E+00 


47 


2.3888E-94 


7.3292E+81 


1 . 8209E+81 


1 .6185E+81 


1 .3161E+01 


49 


2 . 4009E-04 


7 .7437E+8 1 


1 .0700E+0 1 


1 . 6899E+8 1 


1 . 8648E+8 1 


51 


2.5888E-04 


3. 1493E+81 


1 .1180E+81 


1.7569E+01 


1.11 22E+0 1 


53 


2.6000E-84 


8 . 5686E+01 


1 . 1 643E+0 1 


1 . 8247E+8 1 


1 . 1 580E+8 1 


55 


2.7889E-04 


8.9672E+81 


1 . 2097E+0 1 


1 . 8895E+0 1 


1 . 20 27E+0 1 


57 


2.8888E-84 


9.3669E+81 


1 .2528E+0 1 


1.9509E+81 


1 . 2451E+8 1 


59 


2 . 9099E-94 


9.7634E+81 


1 . 291 8E+ 0 1 


2 . 8 1 86E+0 1 


1 . 2835E+8 1 


61 


3.0080E-04 


i . 8 154E+82 


1.3301E+01 


2.0669E+81 


1 . 3209E+ 0 1 


63 


3.1008E-04 


1 .0542E+82 


1 . 3669E+81 


2.1219E+01 


1 . 3578E+0 1 


65 


3.2880E-84 


1 . 8933E+82 


1 . 4026E+0 1 


2.1774E+81 


1 .3918E+01 


67 


3.3988E-84 


1.131 7E+92 


1 . 4379E+9 1 


2.2388E+81 


1 .4261E+81 


69 


3.4888E-84 


1.1781E+02 


1 . 4712E+01 


2.2828E+81 


1 . 4585E+ 0 1 


71 


3.5880E-84 


1 . 2986E+82 


1 . 5834E+ 0 1 


2 . 3353E+0 1 


1 . 4896E+0 1 


73 


3.6980E-84 


1 . 2465E+8 2 


1 . 5354E+0 1 


2.3858E+01 


1 . 528 5E+8 1 


75 


3.7888E-04 


1 . 2846E+82 


1 . 5655E+ 8 1 


2 . 4379E+0 1 


1 .5494E+81 


77 


3.8080E-04 


1 . 3226E+ 82 


1.5948E*81 


2.4892E+81 


1 . 5767E+8 1 



99 



79 

81 

83 

85 

37 

8 ? 

91 

93 

95 

97 

99 

181 

183 

105 

107 

109 

111 

113 

115 

117 

119 

121 

123 

125 

127 

129 

131 

133 

135 

137 

139 

141 

143 

145 

147 

149 

151 

153 

155 

157 

159 

181 

163 

165 

167 

169 

171 

173 

175 

177 

179 



3.9800E-04 


1.3599E+82 


1 . 6222E+8 1 


2.5382E+01 


1 . 6836E+0 1 


4.0033E-04 


1 . 3974E+02 


1 . 6498E+0 1 


2.5834E+81 


1 . 6291 E+8 1 


4.1080E-04 


1.4349E+82 


1.6748E+81 


2.6395E+01 


1.6526E+81 


4.2300E-04 


1 . 471 9E+82 


1 .6937E+01 


2.6885E+81 


1.6759E+01 


4.3008E-04 


1 . 5088E+02 


1.7222E+81 


2.7371E+01 


1 . 6979E+0 1 


4.4080E-04 


1 . 5454E+82 


1.7432E+01 


2.7366E+81 


1.7174E+01 


4.5080E-84 


1.5817E+82 


1.7639E+81 


2.3350E+01 


1.7365E+01 


4.6000E-04 


1 .6175E+82 


1 .7843E+9 1 


2.3822E+01 


1 . 7553E+8 1 


4.7000E-04 


1.6533E+02 


1.8017E+81 


2.9304E+01 


1 .7710E+01 


4.3808E-04 


1 . 6888E+0 2 


1 .8180E+01 


2.9778E+01 


1 . 7858E+0 1 


4.9888E-04 


1.7237E+82 


1.8342E+01 


3.0233E+81 


1 . 8060E+0 1 


5.0800E-84 


1.7588E+82 


1.3487E+81 


3.8677E+81 


1.8127E+01 


5.1008E-84 


1.7928E+82 


1.8681E+81 


3.1125E+01 


1.8223E+01 


5.200BE-04 


1 . 8255E+02 


1.8714E+01 


3.1550E481 


1 . 331 7E+8 1 


5.3088E-84 


1 . 8585E+82 


1.8824E+81 


3.1964E+01 


1.8488E+01 


5.4008E-04 


1.8910E+82 


1 . 8980E+0 1 


3.2371E+81 


1.8463E+01 


5.5008E-84 


1.9229E+02 


1 .8957E+01 


3.2765E+01 


1.8501E+81 


5.6800E-84 


1 .9548E+82 


1.9013E+01 


3 .31 36E+8 1 


1 . 3537E+0 1 


5.7000E-64 


1.9845E+82 


1 . 9863E+81 


3.3489E+01 


1 . 8565E+0 1 


5.8088E-04 


2.8147E+82 


1 .9871E+81 


3.3847E+81 


1 . 8552E+31 


5.9000E-04 


2.0442E+82 


1 . 9074E+81 


3.4182E+01 


1 .8533E+81 


6.0080E-84 


2.0738E+82 


1 . 9876E+8 1 


3.4496E+01 


1 .8512E+81 


6.180BE-84 


2.1010E+82 


1.9057E+01 


3.4789E+01 


1.8471E+01 


6.2880E-04 


2 . 1 283E+82 


1 . 8995E+8 1 


3.5874E+81 


1 . 8388E+0 1 


6.3880E-04 


2 . 1 550E+02 


1.8932E+01 


3.5335E+91 


1.3388E+01 


6.4008E-84 


2.1811E+02 


1 . 8869E+0 1 


3 .5583E+8 1 


1.8213E+01 


6.5000E-04 


2.2067E+02 


1 .3798E+81 


3 . 5823E+8 1 


1 .81 10E+81 


6.6888E-84 


2.2317E+82 


1 . 8681 E+0 1 


3.6856E+81 


1 . 7978E+8 1 


6.7000E-04 


2.2583E+02 


1 . 8572E+0 1 


3.6277E+01 


1.7844E+31 


6.8888E-84 


2.2885E+82 


1 .8463E+81 


3.6484E401 


1 .7711 E+8 1 


6.9008E-04 


2.3041E+82 


1.8340E+01 


3.6682E+01 


-1 . 7564E+01 


7.0000E-04 


2.3271E+82 


1.8189E+01 


3.6873E+01 


1 . 7389E+0 1 


7.1000E-04 


2.3497E+02 


1.8835E+01 


3.7054E+01 


1 . 721 1E+01 


7.2000E-84 


2.3719E+02 


1 . 7382E+8 1 


3.7222E+01 


1 . 7034E+0 1 


7.3008E-04 


2.3936E+82 


1.7723E+01 


3.7379E+81 


1.6851E+01 


7.4888E-04 


2.4147E+82 


1 . 7531E+8 1 


3.7534E+01 


1 .6635E+01 


7.5000E-84 


2.4354E+02 


1 . 7336E+0 1 


3.7878E+01 


1.6415E+81 


7.6888E-84 


2.4558E+82 


1 . 71 46E+8 1 


3.7806E+01 


1 .6201E+81 


7.7000E-04 


2.4757E+02 


1 .6961E+01 


3.792 4E+01 


1.5993E+81 


7.3880E-84 


2.4958E+82 


1.6736E+81 


3.8842E+01 


1.5744E+01 


7.9808E-04 


2.5139E+82 


1 . 6504E+0 1 


3.8147E+01 


1.5489E+01 


8.0000E-84 


2.5325E+02 


1 . 6279E+9 1 


3.8242E+01 


1.5240E+01 


8.1888E-84 


2.5508E+02 


1 . 6862E+9 1 


3.8326E+01 


1.5001E+01 


8.2008E-84 


2 . 5685E+82 


1.5809E+01 


3.8487E+81 


1 . 4725E+01 


8.3800E-84 


2.5857E+82 


1.5546E+81 


3.8481E+01 


1 .4448E+01 


8.4888E-84 


2.6826E+82 


1 . 5298E+8 1 


3.8542E+01 


1.4161E+01 


8.5888E-84 


2.6193E+02 


1.5039E+01 


3.8594E+01 


1 . 3883E+01 


8.6800E-84 


2.6355E+82 


1.4778E+81 


3.8640E+81 


1 . 3606E+0 1 


8.7088E-84 


2.6510E+82 


1.4488E+81 


3.8680E+01 


1.3294E+01 


8.8800E-04 


2.6663E+02 


1 . 4288E+8 1 


3.8718E+81 


1 . 2985E+8 1 


8.9888E-04 


2.6813E+82 


1 .391 7E+0 1 


3.8730E+01 


1.2682E+81 



180 



131 

183 

185 

187 

189 

191 

193 

195 

197 

199 

201 

283 

285 

287 

209 

211 

213 

215 

217 

219 

221 

223 

225 

227 

229 

231 

233 

235 

237 

239 

241 

243 

245 

247 

249 

251 

253 

255 

257 

259 

241 

243 

245 

247 

249 

271 

273 

275 

277 

279 

281 



9.0088E-04 
9.1800E-84 
9.2080E-04 
9.3839E-04 
9.4000E-04 
9.5008E-04 
9.4000E-04 
9.7000E-84 
9.8000E-04 
9.9808E-84 
1.6000E-03 
1.8180E-03 
1 . 0200E-03 
1 . 0300E-03 
1.0480E-03 
1.8508E-83 
1.8488E-03 
1.8780E-03 
1.0800E-03 
1.8988E-83 
1.1889E-83 
1.1188E-83 
1.1 280E-83 
1.1380E-83 
1.1480E-83 
1 . 1 588E-83 
1.1488E-83 
1 . 1 700E-03 
1.1 800E-03 
1 . 1988E-03 
1.2080E-03 
1 . 210&E-03 
1.2288E-83 
1.2308E-83 
1.2488E-03 
1.2508E-03 
1.2400E-83 
1.2780E-83 
1.2888E-83 
1 . 2908E-83 
1.3888E-83 
1.3100E-03 
1.3200E-03 
1.3388E-83 
1 . 3480E-83 
1.3588E-83 
1.3400E-03 
1.3788E-83 
1 . 3300E-03 
1.3988E-33 
1 . 4000E-83 



2.4941E+02 
2 . 71 82E+02 
2.7233E+02 
2.7372E+82 
2.7585E+82 
2 . 7433E+ 02 
2.7755E+02 
2.7873E+82 
2.7991E+82 
2.3104E+82 
2 .821 4E+82 
2.8320E+82 
2.8423E+82 
2 . 8524E+02 
2 .8423E+0 2 
2.8715E+02 
2.8892E+02 
2 . 8887E+ 02 
2.8971E+02 
2 .9054E+02 
2 . 91 30E+82 
2.9281E+02 
2.9272E+02 
2.9342E+02 
2.9411E+82 
2 . 9473E+82 
2.9531E+82 
2.9589E+82 
2.9444E+02 
2 . 9703E+02 
2.9752E+82 
2.9794E+02 
2.9841E+82 
2 . 9884E+ 02 
2.9929E+02 
2 . 9944E+82 
2.9999E+02 
3.8031E+82 
3 . 0843E+8 2 
3.8094E+02 
3.8121E+02 
3.0142E+02 
3.8144E+02 
3.81 85E+02 
3.0204E+02 
3.0225E+02 
3.0234E+02 
3.8244E+82 
3. 8254E+02 
3.0247E+92 
3.0278E+02 



1 .3443E+01 
1 . 3335E+0 1 
1 . 38 19E+0 1 
1 . 2789E+8 1 
1.2484E+81 
1.2180E+81 
1.1745E+81 
1 . 1 429E+01 
1.1108E+01 
1 . 0773E+8 1 
1 . 0442E+01 
1 .8889E+8 1 
9.7338E+98 
9 . 3938E+98 
9.0543E+80 
8.78 44E+99 
8.3494E+99 
7 .9933E+08 
7.4442E+00 
7 . 39 34E+99 
4.9413E+80 
4.5447E+09 
4 . 1 998E+99 
5.841 4E+00 
5 . 4887E+90 
5 . 1 148E+09 
4.7374E+80 
4.3449E+99 
4 . 0038E+90 
3 . 4475E+99 
3.2711E+80 
2 . 891 1E+00 
2.5193E+00 
2 . 1 528E+88 
1 .7894E+88 
1 .4148E+00 
1 .9388E+90 
4.4715E-81 
2.9858E-01 
-4.5745E-02 
-4.3242E-01 
-8.0484E-81 
-1 . 1 744E+00 
-1 .5432E+00 
-1 .9944E+09 
-2. 2444E+08 
-2.4320E+00 
-2 . 9942E+88 
-3.3579E+09 
-3.7125E+00 
-4.0445E+00 



3.8739E+01 
3.3744E+01 
3.8744E+01 
3.8734E+81 
3.8713E+01 
3. 8484E+8 1 
3.9452E*81 
3.8488E+81 
3.8552E+01 
3.8489E+01 
3.8421E+01 
3.3344E+01 
3.8257E+81 
3.8143E+91 
3.8841E+81 
3.7956E+81 
3 . 7849E+0 1 
3.7734E+01 
3.7611 E+0 1 
3.7488E+01 
3.7345E+81 
3.7284E+01 
3.7858E+01 
3.4898E+01 
3.4729E+91 
3.4555E+81 
3.4375E+01 
3.4185E+81 
3.5985E+01 
3.5777E+01 
3.5571 E+9 1 
3.5359E+91 
3.5138E+01 
3.4989E+01 
3.4474E+81 
3.4442E+0 1 
3.4285E+01 
3. 3959E+ 0 1 
3.3707E+81 
3.3451E+01 
3.3198E+01 
3.2925E+01 
3.2447E+81 
3.2343E+01 
3.2870E+01 
3.1 774E+0 1 
3.1478E+01 
3.11 75E+0 1 
3.8849E+01 
3 .8553E+8 1 
3.0233E+01 



1 . 2383E+0 1 
1 . 2860E+0 1 
1 . 1725E+01 
1 . 1 397E*8 1 
1 . 1 077E+8 I 
1 .8754E+0 1 
1 . 0403E+8 1 
1 . 8852E+81 
9.7097E+08 
9.3739E+00 
9 . 0253E+00 
3.4484E+00 
8 . 2991 E+ 00 
7.9448E+00 
7 . 5959E+00 
7.2392E+08 
4.8749E+80 
4.5125E+00 
4.1599E+00 
5.8128E+88 
5.4448E+00 
5.8494E+88 
4.7009E+00 
4.3419E+00 
3.9894E+00 
3.4193E+08 
3.2424E+00 
2.8742E+00 
2.5183E+00 
2.1484E+88 
1 .7998E+00 
1 . 4288E+80 
1 .0473E+00 
7 . 1 232E-0 1 
3.4192E-01 
3.2835E-83 
-3.5943E-81 
-7.1455E-01 
-1.0453E+00 
-1 . 4187E+00 
-1.7571E+08 
-2.1101E+00 
-2.4572E+80 
-2.8000E+00 
-3.1344E+00 
-3 . 4702E+88 
-3.8104E+00 
-4 . 1 442E+00 
-4.4775E+00 
-4.8823E+80 
-5 . 1 225E+00 



181 



283 

285 

287 

289 

291 

293 

295 

297 

299 

301 



1.4189E-03 
1.4288E-83 
1.4308E-83 
1.4438E-03 
1.4500E-03 
1.4488E-83 
1.4700E-03 
1.4800E-93 
1 . 4900E-83 
1.5000E-03 



3 . 8280E+02 
3.0230E+02 
3.8279E+82 
3.8280E+02 
3.8281E+82 
3.0277E+82 
3.8247E+82 
3 , 0257E+ 82 
3.8247E+82 
3 , 0239E+ 02 



-4.4214E+00 
-4.7754E+88 
-5.1242E+80 
-5.4737E+88 
-5.81 80E+80 
-4.14UE+88 
-<4.581 8E+ 80 
-4.8409E+80 
-7 . 1 776E+89 
-7.5121E+08 



2.9917E+01 
2 . 9599E+3 1 
2 . 9276E+ 9 1 
2.8944E+01 
2.3404E+8I 
2 . 3263E+8 1 
2.7934E+81 
2.7599E+81 
2 . 7255E+8 1 
2.4983E+01 



-5.44-43E+08 
-5 . 76<42E+08 
-<4 .0313E+98 
-4.3920E+08 
-<4 . <4982E+8 0 
-7.8821E+08 
-7.3924E+00 
-7.5999E+08 
-7.3941E+00 
-8 . 1 850E+99 



( This history is abbreviated tor illustrative purposes.) 
IVELOCITf DISTRIBUTION ALONG SUBSTRUCTURE 1 



STEP 
1 

3 1. 
5 
7 
9 



TIME VEL< 1,1) VEL( 1,2) 



7.4928E-04 -4.4088E-83 
3.4147E-03 -2.8700E-02 
8.8040E-03 -3 . 3246E-0 2 
1 . 6429E-02 -3 . 7045E-82 
11 5.0000E-05 -2 . 6541 E-02 -2. 68<4<4E-02 



0080E-05 

2.8088E-05 

3. 

4. 



0808E-05 

0080E-85 



13 


4.8008E-85 -3 . 8728E-82 


3.8282E-04 


15 


7.0000E-05 -5, 2889E-02 


4.4051E-82 


17 


8.0000E-05 -4 , 881 9E-82 


1 . 0709E-0 1 


19 


9 . 8008E-05 -8.4191E-82 


1 . 3835E-01 


21 


1.0000E-84 -1 , 8522E-0 1 


2.9143E-01 


23 


1 . 1 909E-84 -1 ,251 4E-81 


4.1250E-01 


25 


1 . 2000E-94 -1 , 4552E-0 1 


5.5122E-81 


27 


1.3000E-04 -1.4745E-01 


7.1054E-01 


29 


1.4000E-04 -1.9027E-81 


8.8589E-01 


31 


1 . 5000E-84 -2.1381E-01 


1.0795E+00 


33 


1.4008E-84 -2.3783E-01 


1 . 2921E+08 


35 


1 , 7888E-04 -2.4205E-81 


1 . 51 73E+08 


37 


1 , 88 00E-04 -2.8485E-01 


1 . 7424E+00 


39 


1.9088E-04 -3.1149E-0I 


2 . 8 222E+8 0 


41 


2.0808E-84 -3.3481E-01 


2.2970E+00 


43 


2.1000E-04 -3, 41 29E-81 


2.5879E+00 


45 


2.2008E-04 -3.8400E-81 


2 . 8949E+00 


47 


2 , 3000E-04 -4.1029E-01 


3.2141E+00 


49 


2.4808E-04 -4.3351E-01 


3.5481E+08 


51 


2.5000E-04 -4.5582E-01 


3.8942E+00 


53 


2.4000E-04 -4 , 7747E-8 1 


4 . 2550E+08 


55 


2.7000E-04 -4.9885E-01 


4 . 4245E+88 


57 


2.8000E-04 -5.1332E-01 


5.0121E+00 


59 


2.9800E-84 -5.3428E-81 


5.411 7E+ 08 


41 


3.0000E-04 -5.5399E-01 


5.8235E+00 


43 


3.1008E-04 -5.7835E-01 


4 . 2448 E+ 00 


45 


3.2000E-04 -5.8550E-81 


4.4781E+88 


47 


3 , 3008E-04 -5.9947E-01 


7.1191E+00 


49 


3.4090E-04 -4 . 1 224E-0 1 


7.5494E+00 



VEL( 

0. 

7.5445E-82 
2.8720E-01 
5.3410E-81 
9.3434E-01 
1.3184E+00 
1 . 7282E+00 
2.1397E+00 
2.5743E+88 
3.0224E+00 
3.4343E+88 
3.9425E+88 
4.4582E+80 
4.9539E+88 
5 . 4650E+8 8 
5.9891 E+88 
4.5147E+88 
7 . 8482E+88 
7.5788E+03 
8.1697E+00 
8.4353E+88 
9,151 9E+88 
9,661 8E+ 38 
1 . 0 1 61 E+ 01 
1 .0648E+81 
1 . 1 1 22E+ 01 
1 . 1 588E+ 81 
1 . 2027E+8 1 
1 .2451E+81 
1 .2835E+81 
1 . 3209E+81 
1 . 3570E+01 
1.3918E+81 
1 .4261E+01 
1.4585E+81 



1,3) VEL( 4,1) VEL< 4,2) 



4.4872E-87 
1 , 9928E-84 
1 . 4589E-0 4 
-3.5152E-08 
-4.3943E-04 
-2 . 8539E-85 
-1 . 21 47E-04 
-3.7008E-84 
-8 , 8429E-04 
-1.7341E-03 
-3.1440E-83 
-5.1171E-03 
-7,71 18E-83 
-1 .181 7E-82 
-1 , 5095E-02 
-1.9994E-82 
-2.5775E-02 
-3,251 4E-82 
-4.8329E-02 
-4.9349E-02 
-5.9491E-02 
-7.1427E-02 
-8.4401E-02 
-9.9240E-02 
-1 , 1534E-0 1 
-1 , 3297E-8 1 
-1 . 5208E-8 1 
-1.7244E-81 
-1 , 9473E-0 1 
-2.1829E-01 
-2 . 4339E-0 1 
-2.7002E-81 
-2.9812E-01 
-3 . 2748E-8 1 



-1.9745E-05 
-3 , 3974E-85 
1 . 2345E-85 
4.4999E-05 
-1 , 4797E-84 
-9 . 1 328E-84 
-2.3448E-83 
-4.2451E-83 
-4 , 0383E-03 
-4.8775E-03 
-5 . 5941 E-93 
-1.2343E-83 
7.0835E-03 
2.8128E-02 
3.8738E-82 
4.3995E-02 
9.7299E-02 
1 , 4034E-0 1 
1 . 9492E-8 1 
2.4245E-01 
3.4499E-01 
4.4315E-81 
5.5824E-01 
4.9154E-01 
8.4451E-01 
1,81 38E+08 
1 . 21 41E+ 00 
1 . 4379E+ 00 
1 . 4854E+ 88 
1.9593E+83 
2.2403E+80 
2.5890E+00 
2.9441E+08 
3.3327E+80 



VEL< 4,3) 

8. 

-4.3583E-14 
1.4181E-12 
3 . 81 35E-1 8 
1.3005E-88 
1 . 8231 E-07 
1.5878E-04 
8.4295E-04 
3.7798E-05 
1 . 3459E-04 
4.8597E-04 
1 , 9480 E-83 
2.5841E-83 
5.3195E-03 
1.0372E-82 
1 , 8757E-02 
3 , 1 738E-02 
5.0420E-02 
7.4405E-82 
1 , 1 045E-0 1 
1.5337E-81 
2.9584E-01 
2 . 4573E-0 1 
3.3527E-01 
4.1357E-01 
5.0041E-01 
5.9453E-81 
7.0155E-81 
8.1594E-01 
9.3985E-01 
1 . 8733E+00 
1 . 21 43E+ 88 
1 . 3485E+00 
1.5299E+00 
1 , 700 1 E+ 00 
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71 

73 

75 

77 

79 

81 

83 

85 

87 

89 

91 

93 

95 

97 

99 

181 

193 

105 

107 

109 

111 

113 

115 

117 

119 

121 

123 

125 

127 

129 

131 

133 

135 

137 

139 

141 

143 

145 

147 

149 

151 

153 

155 

157 

159 

161 

163 

165 

167 

169 

171 



3.5088E-84 
3.6000E-04 
3.7000E-84 
3.3000E-04 
3.9009E-84 
4.0000E-84 
4.1080E-04 
4.2000E-04 
4.3080E-04 
4.4000E-04 
4.5000E-84 
4.6000E-04 
4.7000E-04 
4.8008E-04 
4.9800E-04 
5.0008E-04 
5.1888E-04 
5.2088E-04 
5.3008E-04 
5.4000E-04 
5.5800E-04 
5.6000E-04 
5.7000E-04 
5.8008E-04 
5.9880E-84 
6.8008E-84 
6.1000E-84 
6.2000E-04 
6.3080E-84 
6.4088E-04 
6.5008E-84 
6.6008E-04 
6.7800E-04 
6.8000E-04 
6.9088E-04 
7.0000E-04 
7.1880E-04 
7.2000E-04 
7.3080E-04 
7.4080E-04 
7.5080E-04 
7.6000E-04 
7.7080E-04 
7.8880E-04 
7.9000E-04 
8.8000E-04 
8, 1000E-04 
8.2000E-04 
8.3000E-04 
8.4000E-04 
8.5000E-04 



-6.2460E-8! 
-6.3518E-01 
-6.4496E-01 
-6.5416E-81 
-6.6229E-01 
-6 . 7034E-81 
-6.7756E-81 
-6.8457E-01 
-6.9093E-01 
-6.9747E-01 
-7.8376E-01 
-7.1020E-81 
-7 . 1 622E-8 1 
-7.2251E-01 
-7.2391E-01 
-7.3577E-01 
-7 . 4298E-0 1 
-7.5894E-81 
-7.5834E-01 
-7.6640E-01 
-7.7556E-01 
-7.8511 E-81 
-7.9475E-01 
-8.0498E-01 
-8.1562E-81 
-8.2643E-01 
-8.3804E-01 
-8.5826E-01 
-3.6285E-01 
-8.7428E-01 
-8.8701E-01 
-8.9966E-01 
-9 . 1 225E-0 1 
-9.2519E-81 
-9.3799E-81 
-9.5874E-01 
-9 . 6374E-0 1 
-9 . 7703E-01 
-9.3938E-01 
-1.8821E+00 
-1 . 0 1 58E+00 
-1 . 0277E+08 
-1 .8396E+00 
-1 .8519E+00 
-1.0647E+00 
-1 . 0765E+00 
-1.0881E+00 
-1 . 1 002E+08 
-1.U21E+00 
-1.1 235E+00 
-1 . 1348E+00 



3.0300E-88 
8.4943E+80 
8.9636E+83 
9.4589E+00 
9.9375E+00 
1 .8431E+01 
1 . 9931E+0 1 
1 . 3 437E+0 1 
1 . 1948E+0 1 
1 . 2465E+0 1 
1.2989E+01 
1 . 351 2E+8 1 
1 . 4043E+01 
1 . 4538E+9 1 
1 .51 18E+01 
1 . 5659E+0 1 
1 . 6208E+01 
1 .6760E+81 
1 .7311E+81 
1 .7871E+01 
1 . 8433E+8 1 
1 .8995E+01 
1.9559E+01 
2.0129E+81 
2.0792E+81 
2.1275E+01 
2.1851E+01 
2.2436E+81 
2.3023E+01 
2.3608E+01 
2.4194E+81 
2.4788E+01 
2.5380E+81 
2.5971E+01 
2 . 6562E+8 1 
2 . 7162E+01 
2.7768E+01 
2.8354E+01 
2.8950E+01 
2.9552E+81 
3.0155E+81 
3.0755E+01 
3.1352E+81 
3.1957E+0I 
3.2563E+01 
3.3167E+01 
3.3766E+01 
3.4372E+01 
3.4980E+01 
3.5586E+01 
3.6189E+01 



1.4896E+01 
1 .5285E+0 1 
1 .5494E+8! 
1 .5767E+01 
1 . 6036E+8 1 
1 .6291E+01 
1 . 6526E+81 
1 . 6759E+8 1 
1.6979E+01 
1.7174E+81 
1 . 7365E+81 
1 . 7553E+0 1 
1.771 0E+91 
1.7356E+01 
1 . 8808E+81 
1 . 81 27E+0 1 
1 . 8223E+01 
1.3317E+81 
1 . 8408E+3 1 
1 . 8463E+0 1 
1 . 8501E+81 
1 . 8537E+8 1 
1.8565E+01 
1.8552E+81 
1.8533E+81 
1 .851 2E+81 
1 . 8471 E+01 
1 .8386E+31 
1 . 3300E+ 01 
1 . 821 3E+81 
1.311 8E+01 
1 . 7978E+81 
1 . 7844E+0 1 
1 .771 1E+81 
1 .7564E+01 
1 .7389E+81 
1.721 1E+0 1 
1 . 7834E+0 1 
1.6851E+01 
1.6635E+81 
1 . 6415E+01 
1 . 6281 E+81 
1 . 5993E+0 1 
1 . 5744E+81 
1 . 5489E+0 1 
1 .5240E+01 
1 . 500 1E+0 1 
1 . 4725E+81 
1 . 4448E+81 
1 .4161E+81 
1 . 3888E+0 1 



-3 . 5835E-0 1 
-3.9026E-81 
-4.2323E-81 
-4.5715E-81 
-4.9189E-81 
-5.2736E-81 
-5.6347E-31 
-6.891 3E-01 
-6.3723E-81 
-6.7464E-01 
-7.1217E-81 
-7.4965E-01 
-7 . 3691E-8 1 
-8.2381E-01 
-3.6322E-01 
-8.9682E-81 
-9.3109E-01 
-9.6536E-81 
-9 . 9372E-8 1 
-1 . 031 1 E+00 
-1 . 0623E+88 
-1 . 3923E+00 
-1 . 1299E+00 
-1 .1480E+88 
-1.1 735E+ 09 
-1 .1974E+00 
-1 . 21 97E+09 
-1 .2403E+80 
-1 . 2592E+00 
-1.2764E+08 
-1 .2919E+80 
-1 .3057E+88 
-1 .31 76E+80 
-1 . 3277E+80 
-1 . 3360E+00 
-1 . 3425E+ 00 
-1 . 3473E+00 
-1 . 3584E+88 
-1 . 3519E+80 
-1 .3518E+88 
-1 .3502E+00 
- 1 . 3472E+88 
-1 . 3429E+80 
-1 .3371E+88 
-1 .3301E+00 
-1 .321 9E+88 
-1 . 31 25E+00 
-1 .3021E+90 
-1 . 2909E+80 
-1.2789E+88 
-1 . 2664E+00 



3 . 7494E+ 98 
4.1969E+08 
4.6753E+00 
5 . I 844E+80 
5 . 7236E+88 
6.2922E+08 
6.8896E+00 
7.5155E+80 
8.1693E+08 
8.8583E+00 
9 . 5574E+08 
1 . 8289E+8 1 
1 . 1 044E+ 9 1 
1 . 1819E+81 
1 . 261 4E+ 01 
1 . 3426E+0 1 
1 . 4253E+8 1 
1 . 5895E+0 1 
1 .5950E+01 
1 . 681 4E+01 
1.7686E+31 
1 . 8562E+0 1 
1 . 9442E+9 1 
2 . 8321E+01 
2. 1 199E+81 
2 . 20 74E+0 1 
2.2944E+01 
2.3886E+81 
2.465?E*31 
2.5499E+81 
2.6325E+81 
2.7136E+01 
2.7930E+01 
2. 3786E+01 
2.9464E+01 
3.0201E+01 
3.0916E+81 
3.1608E+01 
3.2277E+01 
3.2922E+81 
3.3543E+01 
3 . 41 40E+8 1 
3.4713E+01 
3.5263E+0! 
3.5739E+01 
3.6292E+01 
3 . 6772E+8 1 
3.7231E+81 
3.7669E+81 
3.8088E+01 
3.8490E+01 



1.879SE+03 
2 . 0679E+88 
2.2639E+98 
2 . 4708E+08 
2.6857E+80 
2 . 9i 1 2E+88 
3. 1470E+00 
3.3937E+88 
3.6513E+89 
3.9219E+00 
4.2843E+80 
4.4993E+89 
4.8869E+08 
5.1266E+00 
5 . 4588E+ 00 
5.8802E+80 
6.1527E+03 
6.5146E*08 
6.8856E+88 
7.2652E+00 
7 . 6533E-* 83 
8.6496E+08 
3.4533E+88 
8.8657E+00 
9 . 2848E+08 
9.7105E*08 
1 . 8 1 42E+8 1 
1 . 0579E+8 1 
I . 1821E+0 J 
1 . 1 468E+8 1 
1 . 1 91 9E+0 1 
1 .2375E+01 
1 . 2836E+9 1 
1 . 3382E+81 
1.3774E+01 
1 . 4250E+0 1 
1 . 4738E+0 1 
1 . 521 5E+81 
1 . 5703E+8 1 
1 .6193E+01 
1 . 6685E+0 1 
1 . 71 77E+0 1 
1 . 7669E+0 1 
1 . 81 59E+8 1 
1 . 3647E+9 1 
1 .9131E+81 
1 .961 1E+01 
2.0935E+81 
2 . 0555E+01 
2. 1819E+81 
2.1477E+81 
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173 

175 

177 

17? 

181 

133 

135 

187 

189 

191 

193 

195 

197 

19? 

281 

203 

285 

287 

28? 

211 

213 

215 

217 

219 

221 

223 

225 

22 ? 

229 

231 

233 

235 

237 

239 

241 

243 

245 

247 

24? 

251 

253 

255 

257 

259 

261 

263 

245 

247 

24? 

271 

273 



8.4808E-04 -l.!4<53E+38 
8.7888E-84 -1.1572E+88 
8.8888E-84 -1.1479E+08 
S.9888E-84 -1.1785E+08 
9.8008E-04 -1.1392E+80 
9.1888E-84 -1.1993E+88 
9.2888E-84 -1 .2890E+08 
9.3888E-04 -1 .2193E+88 
9.4888E-84 -1 .2291E+88 
9.5888E-84 -1 .2383E+88 
9.<5888E-84 -1.2473E+88 
9.7888E-84 -1 .2564E+08 
9.8088E-84 -1 , 2655E+09 
9.9880E-84 -1 .2737E+08 
1.8888E-83 -1.2818E+88 
1.8188E-83 -1 . 2398E+88 
1.8280E-83 -1.2932E+88 
1.0388E-83 -1 .3857E+08 
1.8488E-83 -1.3129E+88 
1.8508E-83 -1 . 3288E+88 
1 .8488E-83 -1.3272E+8B 
1.8788E-83 -1 -3338E+88 
1.0808E-83 -1.3399E+88 
1.8988E-03 -1 .34<58E+88 
1.1888E-83 -1.3528E+88 
1.1188E-83 -1.3579E+88 
1.1208E-03 -1 .3(533E+08 
1 . 1 388E-83 -1 .3483E+88 
1.1488E-83 -1.3730E+88 
1.1588E-83 -1.3777E+88 
1.1688E-83 -1 . 3828E+ 88 
1.1788E-83 -1 .3848E+88 
1.1888E-83 -1.3988E+8B 
1 . 1 988E-83 -1 .3936E+88 
1.2888E-83 -1.3949E+88 
1.2188E-83 -1 .4888E+88 
1.2288E-83 -1.4038E+88 
1.2388E-83 -1.4858E+88 
1.2488E-83 -1 .4888E+88 
1.2588E-83 -1.4188E+88 
1 . 2608E-83 -1 .41 19E+88 
1.2788E-83 -1 -4137E+88 
1.2888E-83 -1.4151E+88 
1.2908E-83 -1.4148E+88 
1 .3888E-83 -! .41(5<5E+88 
1.3188E-83 -1 .4174E+88 
1.3288E-83 -1.4178E+88 
1.3388E-83 -1.4179E+88 
1.3488E-83 -1 .4175E+88 
1.3588E-83 -1.4148E+88 
1.36B8E-83 -1.4143E+88 



3.4798E+81 1 .3686E+81 
3.7488E+81 1.3294E+81 
3.8888E+81 1.2985E+81 
3.8613E+81 1 .2632E+81 
3.9213E+81 1 .2388E+8! 
3.9821E+81 1.2048E+81 
4 .0438E+8 1 1 . 1725E+81 
4.1838E+81 1 . 1397E+3 1 
4.1441E+81 1.1077E+81 
4.2248E+81 1.8754E+81 
4.2849E+81 1 .8483E+81 
4 . 3457E+81 1.B852E+81 
4.4842E+81 9.7897E+88 
4.46<51E+81 9.3739E+88 
4.5243E+81 9.8253E+88 
4.5878E+81 8.<S684E+88 
4.4477E*81 8.2991E+88 
4.7879E+81 7.9440E+08 
4.7474E+81 7.5959E+88 
4 . 8279E+8 1 7.2392E+88 
4.8885E+81 6.8749E+08 
4.9498E+81 4.5125E+88 
5 . 8889E+8 1 (5.1599E+88 
5.8484E+81 5.8128E+08 
5.1 285E+81 5.4443E+88 
5.1 889E+8 1 5.8494E+88 
5.2490E+81 4.7089E+08 
5 . 3885E+8 1 4.3419E+08 
5.347(5E+81 3.9394E+88 
5.4274E+81 3.4193E+88 
5.4874E+81 3.2426E+88 
5.5473E+81 2.8742E+88 
5.(58(5(5E+81 2.5133E+08 
5.4454E+01 2.1<584E+88 
5.7248E+81 1 .7998E+88 
5.7844E+81 1.4288E+88 
5.8437E+81 1 .8673E+88 
5.9825E+81 7.1232E-81 
5.9(58(5E+81 3.4192E-81 
4.8192E+81 3.2835E-83 
(5 .8781E+81 -3.5943E-81 
4.1348E+81 -7.1455E-81 
(5.1958E+81 -1 . 0653E+88 
6 . 2524E+8 1 -1 .4187E+88 
(5.3181E+81 -1 . 757 1E+88 
4.3484E+81 -2.U81E+88 
4.4248E+81 -2.4572E+88 
(5.484(5E+81 -2.8888E+88 
4.5419E+81 -3.1364E+88 
(5.5998E+81 -3.4782E+88 
(5 .(55(58E+81 -3.8184E+88 



-1 . 2532E+88 3.8875E+81 
-1 . 2397E+88 3.9247E+81 
-1 . 2257E+88 3.9685E+81 
-1 .211 5E+88 3.9951E+81 
-1 .1971E+80 4.8287E+81 
-1 . 1 826E+88 4.8415E+81 
-1 . ! 681E+80 4 . 8937E+8 1 
-1 . 1539E+98 4.1257E+81 
-1 . 1 399E+08 4.1575E+81 
-1 . 1264E+ 88 4 . 1 894E+81 
-1.11 35E+88 4.2216E+81 
-1.101 2E+ 80 4.2542E+01 
-1 ,8895E*88 4.2873E+81 
-1 . 8786E+08 4.3213E+81 
-1.8<53<5E+88 4.3563E+81 
-1 . 8594E+88 4.3926E+81 
-1 .8512E+88 4.4383E+81 
-1 . 8441 E+83 4 . 4696E+81 
-1 .0381E+88 4.5185E+81 
-1 . 8334E+88 4.5533E+01 
-1 . 8388E+88 4.5988E+81 
-1 . 8288E+88 4.4447E+8! 
-1 .0273E+88 4.4934E+81 
-1 .8288E+88 4.7447E+81 
-1 .838 1E+88 4 . 7982E+8 1 
-1 . 8335E+ 88 4.8541E+01 
-1 . 8384E+88 4.9124E+81 
-1 .8446E+88 4.9731E+81 
-1.8522E+08 5.0348E+01 
-1.8413E+88 5.1014E+81 
-1 . 8719E+08 5. 1 o98E+81 
-1 . 8838E+08 5.2389E+81 
-1.3971E+88 5.31 11E+81 
-1 .11 18E+88 5.3853E+81 
-1 . 1276E+88 5.4417E+81 
-1.1 445E+88 5.5398E+81 
-1 .1 626E+ 88 5. 6197E+81 
-1 .1817E+88 5.7812E+81 
-1 .291 7E+80 5 .7848E+0 1 
-1 . 2227E+88 5 . 3(581 E+81 
-1 . 2445E+88 5.9533E+81 
-1 . 2<S71 E+80 6 . 0395E+0 1 
-1 .2985E+88 (5.1244E+81 
-1.3144E+08 (5.2138E+81 
-1 .3387E+88 (5.3814E+81 
-1 .3433E+88 4.3898E481 
-1 . 3882E+00 (5.4745E+81 
-1.4131E+88 4.5634E+81 
-1 . 4388E+88 6.6581E+81 
-1 .4429E+88 (5.7360E+01 
-1 .4875E+88 4.8289E+81 



2.1929E+81 
2 . 2374E+8 1 
2.2311E+81 
2 . 3248E+81 
2.3<5<53E+0! 
2 . 4078E+01 
2.4471E+81 
2 . 4862E+81 
2 .5245E+0 1 
2.5617E+81 
2 . 5988E+81 
2 . 6333E+8 1 
2.6(57sE+01 
2.7887E+81 
2.7325E+81 
2.7(529E+01 
2.7918E+01 
2.8192E+81 
2.8448E+81 
2.8<587E+81 
2.8988E+81 
2.9111E+81 
2.9294E+81 
2.9458E+81 
2.9682E+01 
2.9727E+01 
2 . 9831 E+81 
2.9916E+81 
2 . 9988E+8 1 
3 . 8824E+81 
3 . 8848E+8 1 
3. 8858E+81 
3.8831E+81 
2 . 9992E+81 
2.9931E+81 
2 . 9858E+81 
2 . 9749E+01 
2.9427E+81 
2.9434E+81 
2 . 9328E+81 
2 . 91 34E+81 
2 . 8925E+81 
2.8492E+81 
2.8435E+81 
2.8153E+81 
2.7845E+01 
2.7512E+81 
2.7154E+01 
2.4772E+81 
2.6365E+9 1 
2.5933E+81 
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275 1.3700E-93 
277 1.3800E-03 
279 1.3980E-83 
281 1.4888E-83 
233 1 . 41 00E-83 
285 1.420BE-03 
287 1 . 4300E-03 
289 1 . 4488E-83 
291 1.4500E-33 
293 1.4600E-03 
295 1 . 4730E-03 
297 1.4800E-03 
299 1.4980E-03 
301 1.5888E-83 

(This history 



-1.4156E+B8 
-1 .4142E+88 
-1.4128E+00 
-1 . 41 09E+03 
-1.4091E+0B 
-1.4878E+00 
-1 , 4044E+08 
-1 .4818E+88 
-1.3987E+08 
-1 . 3956E+00 
-1.3921E+88 
-1.3884E+80 
-1.3844E+80 
-1.3881E+80 



6.7143E+01 
6. 7718E+01 
6 . S285E+8 1 
6 . 8845E+0 1 
4.9415E+81 
8.9984E+81 
7.8549E+81 
7.111 8E+81 
7,16o4E-*01 
7.2220E+01 
7.2780E+81 
7.3338E+01 
7.3392E+81 
7.4440E+81 



-4.1462E+88 
-4.4775E+08 
-4.8023E+88 
-5. 1225E+00 
-5.4463E400 
-5.7662E+00 
-8,0813E+08 
-4.3928E+80 
-6 . 6982E+88 
-7.0021E+80 
-7.3824E+00 
-7 .5999E+80 
-7.3941E+88 
-8.1850E+00 



-1 . 51 20E+88 
-1 . 5360E+ 80 
-1.5593E+00 
-1 .5825E+88 
-1 .6046E+88 
-1 . 6259E+00 
-1.6461E+88 
-1.M52E+88 
-1 . 6832E+ 08 
-1.8999E+80 
-1.7154E+08 
-1.7295E+00 
-1.7421E+00 
-1 .7533E+08 



6.9847E+01 
6 . 9870E+6 1 
7 . 0678E+8 1 
7.1469E+81 
7 . 224 1 E+0 1 
7 . 2993E+0 1 
7.3724E+81 
7.4433E+01 
7.51 18E+0 1 
7.5779E+01 
7.6414E+01 
7 . 7023E+0 1 
7 . 7605E+0 1 
7.8161 E+0 1 



is abbreviated Tor illustrative purposes.) 



2.5478E+91 
2.5808E+01 
2.4498E+81 
2.3973E+01 
2 . 3426E+ 0 1 
2 . 2856E+0 1 
2.2264E+81 
2 . 1 652E+9 1 
2.1018E+0! 
2 . 8363E+0 1 
1.9688E+81 
1 .8994E+01 
1*. 8280 E+8 1 
1 .7546E+81 
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E. INTEGRATION USING SIMPS1 

The -final step in the ELSHQK analysis of Case I is to 
u se an ap p r op r i ate i n t e g r a t i on a 1 gor i t hm t o c a 1 cul a t e a 
de-flection history -from the velocity i n-f ormat i on provided 
by the US LOB c ode . S i mp son *'* s 1 /3 rule was c h ose n i n t h i s 
work and in this example, implemented in the -form o-f 
S I MRS 1. BAS which appears in appendix B. 

The foil ow i n g is an entire ou t p u t listing f r om t h i s 
program. In the listing, the coordinate directions, are 



defined as : 



X 

Y 

7 



fore and aft 

vertical or normal to the hull 
athwar t sh i ps 



T = 


8.0198 


msec 


ox = 


0.000000 


in. 


DY = 


0.806666 


in. 


DZ = 


-8.308802 


T = 


0.0308 


msec 


DX = 


9.099999 


in. 


DY = 


0.900601 


in. 


DZ = 


-0.000014 


T = 


8.0580 


msec 


DX = 


9.000001 


in. 


DY = 


' 3.000901 


in. 


DZ = 


-0.308848 


T = 


8.0708 


msec 


DX = 


0.909002 


in. 


DY = 


8.990000 


in. 


DZ = 


-B.088383 


T = 


0.8980 


msec 


DX = 


0.960094 


in. 


DY = 


-8.808094 


in. 


DZ = 


-0.883144 


T = 


0.1100 


msec 


DX = 


9.990996 


in. 


DY = 


-0.900012 


in. 


DZ = 


-0.008223 


T = 


8.1380 


msec 


DX = 


0.000099 


in. 


DY = 


-0,808926 


in. 


DZ = 


-8.000322 


T = 


0.1580 


msec 


DX = 


9.099013 


in. 


DY = 


-0.000847 


in. 


DZ = 


-8.888441 


T = 


0.1700 


msec 


DX = 


0.999918 


in. 


DY = 


-0.089075 


in. 


DZ = 


-8.000581 


T = 


8.1988 


msec 


DX = 


0.090923 


in. 


DY = 


-0 .8081 12 


in. 


DZ = 


-9.000741 


T = 


0.2100 


msec 


DX = 


8.000029 


in. 


DY = 


-0.008157 


in. 


DZ = 


-0.800920 


T = 


0.2308 


msec 


DX = 


9.900036 


in. 


DY = 


-0.080210 


in. 


DZ = 


-0.881116 


7 = 


8.2588 


msec 


DX = 


0.090943 


in. 


DY = 


-8.888271 


in. 


DZ = 


-0.081329 


T — 
1 - 


8.2708 


msec 


DX = 


0.000050 


in, 


DY = 


-0.088339 


in. 


DZ = 


-0.881555 


T = 


0.2988 


msec 


DX = 


9.999956 


in. 


DY = 


-0.088414 


in. 


DZ = 


-0.301793 


T = 


0.3100 


msec 


DX = 


0.090963 


in. 


DY = 


-0.880493 


in. 


DZ = 


-0 .002048 


T = 


8.3300 


msec 


DX = 


9.909969 


in. 


DY = 


-0.888577 


in. 


DZ = 


-0.082295 


T = 


0.3508 


msec 


DX = 


0.908074 


in. 


DY = 


-0.808662 


in. 


DZ = 


-8 .882555 


T = 


8.3700 


msec 


DX = 


8.080879 


in. 


DY = 


-0.088748 


in. 


DZ = 


-0 .002819 


T = 


8.3988 


msec 


DX = 


0.000082 


in. 


DY = 


-0.080832 


in. 


DZ = 


-0.083886 


7 = 


0.4109 


msec 


DX = 


0.800084 


in. 


DY = 


-0.088913 


in. 


DZ = 


-0.083354 


7 = 


0.4388 


msec 


DX = 


9.909085 


in, 


DY = 


-0.000988 


in. 


DZ = 


-0.003628 


7 = 


0.4590 


msec 


DX = 


0.000085 


in. 


DY = 


-0.001057 


in. 


DZ = 


-0.803883 


7 = 


0.4780 


msec 


DX = 


8.006684 


in. 


DY = 


-0.001117 


in. 


DZ = 


-0.004141 


7 = 


8.4908 


msec 


DX = 


0.086881 


in. 


DY = 


-0.001167 


in. 


DZ = 


-0.084392 


7 = 


8.5100 


msec 


DX = 


8.888877 


in. 


DY = 


-8.881206 


in. 


DZ = 


-8.084634 
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T = 


0.5388 


msec 


ox = 


8.888873 


in. 


OY = 


-0.081233 


in. 


02 


= 


-3,884364 


T = 


8.5589 


msec 


ox = 


8.008867 


in. 


or = 


-8.381248 


in. 


02 




-8.885881 


T = 


8.5788 


msec 


ox = 


8.888868 


in. 


DY = 


-0.881251 


in. 


02 


= 


-0.885283 


T = 


8.5980 


msec 


ox = 


9.388353 


in. 


OY = 


-8.831241 


in. 


02 


= 


-3.885468 


T = 


8.6188 


msec 


ox = 


0.888846 


in. 


OY = 


-8.881219 


in, 


02 


= 


-8.885634 


T = 


8.6388 


msec 


ox = 


8.888838 


in. 


OY = 


-8.381186 


in. 


02 


= 


-0.085788 


T = 


8.6588 


msec 


ox = 


8.888838 


in. 


DY = 


-8.881144 


in. 


02 


= 


-0.885983 


T = 


8.6790 


msec 


ox = 


8.888021 


in. 


OY = 


-8.801893 


in. 


02 


= 


-8.886884 


T = 


8.6988 


msec 


ox = 


8.888814 


in. 


OY = 


-8.881835 


in. 


02 


= 


-8.00i879 


T = 


8.7188 


msec 


ox = 


8.888886 


in. 


OY = 


-8.808972 


in. 


02 


= 


-8.386129 


T = 


8.7388 


msec 


ox = 


-0.888981 


in. 


DY = 


-8.888985 


in. 


02 


= 


-8.836152 


T = 


8.7588 


msec 


ox = 


-8.388888 


in. 


OY = 


-0.880837 


in. 


02 


= 


-8.886146 


T = 


8.7788 


msec 


ox = 


-8.800014 


in. 


OY = 


-8.988778 


in. 


02 


= 


-8.886113 


T = 


8.7988 


msec 


ox = 


-8.898819 


in. 


OY = 


-8.888786 


in. 


02 


= 


-0.386858 


T = 


8.8188 


msec 


ox = 


-8.888824 


in. 


OY = 


-8.888646 


in. 


02 


= 


-0.805957 


T = 


8.8388 


msec 


ox = 


-8.888827 


in. 


OY = 


-0.888592 


in. 


02 


= 


-3.985835 


T = 


8.8588 


msec 


ox = 


-8.808038 


in. 


OY = 


-0.888546 


in. 


02 


= 


-8.885683 


T = 


8.8788 


msec 


ox = 


-8.838832 


in. 


OY = 


-8.888589 


in. 


02 


= 


-8.005502 


T = 


8.8988 


msec 


ox = 


-0.008832 


in. 


OY = 


-8.888483 


in. 


02 


= 


-8.885291 


T = 


8.9188 


msec 


ox = 


-3.388832 


in. 


OY = 


-0.888467 


in. 


02 


= 


-8.885858 


T = 


8.9388 


msec 


ox = 


-0.808031 


in. 


DY = 


-8.388462 


in. 


02 


= 


-8.884781 


T = 


8.9588 


msec 


ox = 


-8.888829 


in. 


OY = 


-8.988469 


in. 


02 


= 


-8.804484 


T = 


8.9788 


msec 


ox = 


-0.800825 


in. 


OY = 


-8.888488 


in. 


02 


= 


-0.884158 


T = 


8.9988 


msec 


ox = 


-8.808822 


in. 


or = 


-0.888516 


in. 


02 


= 


-8.883886 


T = 


1.8188 


msec 


ox = 


-8.800817 


in. 


OY = 


-8.888555 


in. 


02 


= 


-0.083426 


T = 


1.8388 


msec 


ox = 


-8.088812 


in. 


OY = 


-8.888683 


in. 


02 


= 


-8.883821 


T = 


1.8588 


msec 


ox = 


-8.888886 


in. 


OY = 


-0.880658 


in. 


02 


= 


-8.802592 


T = 


1 .0790 


msec 


ox = 


9.888983 


in. 


OY = 


-8.888719 


in. 


02 


= 


-8.882141 


T = 


1.8988 


msec 


ox = 


8.888886 


in. 


or = 


-8.888783 


in. 


02 


= 


-0.881668 


T = 


1.1188 


msec 


ox = 


8.989913 


in. 


OY = 


-8.888858 


in. 


02 


= 


-8.881175 


T = 


1.1388 


msec 


ox = 


8.808819 


in. 


DY = 


-8.888917 


in. 


02 


S 


-8.888663 


T = 


1.1588 


msec 


ox = 


8.888826' 


in. 


OY = 


-8.888982 


in. 


02 


= 


-8.830135 


T = 


1.1788 


msec 


ox = 


8.888832 


in. 


OY = 


-8.881844 


in. 


02 


= 


8.880488 


T = 


1.1988 


msec 


ox = 


8.888837 


in. 


DY = 


-8.881188 


in. 


02 


= 


8.888965 


T = 


1.2180 


msec 


ox = 


8.888843 


in. 


OY = 


-0.881149 


in. 


02 


= 


e. 881533 


T = 


1.2308 


msec 


ox = 


8.800047 


in. 


DY = 


-8.881189 


in. 


02 




8.982111 


T = 


1.2588 


msec 


ox = 


8.888851 


in. 


DY = 


-8.881219 


in. 


02 


= 


8.882697 


T = 


1.2708 


msec 


ox = 


8.838854 


in. 


OY = 


-0.881238 


in. 


02 


= 


3.883298 


T = 


1.2988 


msec 


ox = 


8.888856 


in. 


OY = 


-8.881246 


in. 


02 


= 


0.803887 


T = 


1.3180 


msec 


ox = 


3.088057 


in. 


OY = 


-8,881242 


in. 


02 


= 


8.834486 


T = 


1.3388 


msec 


ox = 


8.888857 


in. 


OY = 


-8.081226 


in. 


02 


= 


0.805885 


T = 


1.3588 


msec 


ox = 


8.998856 


in. 


DY = 


-8.881199 


in. 


02 


= 


e. 085632 


T = 


1.3708 


msec 


ox = 


8.888854 


in. 


OY = 


-8.881161 


in. 


02 


= 


8.886274 


T = 


1.3988 


msec 


ox = 


0.889951 


in. 


or = 


-8.881113 


in. 


02 


= 


8.386860 


T = 


1.4188 


msec 


ox = 


8.888847 


in. 


OY = 


-8.881857 


in. 


02 


= 


8.887437 


T = 


1.4308 


msec 


ox = 


8.880042 


in. 


OY = 


-8.888993 


in. 


02 


= 


8.988884 


T = 


1.4590 


msec 


ox = 


8.888837 


in. 


OY = 


-8.888924 


in. 


02 


= 


8.888558 


T = 


1.4788 


msec 


ox = 


9.888838 


in. 


or = 


-8.888852 


in. 


02 


= 


8.889898 


T = 


1.4980 


msec 


ox = 


8.888823 


in. 


DY = 


-8.888777 


in. 


02 


= 


8.889621 


T = 


1.5188 


msec 


ox = 


8.998815 


in. 


DY = 


-8.888783 


in. 


02 


= 


8.810127 


T = 


1.5308 


msec 


ox = 


8.088887 


in. 


OY = 


-8.888632 


in. 


02 


= 


8.018612 
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T = 


1.5580 


msec 


ox = 


-0.008881 


in. 


OY = 


-0.808584 


in. 


02 = 


3.011875 


T = 


1.5700 


resec 


ox = 


-0.000018 


in. 


OY = 


-0.900583 


in. 


02 = 


0.011514 


T = 


1.5980 


rasec 


ox = 


-0.800019 


in. 


OY = 


-0.090449 


in. 


02 = 


8.011927 


T = 


1.8100 


resec 


ox = 


-0.088028 


in. 


OY = 


-0.809483 


in. 


02 = 


0.012312 


T = 


1.8300 


resec 


ox = 


-0.008838 


in. 


OY = 


-0.000388 


in. 


02 = 


0.012887 


T = 


1.8500 


msec 


ox = 


-0.308045 


in. 


OY = 


-0.030344 


in. 


02 = 


0.812990 


T _ 

i - 


1.8780 


msec 


ox = 


-0,900853 


in. 


DY = 


-0.000331 


in. 


02 = 


9.013282 


T = 


1.8900 


resec 


ox = 


-8.000080 


in. 


OY = 


-9.800329 


in. 


02 = 


0.013540 


T = 


1.7180 


rasec 


ox = 


-8.008987 


in. 


OY = 


-0.003339 


in. 


02 = 


0.013784 


T = 


1.7308 


resec 


ox = 


-0.000872 


in. 


OY = 


-0.900359 


in. 


02 = 


0.813952 


T = 


1.7580 


resec 


ox = 


-0.000077 


in. 


DY = 


-0.000390 


in. 


02 = 


8.014104 


T = 


1.7708 


resec 


ox = 


-8.000082 


in. 


OY = 


-0.808430 


in. 


02 = 


8.814219 


T = 


1 .7900 


resec 


ox = 


-0.000885 


in. 


OY = 


-0.008477 


in. 


02 = 


0.014297 


T = 


1.3188 


resec 


ox = 


-0.000087 


in. 


DY = 


-9.808531 


in. 


02 = 


0.914338 


T = 


1.8380 


resec 


ox = 


-0.808838 


in. 


OY = 


-0.038588 


in. 


02 = 


8.814337 


T = 


1.8508 


resec 


ox = 


-0.000088 


in. 


OY = 


-0.000848 


in. 


02 = 


0.014299 


T = 


1.8708 


rasec 


ox = 


-0.000887 


in. 


DY = 


-0 .000709 


in. 


02 = 


0.014222 


T = 


1.3900 


resec 


ox = 


-0.000885 


in. 


OY = 


-8,880789 


in. 


02 = 


0.014135 


T = 


1.9180 


rasec 


ox = 


-8.000083 


in. 


DY = 


-0.000824 


in. 


02 = 


0.013950 


T = 


1.9300 


resec 


ox = 


-0.888079 


in. 


DY = 


-0.989875 


in. 


02 = 


3.013755 


T = 


1.9500 


rasec 


ox = 


-0.000075 


in. 


DY = 


-0.000919 


in. 


D2 = 


3.813522 


T = 


1.9700 


resec 


ox = 


-0.000870 


in. 


DY = 


-0.889954 


in. 


02 = 


0.013251 


T = 


1.9900 


resec 


ox = 


-0.000084 


in. 


OY = 


-0 .008980 


in. 


02 = 


0.012943 


T = 


2.0100 


resec 


ox = 


-9.000058 


in. 


DY = 


-0.300995 


in. 


D2 = 


9.012598 


T = 


2.0300 


rasec 


ox = 


-0.000051 


in. 


DY = 


-0.000999 


in, 


02 = 


0.012218 


T = 


2.0580 


rasec 


ox = 


-8.800845 


in. 


OY = 


-8.800991 


in. 


02 = 


0.311304 


7 = 


2.0700 


resec 


ox = 


-9.000038 


in. 


OY = 


-0.090971 


in. 


02 = 


0.011355 


T = 


2.0980 


resec 


ox = 


-0.808032 


in. 


DY = 


-0.080948 


in. 


02 = 


0.913875 


T = 


2.1100 


rasec 


ox = 


-0.000025 


in. 


OY = 


-0.000898 


in. 


02 = 


8.018383 


T = 


2.1380 


resec 


DX = 


-0.080019 


in. 


DY = 


-8.900848 


in. 


02 = 


0.009822 


T = 


2.1500 


rasec 


ox = 


-0.080014 


in. 


DY = 


-0.000788 


in. 


02 = 


0.089252 


T = 


2.1700 


resec 


DX = 


-0.800009 


in. 


OY = 


-0.008718 


in. 


02 = 


0.808854 


T = 


2.1900 


resec 


ox = 


-0.800885 


in. 


OY = 


-0.000845 


in. 


02 = 


8.898033 


T = 


2.2100 


resec 


DX = 


-9.909991 


in. 


OY = 


-0.800588 


in. 


02 = 


0.807383 


T = 


2.2300 


rasec 


ox = 


9.909091 


in. 


OY = 


-0.000489 


in. 


02 = 


8.008713 


T = 


2.2508 


rasec 


ox = 


9.990903 


in. 


DY = 


-8.000410 


in. 


02 = 


9.998824 


T = 


2.2780 


rasec 


ox = 


0.999994 


in. 


OY = 


-0.000333 


in. 


02 = 


0.805318 


T = 


2.2900 


resec 


ox = 


9.800003 


in. 


OY = 


-0.008280 


in. 


02 = 


0.804595 


T = 


2.3100 


rasec 


ox = 


0.909992 


in. 


OY = 


-0.000193 


in. 


02 = 


8.003858 


T = 


2.3300 


rasec 


DX = 


-9.009999 


in. 


OY = 


-0.080133 


in. 


02 = 


9.993111 


T = 


2.3588 


rasec 


ox = 


-0.000004 


in. 


OY = 


-0.088882 


in. 


02 = 


8.802354 


T = 


2.3700 


rasec 


DX = 


-0.000008 


in. 


OY = 


-8.000941 


in. 


02 = 


0.081589 


T = 


2.3900 


rasec 


ox = 


-0.000013 


in. 


OY = 


-0.000011 


in. 


02 = 


0.900820 


T = 


2.4108 


resec 


DX = 


-0.800019 


in. 


DY = 


0.808007 


in. 


02 = 


9.000048 


T = 


2.4300 


rasec 


ox = 


-0.000028 


in. 


OY = 


0.900013 


in. 


02 = 


-8.000724 


T = 


2.4580 


rasec 


DX = 


-8.800833 


in. 


DY = 


3.999098 


in. 


02 = 


-8.881493 


T = 


2.4700 


rasec 


DX = 


-8.800841 


in. 


OY = 


-0.008888 


in. 


02 = 


-8.802258 


T = 


2.4908 


resec 


ox = 


-8.800049 


in. 


OY = 


-8.008835 


in. 


02 = 


-0.003817 


T = 


2.5100 


rasec 


ox = 


-0.000057 


in. 


OY = 


-0.000071 


in. 


02 = 


-0.033785 


T = 


2.5380 


rasec 


ox = 


-8.080885 


in. 


DY = 


-0,000118 


in. 


02 = 


-0.304582 


T = 


2.5500 


rasec 


ox = 


-8.800074 


in. 


OY = 


-0.000187 


in. 


02 = 


-0.885228 
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T = 


2.5700 


msec 


DX = 


-0.098082 


in. 


DY = 


-8.000223 


in. 


DZ 


= 


-8.005933 


T = 


2.5908 


msec 


OX = 


-0.809090 


in. 


DY = 


-8.088282 


in. 


DZ 


= 


-0.086623 


T = 


2.6198 


msec 


DX = 


-8.989098 


in. 


DY = 


-0.089342 


in, 


DZ 


= 


-8.807294 


T = 


2.6308 


msec 


DX = 


-9.890185 


in. 


DY = 


-8.800402 


in. 


DZ 


= 


-9.037942 


T = 


2.6509 


msec 


OX = 


-8.080111 


in. 


DY = 


-0.008453 


in, 


DZ 


- 


-0.908566 


T = 


2.6700 


msec 


OX = 


-0.089117 


in. 


DY = 


-0.009510 


in. 


DZ 


r 


-0.089166 


T = 


2.6980 


msec 


DX = 


-0.000122 


in. 


DY = 


-0.008555 


in. 


DZ 


= 


-8.80973? 


T = 


2.7100 


msec 


OX = 


-0.800126 


in. 


DY = 


-0.009592 


in. 


DZ 


= 


-8.010281 


T = 


2.7300 


msec 


DX = 


-0.90012? 


in. 


DY = 


-0.008628 


in. 


DZ 


= 


-0.010794 


T = 


2.7500 


msec 


OX = 


-0.000131 


in. 


DY = 


-8.000633 


in. 


DZ 


= 


-9.811274 


T = 


2.7700 


msec 


OX = 


-0.000132 


in. 


DY = 


-0.000644 


in. 


DZ 


= 


-0.011720 


T = 


2.7909 


msec 


DX = 


-9.809131 


in. 


DY = 


-0.00063? 


in, 


DZ 


= 


-0.012133 


T = 


2.8100 


msec 


DX = 


-0.000130 


in. 


DY = 


-0.000623 


in. 


DZ 


= 


-0.012510 


T = 


2.8300 


msec 


DX = 


-0.000128 


in. 


DY = 


-8.088595 


in. 


DZ 


= 


-0.012350 


T = 


2.3500 


msec 


DX = 


-8.000125 


in. 


DY = 


-8.000556 


in. 


DZ 


= 


-0.013155 


T = 


2.8700 


msec 


DX = 


-0.000121 


in. 


DY = 


-0.000507 


in. 


DZ 


r 


-8.013421 


T = 


2.8908 


msec 


DX = 


-8.000116 


in. 


DY = 


-0.00044? 


in, 


DZ 


= 


-0.813648 


T _ 
1 - 


2.9108 


msec 


DX = 


-0.800111 


in. 


DY = 


-0.000384 


in. 


DZ 


= 


-0.813837 


T = 


2.9380 


msec 


DX = 


-8.888105 


in. 


DY = 


-8.088313 


in. 


DZ 


= 


-0.013987 


T = 


2.9500 


msec 


DX = 


-0.000098 


in. 


DY = 


-0.000233 


in. 


DZ 


= 


-0.014997 


T = 


2.9700 


msec 


DX = 


-0.800091 


in. 


DY = 


-0.800161 


in. 


DZ 


= 


-0.014166 


T = 


2.9988 


msec 


DX = 


-9.080084 


in. 


DY = 


-8.888885 


in. 


DZ 


= 


-8.914196 


T = 


3.0100 


msec 


DX = 


-9.000977 


in. 


DY = 


-0.000010 


in. 


DZ 


= 


-0.814135 


T = 


3.3300 


msec 


DX = 


-9.099070 


in. 


DY = 


0.008061 


in. 


DZ 


= 


-9.014135 


T = 


3.8580 


msec 


DX = 


-0.000063 


in. 


DY = 


0.000126 


in. 


DZ 




-0.814046 


T = 


3.8780 


msec 


DX = 


-0.989956 


in. 


DY = 


8.800184 


in. 


DZ 




-9.913917 


T = 


3.0980 


msec 


DX = 


-0.080050 


in. 


DY = 


0.880233 


in. 


DZ 


= 


-0.813751 


T = 


3.1190 


msec 


DX = 


-0.800944 


in. 


DY = 


0.000272 


in. 


DZ 


= 


-8.013548 


T = 


3.1300 


msec 


DX = 


-0.00003? 


in. 


DY = 


0.080299 


in. 


DZ 




-8.913309 


T = 


3.1500 


msec 


DX = 


-0.090035 


in. 


DY = 


0.099315 


in. 


DZ 


= 


-8.813034 


T = 


3.1788 


msec 


DX = 


-0.080831 


in. 


DY = 


8.888320 


in. 


DZ 


= 


-0.012726 


T = 


3.1990 


msec 


DX = 


-8.003028 


in. 


DY = 


8.009312 


in. 


DZ 




-9.012385 


T = 


3.2100 


msec 


DX = 


-0.900027 


in. 


DY = 


0.000294 


in. 


DZ 


= 


-8.01201! 


T = 


3.2388 


msec 


DX = 


-9.900026 


in. 


DY = 


8.088264 


in. 


DZ 


= 


-9.011688 


T = 


3.2508 


msec 


DX = 


-0.000026 


in. 


DY = 


0.080225 


in. 


DZ 


r 


-0.811176 


T = 


3.2709 


msec 


DX = 


-0.009027 


in. 


DY = 


0.009177 


in. 


DZ 


= 


-8.019715 


T = 


3.2990 


msec 


DX = 


-0.000030 


in. 


DY = 


0.000123 


in. 


DZ 




-8.010231 


T = 


3.3199 


msec 


DX = 


-8.099033 


in. 


DY = 


9.080963 


in. 


DZ 


= 


-9.989723 


T = 


3.3389 


msec 


DX = 


-0.000037 


in. 


DY = 


0.880000 


in. 


DZ 


= 


-0.089192 


T = 


3.3590 


msec 


DX = 


-0.000042 


in. 


DY = 


-0.890064 


in. 


DZ 


= 


-0.998642 


T = 


3.3798 


msec 


DX = 


-0.908047 


in. 


DY = 


-0.000127 


in. 


DZ 


= 


-0.808076 


T _ 
1 - 


3.3980 


msec 


DX = 


-0.090054 


in. 


DY = 


-0.000138 


in. 


DZ 


= 


-0.087492 


T = 


3.4100 


msec 


DX = 


-9.998869 


in. 


DY = 


-8.000244 


in. 


DZ 


= 


-8.886897 


T = 


3.4388 


msec 


DX = 


-0.000068 


in. 


DY = 


-0.908294 


in. 


DZ 


= 


-0.086299 


T = 


3.4599 


msec 


DX = 


-0.000075 


in. 


DY = 


-0.088336 


in. 


DZ 


= 


-0.005673 


T = 


3.4700 


msec 


DX = 


-0.090983 


in. 


DY = 


-0.909368 


in. 


DZ 


= 


-0.905050 


T = 


3.4999 


msec 


DX = 


-9.900991 


in. 


DY = 


-8.080398 


in. 


DZ 


r 


-8.004422 


T = 


3.5109 


msec 


DX = 


-0.000098 


in. 


DY = 


-0.000409 


in. 


DZ 


= 


-0.003791 


T = 


3.5300 


msec 


DX = 


-8.099195 


in. 


DY = 


-8.08839? 


in. 


DZ 


= 


-0.003159 


T = 


3.5500 


msec 


DX = 


-0.080112 


in. 


DY = 


-0.009387 


in. 


DZ 


= 


-0.092539 


T = 


3.5700 


msec 


DX = 


-0.00911? 


in. 


DY = 


-9.080362 


in. 


DZ 


= 


-9.801905 
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T = 


3.5900 


msec 


DX 


= 


-0.080125 


in. 


DY = 


-0.308327 


in. 


DZ = 


-0.081286 


T = 


3.5109 


msec 


OX 


= 


-3.003133 


in. 


DY = 


-8.033281 


in. 


DZ = 


-0.808677 


T = 


3.6300 


msec 


OX 


= 


-0.000134 


in. 


DY = 


-8.330227 


in. 


DZ = 


-0.000079 


T = 


3.6533 


msec 


DX 




-0.300138 


in. 


DY = 


-8.308165 


in, 


DZ = 


0.380583 


T = 


3.6700 


msec 


OX 




-0.088140 


in. 


DY = 


-0.000097 


in. 


DZ = 


8.081878 


T = 


3.6900 


msec 


ox 




-3.008142 


in. 


DY = 


-0.088024 


in. 


DZ = 


8.001632 


T = 


3.7100 


msec 


DX 


= 


-8.088142 


in. 


DY = 


8.800850 


in. 


DZ = 


8.882167 


T = 


3.7300 


msec 


DX 


= 


-8.888142 


in. 


DY = 


8.888125 


in. 


DZ = 


0.802681 


T = 


3.7500 


msec 


DX 


= 


-0 .088148 


in. 


DY = 


0.880197 


in. 


DZ = 


0.883174 


T = 


3.7700 


msec 


DX 


- 


-8.888137 


in. 


DY = 


8.338266 


in. 


DZ = 


0.333643 


T = 


3.7900 


msec 


DX 


= 


-8.888134 


in. 


DY = 


8.888338 


in. 


DZ = 


0.804886 


T = 


3.8130 


msec 


DX 


= 


-8.888138 


in. 


DY = 


8.888387 


in. 


DZ = 


3.884501 


T = 


3.8300 


msec 


DX 


= 


-0.000124 


in. 


DY = 


8.000434 


in. 


DZ = 


0.094889 


T = 


3.3500 


msec 


DX 


= 


-8.388118 


in. 


DY = 


8.838472 


in. 


DZ = 


8.885246 


T = 


3.8700 


msec 


DX 


= 


-0.080112 


in. 


DY = 


8.080499 


in. 


DZ = 


8.885571 


T = 


3.8900 


msec 


DX 


= 


-8.880185 


in. 


DY = 


8.888514 


in, 


DZ = 


0.805865 


T = 


3.9100 


msec 


DX 


= 


-0.800897 


in. 


DY = 


8.800518 


in. 


DZ = 


8.886125 


T = 


3.9300 


msec 


DX 


= 


-8.388889 


in. 


DY = 


8.808589 


in. 


DZ = 


0.386351 


T = 


3.9500 


msec 


DX 


= 


-8.883881 


in. 


DY = 


8.888489 


in. 


DZ = 


8.806543 


T = 


3.9703 


msec 


DX 


= 


-8.888873 


in. 


DY = 


8.888458 


in. 


DZ = 


9.836699 


T = 


3.9900 


msec 


DX 


= 


-0.000065 


in. 


DY = 


0.000418 


in. 


DZ = 


0.006819 


T = 


4.0100 


msec 


DX 


= 


-8.008858 


in. 


DY = 


3.388369 


in. 


DZ = 


0 .006985 


T = 


4.0300 


msec 


DX 


= 


-0.000051 


in. 


DY = 


8.000313 


in. 


DZ = 


8.086953 


T = 


4.0500 


msec 


DX 


= 


-8.888844 


in. 


DY = 


3.888252 


in. 


DZ = 


8.886966 


T = 


4.0700 


msec 


DX 


= 


-0.000838 


in. 


DY = 


0.008187 


in. 


DZ = 


0.886944 


T = 


4.0900 


msec 


DX 


= 


-0.388832 


in. 


DY = 


0.088121 


in. 


DZ = 


8.886384 


T = 


4.1130 


msec 


DX 


= 


-0.000028 


in. 


DY = 


8.000056 


in. 


DZ = 


8.806789 


T = 


4.1300 


msec 


DX 


= 


-0.380024 


in. 


DY = 


-8.800087 


in. 


DZ = 


8.386659 


T = 


4.1500 


msec 


DX 


= 


-8.800821 


in. 


DY = 


-0.080065 


in. 


DZ = 


8.006494 


T = 


4.1730 


msec 


DX 


= 


-8.880020 


in. 


DY = 


-0.888113 


in. 


DZ = 


8.806294 


T = 


4.1900 


msec 


DX 


= 


-0.000019 


in. 


DY = 


-0.000162 


in. 


DZ = 


8.006062 


T = 


4.2100 


msec 


DX 


= 


-0.000019 


in. 


DY = 


-0.880197 


in. 


DZ = 


8.005796 


T = 


4.2300 


msec 


DX 


= 


-0.080020 


in. 


DY = 


-8.000222 


in. 


DZ = 


0.005500 


T = 


4.2500 


msec 


DX 


= 


-0.008822 


in. 


DY = 


-8.008236 


in. 


DZ = 


8.885174 


T = 


4.2700 


msec 


DX 


= 


-0.000025 


in. 


DY = 


-0.000238 


in. 


DZ = 


8.004818 


T = 


4.2900 


msec 


DX 


= 


-8.888829 


in. 


DY = 


-0.000228 


in. 


DZ = 


0.804436 


T — 
1 - 


4.3100 


msec 


DX 


= 


-0.800034 


in. 


DY = 


-0.808207 


in. 


DZ = 


0.084027 


T = ' 


4.3300 


msec 


DX 


= 


-0.000039 


in. 


DY = 


-0.338175 


in. 


DZ = 


0.883593 


T = 


4.3500 


msec 


DX 


= 


-8 .880045 


in. 


DY = 


-8.000132 


in. 


DZ = 


0.003135 


T = 


4.3700 


msec 


DX 


= 


-8.000052 


in. 


DY = 


-0.00008! 


in. 


DZ = 


0.002657 


T = 


4.3900 


msec 


DX 


= 


-8.000058 


in. 


DY = 


-0.000021 


in. 


DZ = 


8.002158 


T = 


4.4100 


msec 


DX 




-8.888865 


in. 


DY = 


8.808044 


in. 


DZ = 


0.001641 


T = 


4.4300 


msec 


DX 


= 


-0.000072 


in. 


DY = 


0.000114 


in. 


DZ = 


0.001109 


T = 


4.4500 


msec 


DX 


= 


-8.838879 


in. 


DY = 


0.000187 


in. 


DZ = 


8.008562 


T = 


4.4700 


msec 


DX 


= 


-0 .000886 


in. 


DY = 


0.000260 


in. 


DZ = 


0.090805 


T = 


4.4903 


msec 


DX 


= 


-8.800093 


in. 


DY = 


0.000331 


in. 


DZ = 


-8.000563 


T = 


4.5100 


msec 


DX 


= 


-8.880099 


in. 


DY = 


0.000399 


in. 


DZ = 


-0.081138 


T = 


4.5300 


msec 


DX 


= 


-0 .880104 


in. 


DY = 


0.000462 


in. 


DZ = 


-8.001719 


T = 


4.5500 


msec 


DX 


= 


-0.000109 


in. 


DY = 


0.000517 


in. 


DZ = 


-0.302303 


T = 


4.5700 


msec 


DX 


= 


-8.888113 


in. 


DY = 


0.000565 


in. 


DZ = 


-0.092888 


T = 


4.5900 


msec 


DX 


= 


-0.000116 


in. 


DY = 


0.000602 


in. 


DZ = 


-8.883472 



110 



T = 


4.6188 


resec 


T = 


4.6388 


resec 


T = 


4.6500 


resec 


T = 


4.6788 


resec 


T = 


4.6900 


resec 


T = 


4.7188 


resec 


T = 


4.7308 


resec 


T = 


4.7588 


resec 


T = 


4.7700 


resec 


T = 


4.7900 


resec 


T = 


4.8100 


resec 


T = 


4.8388 


resec 


T = 


4.8500 


resec 


T = 


4.8790 


resec 


T = 


4.8908 


resec 


T = 


4.9180 


resec 


T = 


4.9388 


resec 


T = 


4.9588 


resec 


T = 


4.9708 


resec 


T = 


4.9988 


resec 


T = 


5.0100 


resec 


T = 


5.0300 


resec 


T = 


5.0508 


resec 


T = 


5.8788 


resec 


T = 


5.0988 


resec 


T = 


5.1108 


resec 


T = 


5.1380 


resec 


T = 


5.1588 


resec 


T = 


5.1708 


resec 


T = 


5.1988 


resec 


T = 


5.2180 


resec 


T = 


5.2300 


resec 


T = 


5.2508 


resec 


T — 
1 - 


5.2708 


resec 


T = 


5.2988 


msec 


T = 


5.3188 


resec 


T = 


5.3388 


resec 


T = 


5.3588 


resec 


T = 


5.3780 


resec 


T = 


5.3980 


resec 


T = 


5.4188 


resec 


T = 


5.4388 


resec 


T = 


5.4580 


resec 


T = 


5.4788 


resec 


T = 


5.4900 


msec 


T = 


5.5108 


resec 


T = 


5.5300 


resec 


T = 


5.5508 


resec 


T = 


5.5700 


resec 


T = 


5.5988 


resec 


T = 


5.6108 


resec 



DX = 


-B .880113 


in. 


DX = 


-0.088119 


in. 


DX = 


-8.000120 


in. 


DX = 


-0.880119 


in. 


DX = 


-0.00011? 


in. 


DX = 


-0.888114 


in. 


DX = 


-0.900110 


in. 


DX = 


-8.809185 


in. 


DX = 


-8.800100 


in. 


DX = 


-8.880893 


in. 


DX = 


-8.080086 


in. 


DX = 


-8.88807? 


in. 


DX = 


-0.000071 


in. 


DX = 


-0.800862 


in. 


DX = 


-0.000054 


in. 


DX = 


-8.808845 


in. 


DX = 


-8.000037 


in. 


DX = 


-0.900028 


in. 


DX = 


-0.800020 


in. 


DX = 


-0.880813 


in, 


DX = 


-0.809006 


in. 


DX = 


0.000000 


in. 


DX = 


0.000006 


in. 


DX = 


8.880811 


in, 


DX = 


0.000814 


in, 


DX = 


8.908817 


in. 


DX = 


0.000019 


in. 


DX = 


8.00881? 


in, 


DX = 


0.80001? 


in. 


DX = 


0.888018 


in. 


DX = 


8.008016 


in, 


DX = 


0.000013 


in. 


DX = 


8.000809 


in. 


DX = 


8.000084 


in. 


DX = 


-0.800801 


in, 


DX = 


-8.009097 


in. 


DX = 


-8.888813 


in. 


DX = 


-8.808819 


in. 


DX = 


-8.880826 


in. 


DX = 


-8.880832 


in. 


DX = 


-0.888838 


in. 


DX = 


-0.888844 


in. 


DX = 


-0.000850 


in, 


DX = 


-8.888055 


in. 


DX = 


-8 .808059 


in. 


DX = 


-8.888063 


in. 


DX = 


-8.880866 


in. 


DX = 


-0.880868 


in. 


DX = 


-0.000069 


in. 


DX = 


-8.000869 


in. 


DX = 


-8.008068 


in. 



DY = 


8.80062? 


in. 


DY = 


0.088644 


in. 


DY = 


0.888643 


in. 


DY = 


0.000640 


in, 


DY = 


0.800620 


in. 


DY = 


0.800590 


in. 


DY = 


8.008558 


in. 


DY = 


8.000591 


in. 


DY = 


8.000445 


in. 


DY = 


0.088383 


in. 


DY = 


0.000318 


in, 


DY = 


8.880251 


in. 


DY = 


0.000184 


in. 


DY = 


0.880120 


in. 


DY = 


0.000060 


in, 


DY = 


0.008005 


in. 


DY = 


-0.000042 


in. 


DY = 


-8.808879 


in. 


DY = 


-0.08010? 


in. 


DY = 


-8.800124 


in. 


DY = 


-0.00012? 


in . 


DY = 


-0.000123 


in. 


DY = 


-0.000105 


in. 


DY = 


-8.888876 


in, 


DY = 


-0.009037 


in . 


DY = 


0.000812 


in. 


DY = 


0.001498 


in. 


DY = 


8.88156! 


in. 


DY = 


0.091629 


in. 


DY = 


8.001709 


in. 


DY = 


•0.801771 


in. 


DY = 


8.881842 


in. 


DY = 


0.00190? 


in. 


DY = 


0.001971 


in. 


DY = 


8.002027 


in. 


DY = 


0.892974 


in. 


DY = 


8.002112 


in. 


DY = 


0.88213? 


in. 


DY = 


8.002155 


in. 


DY = 


0.002168 


in. 


DY = 


0.092152 


in . 


DY = 


0.002133 


in. 


DY = 


8.002104 


in. 


DY = 


8.882864 


in. 


DY = 


8.082015 


in. 


DY = 


8.001959 


in. 


DY = 


8.801898 


in. 


DY = 


8.801832 


in. 


DY = 


8.001764 


in. 


DY = 


0.001696 


in. 


DY = 


0.001638 


in. 



Ill 



DZ = -0.804853 
DZ = -0.89462? 
OZ = -0.365196 
OZ = -0.885754 
DZ = -9.08629? 
DZ = -8.88683! 
DZ = -0.897345 
DZ = -8.337341 
DZ = -0.008318 
DZ = -0.088772 
DZ = -8.009282 
DZ = -0.009693 
DZ = -8.009987 
DZ = -8.018337 
DZ = -0.018658 
DZ = -0.813943 
DZ = -0.811285 
DZ = -0.81142? 
DZ = -0.011618 
DZ = -0.811771 
DZ = -0.011883 
DZ = -8.81196? 
DZ = -0.012812 
DZ = -8.812017 
DZ = -0.011985 
DZ = -8.911914 
DZ = -8.811806 
DZ = -0.811659 
DZ = -8.011474 
DZ = -0.811252 
DZ = -0.810992 
DZ = -8.810696 
DZ = -8.818363 
DZ = -0.009996 
DZ = -8.089594 
DZ = -0.089160 
DZ = -8.088694 
DZ = -8.888198 
DZ = -8.087674 
DZ = -0.807122 
DZ = -0.086544 
DZ = -0.805943 
DZ = -0.085318 
DZ = -8.884674 
DZ = -0.084018 
DZ = -9.00332? 
DZ = -0.092634 
DZ = -8.881927 
DZ = -8.801209 
OZ = -0.000483 
DZ = 0.888248 



T = 


5.(5300 


msec 


DX 


= 


-0.000066 


in. 


DY = 


8.831567 


in. 


D2 


- 


8.800933 


T = 


5.6500 


msec 


DX 


= 


-0.003063 


in. 


DY = 


0.801510 


in. 


DZ 


= 


0.301728 


T = 


5.6700 


msec 


DX 


= 


-0.888359 


in. 


II 


8.801461 


in. 


DZ 


= 


8.002456 


T = 


5.6908 


msec 


DX 


= 


-0.000854 


in. 


DY = 


8.081420 


in. 


DZ 


= 


0.083189 


T = 


5.7130 


msec 


DX 


= 


-8.888849 


in. 


DY = 


8.881389 


in. 


DZ 


= 


8.803917 


T = 


5.7300 


msec 


DX 


= 


-0.333042 


in. 


DY = 


8.001369 


in. 


DZ 




0.034638 


T = 


5.7500 


msec 


DX 


= 


-8.080035 


in. 


DY = 


0.001360 


in. 


DZ 


= 


0.805358 


T = 


5.7708 


msec 


DX 


= 


-0.000827 


in. 


S 

ii 


0.801363 


in. 


DZ 


S 


0.806051 


T = 


5.7988 


msec 


DX 


= 


-8.088019 


in. 


DY = 


8.881377 


in. 


DZ 


= 


0.886737 


T = 


5.8180 


msec 


DX 


= 


-0.800818 


in. 


DY = 


0.001403 


in. 


DZ 


= 


0.037408 


T = 


5.8308 


msec 


DX 


= 


-8.000001 


in. 


DY = 


0.001439 


in. 


DZ 




0.308861 


T = 


5.8500 


msec 


DX 


= 


0.000008 


in. 


DY = 


0.381485 


In. 


DZ 


= 


9.008694 


T = 


5.8788 


msec 


DX 


3 


8.888017 


in. 


DY = 


8.301539 


in. 


DZ 


= 


8.889305 


T = 


5.8908 


msec 


DX 


= 


8.000025 


in. 


DY = 


8.881599 


in. 


DZ 


= 


3.009892 


T = 


5.9188 


msec 


DX 


= 


8.000834 


in. 


DY = 


0.001665 


in. 


DZ 


= 


8.018454 


T = 


5.9300 


msec 


DX 


= 


8.003042 


in. 


DY = 


8.081734 


in. 


DZ 


= 


9.819998 


T = 


5.9588 


msec 


DX 


= 


0.830049 


in. 


DY = 


0.081804 


in. 


DZ 


= 


8.011497 


T = 


5.9780 


msec 


DX 


= 


3.000056 


in. 


DY = 


0.001873 


in. 


DZ 


= 


0.011974 


T = 


5.9988 


msec 


DX 


= 


0.000362 


in. 


DY = 


0.081939 


in. 


DZ 


r 


8.012420 


TMAX(X) 


' = 3.7100 


nsec 


DtttX(X) = -1 


3.088142 


in. 












TMAX(Y) 


= 5.3988 


nsec 


DMAX(Y) = 3.382168 


in. 












TMX<Z) 


> = 1.8303 


nsec 


DMAX(Z) = I 


3.814337 


in. 
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